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Abstract 

This  research  examines  the  performance  of  a  spectrally  encoded,  multi-carrier, 
phase  shift  keying  communications  system  in  a  frequency-selective,  slowly-fading  mul¬ 
tipath  channel.  The  specific  communications  system  modeled  is  the  transform  domain 
communication  system  (TDCS)  originally  researched  as  an  interference  avoidance 
technique.  Previous  TDCS  research  assumed  an  additive  white  Gaussian  noise  chan¬ 
nel,  which  is  not  representative  of  a  realistic  environment. 

This  thesis  presents  overviews  of  previous  TDCS  research,  the  multipath  fading 
channel,  and  the  RAKE  receiver.  Analysis  and  Matlab®  simulations  compare  the 
performance  of  spectrally  encoded  and  un-encoded  signals  through  a  multipath  fading 
channel  using  an  L-diversity  TDCS  RAKE  receiver.  Encoded  signals  take  on  the 
spectral  shape  of  the  multipath  fading  channel  transfer  function.  Un-encoded  signals 
have  a  flat  magnitude  spectrum.  The  research  also  evaluates  the  interference  rejection 
capability  of  spectrally  encoded  signals  in  a  multipath  channel. 

Research  results  indicate  for  diversities  (L)  ranging  between  2  and  50,  spectrally 
encoded  signals  need  1.0  to  2.75dB  less  transmitted  normalized  bit  energy  to  noise 
power  spectral  density  ratios  ( E^/N0 )  to  achieve  the  same  probability  of  bit  error 
(Pb)  as  un-encoded  signals.  Results  also  demonstrate  that  spectrally  encoded  TDCS 
signals  retain  the  original  interference  rejection  capability. 
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Performance  of  a  Spectrally  Encoded  Multi-Carrier 
Phase  Shift  Keying  Communication  System  in  a 
Frequency-Selective,  Slowly-Fading  Multipath  Channel 

I.  Introduction 

1.1  Background 

As  early  as  the  5th  century  B.C.,  Sun  Tzu  extolled  the  principle  of  knowing  the 
enemy  in  combat  in  the  following  proverb: 

If  you  know  the  enemy  and  know  yourself,  you  need  not  fear  the  result 
of  a  hundred  battles.  If  you  know  yourself  but  not  the  enemy,  for  every 
victory  gained  you  will  also  suffer  a  defeat.  If  you  know  neither  the  enemy 
nor  yourself,  you  will  succumb  in  every  battle  [2], 

While  the  application  has  changed  after  2,500  years  of  history,  this  principle 
still  holds  true.  Today’s  warrior  desires  an  integrated  picture  of  the  battlespace  to 
1)  reveal  the  enemy’s  every  move  and  intention,  and  2)  understand  the  status  and 
capabilities  of  friendly  forces.  The  fundamental  assumption  behind  the  existence  of 
this  integrated  picture  is  the  ubiquitous  network.  A  global  network  that  connects 
air,  space,  land  and  sea  assets  around  the  world  is  the  foundation  required  to  deliver 
actionable  information  to  the  commander  on  the  battlefield. 

The  need  for  a  pervasive  network  is  expressed  clearly  in  the  vision  documents  for 
all  of  the  United  States  armed  services.  In  the  US  Army,  it  is  known  as  LandWarNet. 
It  is  the  Navy  and  Marine  Corps  FORCEnet.  It  is  encompassed  in  the  Air  Force’s  core 
competency  of  Information  Superiority.  The  worldwide  defense  network  is  evolving 
toward  the  vision  of  the  Global  Information  Grid  (GIG)  which  consists  primarily  of 
fiber,  wire,  and  satellite  communications.  However,  near  the  forward  edge  of  the 
battle  area,  at  the  edge  of  the  network,  the  electronic  lifeline  from  the  troops  to  the 
GIG  travels  over  wireless  communications  systems.  This  lifeline  brings  situational 
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awareness,  command  and  control  and  intelligence  information  to  commanders  at  the 
very  tip  of  the  spear. 

The  wireless  communications  system  is  often  called  upon  to  provide  connections 
in  the  harshest  of  environments.  In  addition  to  the  natural  communications  barriers 
of  weather,  atmospheric  effects  and  terrain  blockages,  the  combat  environment  brings 
its  own  unique  difficulties.  The  adversary  brings  electronic  warfare  efforts  to  bear  on 
the  already  crowded  electromagnetic  spectrum  and  friendly  interference  sources.  As  if 
these  were  not  enough,  consider  that  the  battlefield  of  today’s  soldiers,  sailors,  airmen 
and  marines  is  the  city  block.  In  Bosnia,  Somalia,  and  now  in  Iraq,  urban  warfare 
has  become  the  adversary’s  answer  to  the  asymmetrical  size  and  power  of  the  US 
military.  The  urban  environment  brings  its  own  adverse  effects  on  wireless  communi¬ 
cations.  The  US  Army  Field  Manual  3-06.11,  “COMBINED  ARMS  OPERATIONS 
IN  URBAN  TERRAIN”  in  addition  to  other  topics,  dedicates  an  appendix  to  advis¬ 
ing  soldiers  on  how  to  provide  communications  capability  in  the  “complexities”  of  the 
urban  environment  [1]. 

The  urban  environment  is  not  a  new  challenge  to  the  communications  world. 
In  the  civilian  sector,  the  desire  for  convenience  has  birthed  numerous  wireless  con¬ 
nectivity  technologies.  Cellular  phones  have  provided  voice  and  data  connectivity  to 
customers  on  the  move  for  decades.  Wireless  networks  are  now  beginning  to  cover  the 
urban  landscape,  providing  personal  digital  assistants  and  laptop  computers  internet 
connectivity  at  home,  work,  restaurants  and  coffee  shops.  Connectivity  and  conve¬ 
nience  are  not  without  challenge.  Obstructions,  reflections  and  interference  have  been 
the  bane  of  reliable  service,  as  seen  in  a  family  of  cellular  phone  service  advertise¬ 
ments.  In  the  commercials,  a  well-dressed  man  roams  the  nation  with  a  cell  phone  to 
his  ear  and  asks  “Can  you  hear  me  now?” 

Whether  in  the  military  or  the  civilian  sector,  overcoming  the  reality  of  the 
harsh  communications  environment  is  topic  of  a  diverse  body  of  research.  Part  of 
that  body  is  the  application  of  spectral  encoding  techniques.  The  transform  domain 
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communication  system  (TDCS)  has  historically  applied  spectral  encoding  techniques 
to  the  problem  of  interference  avoidance. 

The  history  of  TDCS  research  by  that  name  reaches  at  least  back  to  the  mid- 
1990s  when  Radcliffe  published  his  thesis  entitled,  “Design  and  Simulation  of  a 
Transform  Domain  Communications  System”  [9].  His  research  spurred  a  body  of 
work  at  AF1T  exploring  the  application  of  spectral  encoding  to  interference  suppres¬ 
sion.  Swackhammer  researched  and  simulated  a  multiple  access  scheme  using  pseudo¬ 
random  (PR)  phase  encoding  and  cyclic  shift  keying  (CSK)  for  modulation  [13]. 
Roberts  researched  synchronization  issues  with  TDCS,  and  most  recently  Nunez  ap¬ 
plied  M-ary  phase  shift  keying  (MPSK)  to  multiple  access  TDCS  systems  [11]  [6]. 

This  thesis  continues  the  research  on  the  application  of  spectral  encoding  using 
the  TDCS.  The  specific  application  of  spectral  encoding  in  this  research  is  to  over¬ 
come  another  reality  of  the  communications  environment:  the  reflections  that  cause 
multipath  fading. 

1.2  Problem  Statement 

The  multipath  fading  channel  adversely  impacts  wireless  communications.  The 
receipt  of  out-of-phase  reflections  at  the  receiver  degrades  performance  in  the  com¬ 
munications  system.  This  research  applies  spectral  encoding  techniques  to  improve 
performance  over  a  multipath  fading  channel. 

1.3  Goals 

The  specific  overarching  research  goals  are: 

1.  To  compare  the  bit  error  rate  performance  of  spectrally  encoded  “shaped  spec¬ 
trum”  TDCS  signals  to  un-encoded  “flat  spectrum”  TDCS  signals  in  a  frequency- 
selective,  slowly-fading  (FSSF)  multipath  channel  using  an  L-diversity  TDCS 
RAKE  receiver. 
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2.  To  evaluate  the  interference  rejection  capability  of  an  L-diversity  TDCS  RAKE 
receiver  using  spectrally  encoded,  shaped  spectrum  TDCS  signals  in  a  FSSF, 
multipath  channel. 

1 . 4  Scope 

This  research  is  restricted  to  examining  the  spectrally  encoded  TDCS  and  the 
FSSF  multipath  channel  with  and  without  a  narrow  band  interference  source.  The 
research  does  not  delve  into  channel  estimation,  spectral  estimation,  different  types 
of  interference  sources,  different  modulation  techniques  or  multiple  access  techniques. 

1.5  Assumptions 

This  research  is  based  on  the  following  assumptions: 

1.  The  TDCS  transmitter /receiver  pair  use  equal  energy,  antipodal,  binary  phase 
shift  keying.  (BPSK)  signaling. 

2.  The  TDCS  transmitter  is  continuously  transmitting. 

3.  The  real-world  effects  of  radio  frequency  (RF)  hardware  in  the  TDCS  transmit¬ 
ter/receiver  pair  are  ignored.  All  amplifiers,  filter,  and  antennas  are  assume  to 
have  unity  gain. 

4.  The  effects  of  transmission  loss  are  accounted  for  in  signal-to-noise  ratio  (SNR). 

5.  The  TDCS  transmitter /receiver  pair  achieve  perfect  spectral  estimates. 

6.  The  TDCS  transmitter /receiver  pair  perfectly  estimate  the  channel  transfer 
function/impulse  response. 

7.  The  TDCS  transmitter  and  receiver  are  perfectly  synchronized  with  respect  to 
timing  and  spectral  coding  characteristics. 

8.  The  multipath  fading  channel  is  modeled  as  a  Rayleigh  fading  channel. 

9.  The  multipath  fading  channel  model  returns  on  average  100%  of  the  energy 
transmitted  through  the  channel. 
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10.  All  dopplcr  effects  are  ignored. 

11.  Since  all  modelling  occurs  at  baseband,  with  low-pass  equivalent  channel  models, 
the  additive  white  Gaussian  noise  (AWGN)  is  generated  with  a  one-sided  power 
spectral  density  (PSD)  of  N0. 

1 . 6  Methodology 

To  achieve  the  research  goals,  both  analysis  and  simulation  are  used.  Mathe¬ 
matical  analysis,  including  closed-form  and  numerical  solutions  are  used  to  predict 
the  performance  of  the  spectrally  encoded  TDCS  in  various  situations.  The  analysis 
is  verified  using  Monte  Carlo  simulation  techniques. 

1.7  Materials  and  Equipment 

All  simulations  are  run  using  Matlab®  Version  7.0.0.19920  (R14).  Simulations 
are  run  on  stand  alone  PCs  running  Windows  XP  as  well  as  the  AFIT  Unix  Clusters. 

1 . 8  Overview 

This  research  is  organized  as  follows.  Chapter  II  contains  the  literature  search. 
It  provides  the  mathematical  background  on  TDCS  performance  prediction,  multipath 
fading  channel  effects  on  communications  system  performance  prediction  and  the  L- 
diversity  RAKE  receiver.  Chapter  III  contains  the  analysis  and  analytical  results  of 
this  research.  The  background  information  presented  in  the  Literature  Search  is  used 
in  the  analysis  to  predict  the  performance  of  spectrally  encoded  TDCS  signals  over 
multipath  fading  channels.  The  methodology  in  Chapter  IV  describes  the  system 
under  test  and  the  design  of  experiments  for  the  simulation  portion  of  the  research. 
Chapter  V  presents  the  results  of  both  the  analysis  and  simulation.  A  comparison 
of  the  analytical  predictions  to  the  simulated  results  is  discussed  for  each  type  of 
simulation.  Finally,  Chapter  VI  contains  the  conclusions  of  this  research  and  potential 
directions  for  future  studies. 
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II.  Literature  Review 


This  chapter  provides  the  analytical  background  required  for  this  research.  In¬ 
troductions  to  the  TDCS,  the  fading  channel  and  L-cliversity  are  presented. 

2.1  The  Function  of  a  TDCS 

The  basic  premise  of  using  a  TDCS  to  combat  interference  is  that  the  commu¬ 
nications  waveform  can  be  designed  to  have  spectral  characteristics  that  completely 
avoid  the  spectrum  occupied  by  the  interference  source.  It  is  this  spectral  flexibility 
that  is  used  to  enhance  the  performance  of  the  communication  system  in  a  multipath 
environment.  The  basic  function  of  a  TDCS  described  below  shows  how  an  MPSK 
TDCS  generates  communications  waveforms  that  have  desirable  spectral  characteris¬ 
tics. 

Diagrams  of  a  typical  TDCS  transmitter /receiver  pair  are  shown  in  Figures  2.1 
and  2.2.  The  basic  operation  of  a  TDCS  begins  as  a  Discrete  Fourier  Transform  (DFT) 
or  other  spectral  estimation  technique  is  applied  to  determine  the  current  spectral 
environment.  The  measured  magnitude  spectrum  is  processed  to  identify  frequencies 
with  potential  interference  sources  using  one  of  many  algorithms  (e.g.  thresholding). 
Frequency  components  Ap  identified  as  potentially  containing  interference  sources  are 
notched  in  the  construction  of  the  communications  waveform,  effectively  filtering  out 
the  interferer’s  power  over  the  notched  frequencies.  Out  of  the  original  P  frequencies, 
those  which  remain  after  processing  are  phase  encoded  to  generate  the  M  =  2k 
fundamental  signaling  waveforms  (FSW). 

Since  each  of  the  frequency  components  that  are  used  to  construct  the  FSWs 
can  be  assigned  an  independent  phase,  the  information  content  in  the  waveform  is 
encoded  in  the  phase  spectrum.  There  are  multiple  methods  cited  in  the  literature  to 
encode  both  data  symbol  and  multiple  access  information  within  the  phase  [6,9,13]. 
Both  cyclic  shift  keying  and  phase  shift  keying  have  been  used  to  encode  symbol  data. 
For  multi-access  encoding,  pseudo-random  and  linear  phase  codes  have  been  used. 
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Figure  2.1:  Conceptual  TDCS  Transmitter 


Once  the  magnitude  and  phase  spectrum  for  the  FSWs  are  established,  the  spec¬ 
trum  magnitudes  are  scaled  to  maintain  a  constant  signal  energy.  Then,  an  Inverse 
Discrete  Fourier  Transform  (IDFT)  is  performed  on  the  notched,  phase-encoded  and 
scaled  spectra  to  construct  time-domain  representations  of  the  M  FSWs. 

To  modulate  data,  the  TDCS  transmitter  converts  k  data  bits  d*  to  one  of 
M  symbols.  The  FSW  for  each  symbol  is  then  converted  from  digital  to  analog  to 
modulate  each  symbol.  There  is  also  the  potential  to  carrier  modulate  each  symbol 
before  amplification  and  transmission  through  the  channel. 

A  scenario  in  which  a  TDCS  constructs  communications  signals  to  avoid  an 
interference  source  can  be  easily  illustrated.  If  a  typical  interfering  source  has  the 
spectral  characteristics  shown  in  Figure  2.3  [6,  4-3,4],  then  after  utilizing  a  simple 
thresholding  algorithm,  the  magnitudes  of  the  TDCS  frequency  components  Ap  are 
shown  in  Figure  2.4.  The  notched  TDCS  magnitude  spectrum  has  no  energy  in 
frequency  bands  where  the  interferer’s  energy  is  highest.  Note  that  the  example 
spectra  are  two-sided. 
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d(t)  Data 


Figure  2.2:  Conceptual  TDCS  Receiver 

To  detect  and  estimate  transmitted  TDCS  signals,  the  receiver  is  designed  to 
perform  coherent  demodulation.  From  the  diagram  in  Figure  2.2,  it  can  be  seen  that 
the  receiver  uses  the  same  process  as  the  transmitter  to  generate  the  fundamental 
signaling  waveforms.  While  the  receiver  may  or  may  not  ’’see”  the  same  spectral 
environment  as  the  transmitter,  it  is  assumed  that  the  transmitter  and  receiver,  after 
initial  synchronization,  can  communicate  the  required  spectral  characteristics  of  the 
FSWs  to  each  other.  The  receiver  generates  the  M  FSWs  using  the  same  magnitudes 
and  phase-encoding  scheme  as  the  transmitter,  which  then  allows  the  receiver  to 
construct  the  complex  conjugate  of  the  FSWs  and  use  a  matched  filter  demodulator  to 
coherently  detect  and  estimate  symbols.  Nunez’  research  showed  that  using  the  TDCS 
transmitter/reciever  architecture  described  above  allows  accurate  probability  of  bit 
error  prediction  using  standard  MPSK  analytical  models  [6,  4-6].  This  research  uses 
strictly  BPSK,  for  which  the  standard  probability  of  bit  error  (Pb)  relationship  to  the 
normalized  bit  energy  to  noise  power  spectral  density  ( Eb/N0 )  is  given  by  (2.1)  [8,  268]. 
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Figure  2.3:  Conceptual  Interference  Source 


Figure  2.4:  Conceptual  TDCS  Notched  Signaling  Spectra 


Pb 


where  Q(x)  is  defined  as  [8,  40]: 


Q(x)  = 


x  >  0 


(2.1) 


(2.2) 


The  analytical,  continuous-time  model  for  a  baseband  MPSK  TDCS  symbol 
Sk(t)  is  given  by  (2.3). 


p- i 

Sk(t)  =  2  Ap  cos  (27 rfsbpt  +  (j)p  +  9k)  tQ  <  t  <  tQ  +  Tsh  (2.3) 

p= o 
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where  the  following  parameters  are  defined: 


P  is  the  number  of  frequency  components  available  to  construct  the  TDCS  waveform. 

fsb  is  the  fundamental  frequency  of  the  waveform.  All  other  components  are  in¬ 
teger  multiples  of  the  fundamental  frequency  to  ensure  orthogonality.  Note 
also  that  the  period  Tsb  is  equal  to  the  symbol  duration  Tsym  and  therefore, 
Tsb  =  1  /  fsb  =  Tsym. 

Ap  is  the  magnitude  of  the  pth  frequency  component  in  waveform. 

(j)p  is  the  random  phase  for  each  frequency  component  of  the  waveform.  It  is  used 
to  make  the  time-domain  representation  of  the  signal  more  featureless  and  has 
been  used  in  previous  work  to  encode  the  multi-access  phase  assignment.  [6,  3-5]. 

6k  is  the  data  modulation  phase  for  the  kth  symbol  and  is  independent  of  p.  The 
data  phase  is  consistent  across  all  frequency  components  of  the  symbol  for  PSK. 

Assuming  the  system  sampling  frequency  is  an  integer  multiple  of  fsb,  and  equal 
to  2  fsbP  to  avoid  aliasing,  the  discrete-time  model  for  Sk(t)  in  (2.3)  is  shown  in  (2.4)  [6, 
3-2], 

2  P_1  n 

Sk[n\  =  —  Ap  cos  (^2np—  +  <pp  +  0  <  n  <  N  —  1  where  N  =  2 P  (2.4) 

p= o 

From  Equations  2.3  and  2.4,  we  can  determine  the  frequency  domain  representations 
of  the  symbols.  Noting  that  the  continuous-time  representation  of  the  symbol  is  a 
sum  of  cosines  windowed  by  a  rect  function,  the  continuous  frequency  representation 
is  found  in  (2.5). 


p- 1 

St(f)  =  J2Ac 

p= o 


sine 


TsbU-pfsb ) 


,+Mv+Ok)  +  s}nc 


Tsb(f+pfsb) 


(2.5) 


where: 


sinc(a;) 


sin(7nc) 

71X 


(2.6) 
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Making  the  same  sampling  assumption  as  above,  the  discrete-frequency  representation 
of  (2.3)  is  found  in  (2.7)  [6,  3-2], 


Sk[m\  =  { 


^q^g+i(^p+0fc) 
Apt?—  j(.4>p+9k) 


if  m  =  p  Vp  G  (0,  P  —  1) 

if  m  =  N  —  1  —  p  Vp  G  (0,  P  —  1) 

otherwise 


(2.7) 


2.1.1  TDCS  Performance  Against  Interference.  Now  that  the  basic  function 
and  analytical  models  for  the  TDCS  have  been  presented,  the  performance  of  a  such 
a  system  against  an  interference  source  is  examined.  There  are  a  number  of  different 
interferer  models  that  have  been  simulated  against  a  TDCS.  The  ability  of  a  TDCS  to 
suppress  narrow-band  and  partial-band  interference  sources  has  been  demonstrated 
in  both  [6,9].  The  basic  model  for  an  interferer  and  its  effect  on  TDCS  performance 
is  demonstrated  as  follows. 

Using  the  information  presented  earlier  on  the  function  of  a  TDCS,  assume  a 
generic  interference  source  is  received  by  a  TDCS  receiver.  The  discrete-frequency 
representation  of  the  interference  source  is  shown  in  (2.8). 


J[m]  =  Jmejern  Vm  G  [0,  iV  —  1] 


(2.8) 


In  the  discrete  representation  above,  Jm  represents  spectral  magnitude  and  0m  is  the 
spectral  phase.  It  can  also  be  shown  that  the  amount  of  energy  from  the  interference 
source  that  is  received  by  the  TDCS  receiver  in  a  single  symbol  duration  J-  is  given 
by  (2.9). 


a 


N- 

£ 


(2.9) 


However,  as  discussed  earlier,  the  TDCS  is  designed  to  mitigate  the  effects  of  inter¬ 
ference  sources  in  the  frequency  domain.  Several  methods  for  excising  interference 
using  a  TDCS  are  reviewed  in  [9,  6-7].  For  simplicity,  this  research  assumes  a  simple 
thresholding  algorithm  is  used  to  avoid  interference.  The  algorithm  creates  a  notch  in 
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the  TDCS  spectrum  by  setting  Ap  =  0  for  all  components  where  the  spectral  estimate 
exceeds  a  given  threshold.  If  the  system  creates  a  notch  in  the  FSWs  over  spectral 
indices  m  G  ( ji,jh )  and  m  G  (N  —  1  —  jh,  N  —  1  —  ji),  then  the  amount  of  energy  from 
the  interference  source  present  in  the  test  statistic,  Ji  is  given  by  (2.10). 


Ji 


1 

N 


31- 1  N-l-(jh+l)  TV— 1 

E-t2  +  E  E  n 


P= 0  P=jh+ 1  P=N-l-(ji-l) 


(2.10) 


If  the  phases  of  the  interfering  signal  are  assumed  to  be  uniformly  random,  then 
the  energy  in  the  test  statistic  J*  can  now  be  treated  as  AWGN.  Therefore  (2.11) 
is  the  theoretical  prediction  for  the  probability  of  bit  error  P^  in  the  presence  of  an 
interfere!'  [12,  760]. 


Pb 


Q 


2  Eb  \ 

n0  +  nJ 


(2.11) 


where  N,J2  is  the  average  power  of  the  interferer  Pi  divided  by  the  bandwidth  of  the 
TDCS  signal  and  given  by: 


Ni  _  Pi  _  Jj  1  _  Jifsb  _  Ji  (2  12> 

2  2Wbb  Tsym  2Wbb  2  P  fsb  2  P 

The  review  of  the  TDCS  and  the  associated  interference  avoidance  capability  using 
spectral  encoding  techniques  is  now  complete.  Now  the  literature  survey  turns  to  the 
central  focus  of  this  research,  the  fading  channel. 


2.2  Introduction  to  the  Fading  Channel 

One  of  the  basic  assumptions  behind  all  of  the  previous  TDCS  research  at  AFIT 
has  been  that  the  channel  is  strictly  AWGN.  While  this  assumption  is  a  valid  start¬ 
ing  point  for  academic  research,  it  is  not  reflective  of  more  realistic  communication 
environments.  Ionospheric  and  tropospheric  scattering  occur  due  to  constant  physi¬ 
cal  motion  of  ions  in  the  atmosphere  [8,  800].  Additionally,  built-up  structures  and 
terrain  features  (e.g.  mountains),  can  cause  reflections  or  obstructions  in  the  channel. 
The  overall  result  of  scattering,  reflection  and  obstruction  is  that  a  number  of  replicas 
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of  the  transmitted  signal,  called  multipath  waves,  arrive  at  the  receiver  instead  of  the 
single  received  signal  modeled  by  the  AWGN  channel  [10,  139].  Each  of  these  multi- 
path  waves  have  distinct  and  randomly  distributed  amplitude  and  phases,  which  when 
combined  at  the  receiver  cause  constructive  or  destructive  interference  [10,  140]. 

An  added  complication  of  multipath  environments  is  that  they  vary  in  time. 
Ionospheric  scatter  varies  by  the  nature  of  the  motion  of  the  particles  suspended  in 
the  atmosphere  [8,  800].  Mobility  in  urban  communications  environments,  either  by 
the  transmitter/receiver  pair  or  the  obstructions  and  reflectors,  causes  changes  in  the 
channel  over  time  [10,  139].  Overall,  the  time  variation  leads  to  the  characteristics 
of  the  channel  appearing  random  to  the  user  [8,  800].  Next  an  analytical  model  of  a 
multipath  fading  channel  is  introduced. 

To  characterize  the  multipath  channels,  we  look  first  to  the  time-varying  im¬ 
pulse  response.  Since  the  source  analysis  applies  to  bandpass  systems,  the  low-pass 
equivalent  channel  is  analyzed.  The  low-pass  channel  equivalent  can  be  represented 
by  a  sum  of  complex- valued,  time- varying,  delayed  impulse  responses  [8,  802],  As  in 
the  impulse  response  in  (2.13)  and  for  the  remainder  of  this  document,  the  subscript 
l  indicates  that  the  complex  low-pass  equivalent  is  represented. 

L  L 

hi(r,t )  =  y ^cn(t)S(T  -  Tn(t ))  =  y^an(f)eJ0nW<5(f  -  Tn(t))  (2.13) 

n= 1  n=l 


where,  5{t)  denotes  the  Dirac  delta  function  described  by  [15]: 


—  a)  dt 


8(t  —  a) 


/(a) 

0  t  /  a 


(2.14) 

(2.15) 


In  (2.13),  cn(t),  an(t),  @„(t)  and  rn(t)  are  the  time-varying  complex  coefficient,  mag¬ 
nitude,  phase,  and  delay  of  the  nth  multipath  wave  respectively.  In  (2.13),  the  variable 
t  denotes  absolute  time  and  r  is  the  delay  beyond  the  absolute  time.  L  is  the  num¬ 
ber  of  multipath  waves  received.  The  statistics  for  the  coefficients  cn  of  the  channel 
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impulse  response  in  (2.13)  can  be  modeled  using  different  distributions  [8,  803].  To 
limit  the  scope  of  this  research,  the  following  model  is  adopted  for  this  research.  The 
cn(t)  are  assumed  to  be  zero-mean,  complex- valued  Gaussian  random  processes.  Since 
the  magnitude  of  a  such  a  random  process  is  Rayleigh  distributed,  this  model  rep¬ 
resents  a  Rayleigh  Fading  Channel  [8,  803].  Equations  (2.16)  and  (2.17)  show  the 
corresponding  distributions  for  an(t )  and  ©„(£). 

an(t )  ~  Rayleigh(cr)  (2-16) 

Qn(t)  ~  Uniform(— 7T,  7r)  (2-17) 

A  Rayleigh  fading  model  is  typically  used  to  model  situations  where  there  is  no 
dominant  stationary  signal  component  (i.e.  line-of-sight  (LOS)  propagation  path). 
For  example,  in  an  urban  environment,  the  LOS  path  may  arrive  with  a  signal  strength 
on  the  same  order  of  magnitude  as  the  multipath  waves.  This  effect  results  in  a 
Rayleigh  type  fading  channel  [10,  174], 

There  are  a  number  of  channel  characteristics  that  can  help  the  reader  under¬ 
stand  the  effects  of  a  multipath  fading  channel  on  a  communication  system.  These 
characteristics  can  be  derived  from  the  autocorrelation  of  the  channel  impulse  response 
hi(r,t )  and  time- variant  transfer  function  Hi(f,t)  defined  as: 

/OO 

hir^e-^dr  (2.18) 

-OO 

An  extensive  treatment  of  the  autocorrelation  functions  can  be  found  in  Proakis’  text 
[8,  804],  The  most  notable  characteristics  that  describe  the  channel  are  listed  and 
described  here  along  with  the  manner  in  which  they  effect  communications  signals 
traversing  the  channel.  The  first  characteristic  of  interest  is  the  multipath  spread  of 
the  channel,  given  the  label  Tm.  The  multipath  spread  of  the  channel  defines  the 
maximum  length  of  the  delay  observed  in  the  channel.  In  a  communications  system, 
it  is  a  measure  of  the  amount  of  intersymbol  interference  induced  by  the  channel  and 
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is  also  related  to  the  amount  of  diversity  in  the  channel.  The  second  characteristic 
is  the  coherence  ba7idwidth  of  the  channel,  labelled  (A f)c.  The  coherence  bandwidth 
indicates  the  bandwidth  over  which  the  channel  treats  signals  approximately  the  same. 
If  the  baseband  bandwidth  of  a  signal  Wbb  is  greater  than  (A f)c,  then  the  channel 
must  be  considered  frequency-selective.  If  the  channel  is  frequency-selective,  then  the 
analysis  must  consider  how  the  channel  affects  the  communications  system  differently 
in  distinct  parts  of  the  frequency  domain.  If  the  opposite  is  true  ( Wbb  <  (A/)c), 
then  the  channel  is  considered  frequency  non- selective,  and  the  analysis  can  consider 
all  parts  of  the  communications  signal  affected  equally.  The  third  characteristic  is 
the  Doppler  spread,  labelled  B The  Doppler  spread  of  the  channel  is  a  measure  of 
the  amount  of  frequency  spreading  due  to  time  variation  in  the  channel.  The  fourth 
characteristic  is  the  reciprocal  of  the  Doppler  spread,  named  the  coherence  time  of 
the  channel,  and  labelled  (A t)c.  The  coherence  time  is  a  measure  of  the  amount 
of  time  over  which  the  channel  can  be  considered  to  have  time- invariant  behavior. 
Equations  (2.19)  and  (2.20)  show  the  relationships  between  the  four  characteristics 
mentioned  above  [8,  804-807]: 

(A  f)c  «  fr  (2.19) 

(Ai)c  «  2-  (2.20) 

JOd 

2.2.1  How  Multipath  Parameters  Impact  the  Channel  Model.  The  relation¬ 
ship  of  the  four  channel  characteristics  to  the  communications  system  parameters 
define  the  ways  in  which  the  analytical  model  for  the  fading  multipath  channel  may 
be  simplified.  Consider  a  communication  system  with  symbol  duration,  Tsym  and 
baseband  bandwidth,  Wbb,  where  bandwidth  is  defined  as  positive  frequencies  occu¬ 
pied.  If  Tsyrn  »  Tm,  then  in  the  analysis,  intersymbol  interference  can  be  considered 
negligible.  The  selection  of  the  symbol  duration  also  affects  the  consideration  of  time 
variations  in  the  channel.  If  Tsym  «  (A t)c,  then  the  channel  is  considered  slowly- 
fading,  i.e.  the  channel  impulse  response  and  transfer  function  are  considered  to  be 
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fixed  over  the  entire  communications  symbol  duration.  If  this  is  true,  the  time  de¬ 
pendency  can  be  removed  from  the  channel  impulse  response  and  transfer  function 
expressions,  as  shown  in  (2.21)  and  (2.22). 

L  L 

hi(r )  =  y>n<Hr  -  rn)  =  ^aneJ0,A(f  -  rn)  0  <  r  <  Tsym  (2.21) 

n= 1  n= 1 

/oo 

h(r)e-^Tdr  (2.22) 

-CO 

As  mentioned  above,  the  selection  of  the  communication  system  bandwidth  deter¬ 
mines  whether  the  channel  is  considered  frequency-selective  or  nonselective.  If 
Wbb  «  (A f)c,  then  the  channel  is  frequency  non-selective.  A  consequence  of  a 
frequency-nonselective  channel  is  that  the  multipath  components  (reflections)  arriv¬ 
ing  at  the  receiver  are  not  resolvable  [8,  815].  Therefore  the  channel  can  be  modeled 
as  a  single  amplitude  and  phase  term  with  a  random  distribution.  In  contrast,  if 
Wbb  »  (A f)c,  then  the  channel  is  frequency-selective  and  the  individual  reflections 
from  the  channel  are  resolvable.  In  this  case  the  channel  may  be  modeled  as  a  tapped 
delay  line  with  independent  random  amplitudes  and  phases  at  each  tap.  There  is  one 
last  simplification  of  the  channel  model  that  can  be  made  depending  on  the  band¬ 
width  of  the  signal.  Provided  Wbb  »  Bd,  the  effects  of  doppler  shift  in  the  received 
signal  are  negligible  in  the  receiver  [10,  165]. 

2.2.2  Typical  Cellular  Channel  Models.  Multipath  fading  is  a  fact  of  life 
in  cellular  phone  channels  and  manifests  itself  in  many  ways.  However,  the  goal  in 
this  section  is  not  to  perform  an  exhaustive  study  of  channel  characteristics,  but  to 
refer  to  some  typical  measured  values  for  systems  operating  in  the  900MHz  band.  The 
purpose  of  this  quick  overview  is  to  show  that  a  combination  of  channel  characteristics 
and  TDCS  parameters  can  be  conceived  such  that  the  multipath  fading  channel  can  be 
modeled  as  a  frequency-selective,  slowly-fading  multipath  channel.  The  characteristics 
of  multipath  channels  have  been  studied,  measured  and  characterized  by  many  authors 
[8]  [10].  Rappaport  cites  several  characteristics  for  urban  and  suburban  environments 
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measured  in  the  900MHz  frequency  band  [10,  162],  Specifically,  he  presents  measured 
values  for  the  root  mean  square  (RMS)  delay  spread  Tm  as  shown  in  Table  2.1.  These 
statistics  are  shown  with  the  corresponding  coherence  bandwidths  (A f)c  using  (2.19). 


Table  2.1:  Typical  Measured  RMS  Delay  Spread  and  Corre¬ 
sponding  Coherence  Bandwidths 


Environment 

Tm  in  {is  [10,  162] 

(A f)c  in  kHz 

Suburban:  Typical 

0.2-0. 3 

3,300-5,000 

Suburban:  Worst 

1.9-2. 1 

475-525 

Urban:  Typical 

1.3-3. 5 

286-769 

Urban:  Worst 

10-25 

40-100 

Using  these  statistics,  the  TDCS  parameters  can  now  be  determined  so  that  the 
channel  can  be  modeled  as  both  a  frequency-selective  and  slowly-fading  channel. 

First,  the  criteria  for  frequency  selectivity  is  determined.  Since  the  bandwidth 
of  the  system,  Wbb  must  be  much  larger  than  the  coherence  bandwidth,  (A f)c. 
Using  (2.5),  it  can  be  shown  that  the  zero  to  null  baseband  bandwidth,  Wbb  of 
a  TDCS  signal  is  Pfsb  Hz.  The  IS-95  cellular  phone  system  data  rate  is  9.6kHz, 
therefore  as  a  means  of  comparison,  it  is  assumed  that  the  TDCS  fsb  is  10kHz  [4], 
Based  on  this  assumption,  P  in  the  TDCS  symbols  must  be  chosen  such  that  Wbb 
is  greater  than  (A f)c.  Using  the  information  in  Table  2.1,  the  Wbb  for  the  TDCS 
signals  must  exceed  0.7-5MHz  for  the  channel  to  be  considered  frequency-selective. 
Therefore,  if  fsb  =  10kHz,  P  must  greater  than  500  to  meet  this  criteria  in  suburban 
environments,  and  greater  than  76.9  for  urban  environments. 

Next,  the  criteria  for  a  slowly-fading  channel  is  addressed.  Recall  that  this  cri¬ 
teria  is  related  to  the  coherence  time  (A t)c,  which  is  inversely  related  to  the  doppler 
spread  B d-  In  most  cases,  the  significant  amount  of  doppler  spread  in  the  channel 
is  due  to  the  relative  motion  between  the  transmitter,  the  receiver  and  the  entities 
causing  the  delayed  reflections  in  the  environment.  In  this  case,  B ^  can  be  derived 
from  the  relative  velocity  v  between  moving  entities  and  the  wavelength  A  of  the 
communications  signal  using  B<i  =  v/X.  In  Section  2.2,  it  was  shown  that  (Af)c  is 
the  simple  reciprocal  of  B however,  in  practice  this  relationship  yields  the  maxi- 
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mum  (A t)c.  Equations  (2.23)  through  (2.25)  show  two  other  approximations  used  in 
practice,  which  provide  a  range  of  values  for  (A t)c  [10,  165]. 


(A t)c  « 
(A t)c  « 

(A t)c  « 


—  Maximum 
Bd 

9 

-  Minimum 

167T  Bd 


geometric  mean  of  max  and  min  above 


(2.23) 

(2.24) 

(2.25) 


To  meet  the  requirements  for  a  slowly-fading  channel,  (At)c  must  be  many  times 
longer  than  the  symbol  duration,  Tsym  =  Tsb.  Establishment  of  this  criteria  also 
dictates  the  fundamental  frequency  fsi,  of  the  waveform,  as  defined  earlier  in  (2.3). 
To  show  that  the  slowly  fading  assumption  is  reasonable  for  a  typical  cellular  channel 
frequency  and  mobility,  Table  2.2  shows  the  number  of  symbol  durations  contained 
in  (A t)c  given  /sft  is  10kHz,  center  frequency  fc  is  900MHz,  and  relative  velocity  is 
lOOkm/h  for  the  three  different  approximations  given  in  Equations  (2.23),  (2.24)  and 
(2.25).  These  numbers  indicate  that  a  range  of  21-120  symbols  can  be  transmitted  at 
typical  cell  phone  symbol  rates  while  the  channel  characteristics  remain  stable. 


Table  2.2:  Number  of  Symbols  in  Coherence  Time  at 
fsb  =  10kHz,  fc  =  900MHz  and  v  =  lOOkmph 

(A t)c  Approximation  Number  of  Symbols 

Maximum  120 

Mean  50 

Minimum  21 


Using  the  assumptions  above,  the  doppler  spread  for  the  typical  cellular  phone 
channel  at  900MHz  is  computed  to  be  approximately  83Hz.  Given  the  TDCS  band- 
widths  considered  above  (0.7-5MHz)  for  the  frequency  selectivity  discussion,  Rap- 
paport’s  assumption  that  doppler  spread  effects  can  be  ignored  when  the  baseband 
bandwidth  is  much  greater  than  the  coherence  bandwidth  is  reasonable;  therefore, 
the  effects  of  doppler  spread  are  ignored  in  the  channel  model. 
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z(t) 


AWGN 

Figure  2.5:  Tapped  Delay  Line  Model  of  a  Frequency- 
Selective,  Slowly-Fading  (FSSF)  Multipath  Channel 

Based  on  the  information  presented  above,  for  the  remainder  of  this  research, 
it  is  assumed  that  the  channel  behaves  as  a  FSSF  multipath  channel.  Now  the  model 
for  the  channel  can  be  introduced. 

2.2.3  Frequency- Selective,  Slowly-Fading  (FSSF)  Channel  Model.  The  ana¬ 
lytical  model  for  a  FSSF  multipath  channel  is  a  tapped  delay  line  [8,  842],  The  model 
of  the  channel  is  graphically  depicted  in  Figure  2.5,  and  the  time  and  frequency  do¬ 
main  analytical  models  in  (2.26)  and  (2.27)  are  used. 


L  L 

hi(r)  =  ^  cn6(r  -  rn)  =  ^  ane,e"<S(t  -  rn)  0  <  r  <  Tsym  (2.26) 

n= 1  n= 1 

/oo 

h^e-^dr  (2.27) 

-OO 

If  it  is  assumed  that  the  system  is  operating  in  the  digital  domain  in  accordance  with 
the  sampling  assumption  stated  in  Section  2.1,  then  the  resolution  of  the  multipath 
profile  in  the  time  domain  is  as  fine  as  the  sampling  interval,  Tsamp  =  1/ fsamp.  There¬ 
fore,  the  Tn  terms  in  (2.26)  can  be  replaced  with  the  expression  nTsamp.  Here,  L  is 
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related  to  Tm  and  Tsamp  as  stated  in  (2.28)  [8,  841], 


L 


L  Tm/ Tsarnp 


J+l 


(2.28) 


where  |_xj  denotes  the  largest  integer  not  to  exceed  x. 

Now  the  discrete  representation  of  the  impulse  response  hi  [n]  can  be  expressed 
as  shown  in  (2.29). 

L  L 

hi[n]  =  Ci5[n  —  i]  —  Q'ieJ'0i<5[n  —  i\  n  e  (0,  iV  —  1)  (2.29) 

i— 1  i=l 

which  is  given  in  the  following  abbreviated  notation  for  the  remainder  of  the  thesis: 


hi  [n 


cn  =  aneJ~n  n  G  (1,  L), 


(2.30) 


To  express  the  channel  transfer  function  in  the  discrete  frequency  domain  Hi[m],  the 
definition  of  the  discrete  Fourier  transform  (DFT)  is  applied  in  (2.31).  Assume  that 
the  sequence  cn  is  zero-padded  to  length  N  and  that  an  Appoint  DFT  is  calculated. 


Hi[m 


L 


0j2'Knm/N 


cneJ 

n=  1 

r  —  H  pi9” 

rn  —  1  mT 


m  G  (0,  —  1) 


(2.31) 

(2.32) 


where: 


Hm  =  \Cm\  and  6m  =  arg  Crn 


(2.33) 


For  the  remainder  of  the  research,  the  abbreviated  notation  for  the  discrete  channel 
transfer  function  is  Cm  along  with  the  magnitude  and  phase  terms  H.m  and  9m  as 
defined  in  Equations  2.33.  The  stochastic  characterization  of  the  channel  transfer 
function  is  located  in  Section  A.l  of  Appendix  A.  Here  the  result  of  this  derivation, 
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the  distribution  of  Hfn  is  given  in  (2.34). 


- Hl/2Lo 2 

=  2J-g2  (2.34) 

The  distribution  in  (2.34)  is  characterized  as  a  scaled  Chi-squared  random  variable 
with  2  degrees  of  freedom.  The  derivation  of  this  distribution  lends  insight  that  is 
valuable  in  the  analysis  of  bit  error  rate  prediction  over  the  FSSF  multipath  channel; 
however,  even  more  valuable  at  this  point  is  a  graphical  representation  of  the  chan¬ 
nel  impulse  response  and  transfer  function.  A  representative  discrete-time  impulse 
response  for  a  Rayleigh  distributed  channel  with  L  =  16  taps  is  shown  in  Figure  2.6. 
The  associated  channel  transfer  function  computed  by  a  zero-padded  128-point  DFT 
is  shown  in  Figure  2.7.  The  notable  aspect  of  the  channel  transfer  function  plot  is 
the  presence  of  “valleys”  in  the  transfer  function.  These  features  cause  a  frequency 
specific  filtering  of  signals  that  traverse  the  channel.  While  the  plots  included  in 
Figures  2.6  and  2.7  are  only  one  instantiation  of  a  random  channel,  the  filtering  be¬ 
havior  is  a  typical  feature  of  fading  channels,  as  we  can  infer  from  the  distribution 
of  the  channel  transfer  function  derived  in  (2.34).  It  is  this  filtering  behavior  that 
causes  degradation  of  performance  in  typical  communication  systems  and  that  the 
application  of  spectral  encoding  in  the  TDCS  can  mitigate. 

This  completes  the  discussion  of  the  channel  model.  The  graphical  representa¬ 
tion,  time  and  frequency  domain  analytical  models,  associated  statistical  distributions 
and  anecdotal  illustration  of  behavior  have  been  shown. 

2.2.4  Channel  Estimation.  One  of  the  key  assumptions  for  this  research  is 
that  the  channel  is  perfectly  estimated,  therefore  the  model  tap  weights  are  known 
exactly.  As  before,  the  purpose  of  this  section  is  not  to  provide  a  survey  of  all  channel 
estimation  research;  however,  enough  is  presented  to  provide  a  foundation  for  the 
assumptions  made  in  the  research.  The  literature  suggests  it  may  be  assumed  that 
the  channel  coefficients  can  be  accurately  estimated,  provided  the  coherence  time 
of  the  channel  exceeds  100  symbol  durations,  i.e.  (A t)c  >  100 Tsym  [8,  847].  The 
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Channel  Impulse  Response  Magnitude,  L  =  16 


Figure  2.6:  Typical  Multipath  Channel  Impulse  Re¬ 
sponse.  The  impulse  response  above  is  one  instantiation  of 
a  fading  channel.  The  magnitudes  an  are  independently  and 
identically  Rayleigh  distributed. 


Channel  Transfer  Function  Magnitude,  1 28— pt  DFT 


Figure  2.7:  Typical  Multipath  Channel  Transfer  Func¬ 
tion.  The  transfer  function,  computed  using  a  zero-padded, 
128-point  DFT  shows  the  “filtering”  effects  of  the  fading  chan¬ 
nel. 
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literature  also  provides  methods  for  estimating  the  channel  coefficients  when  using 
orthogonal  and  antipodal  signaling  [8,  847]. 

In  addition  to  the  methods  given  in  the  literature,  it  is  postulated  that  a  train¬ 
ing  signal  could  be  transmitted  periodically  that  would  approximate  an  impulse.  A 
transmitted  TDCS  symbol  with  equal-amplitude  frequency  components  and  flat  phase 
response  would  approximate  the  channel  impulse  response  cn  at  the  receiver. 

Now  that  the  multipath  fading  channel  model  and  the  applicable  channel  as¬ 
sumptions  for  this  research  have  been  supported,  a  method  for  predicting  communi¬ 
cation  system  performance  through  the  channel  is  analyzed. 

2.2.5  Predicting  Performance  over  Multipath  Fading  Channels.  To  under¬ 
stand  how  probability  of  bit  error  over  a  multipath  fading  channel  is  predicted,  the 
most  simple  case  is  analyzed.  As  an  initial  step,  the  channel  is  assumed  to  be  fre¬ 
quency  non-selective  and  slowly-fading,  therefore  the  multipath  components  are  not 
resolvable.  In  a  frequency  non-selective  channel,  the  received  signal  appears  as  a 
single  component  with  Rayleigh  distributed  amplitude  a  and  uniformly  distributed 
phase  0  in  addition  to  AWGN  z(t)  [8,  817].  The  analytical  model  for  this  received 
signal  is  provided  in  (2.35). 


ri(t)  =  ae  si(t )  +  z(t)  (2.35) 

While  the  complete  derivation  can  be  found  in  [8],  an  abbreviated  version  is  discussed 
here  as  a  foundation  for  further  analysis.  For  a  matched  filter  demodulator  operating 
in  an  AWGN  channel,  performance  is  dictated  by  the  signal-to-noise  ratio  (SNR)  in  the 
demodulator  [8,  255].  Because  of  the  multipath  channel  effects,  a  distinction  between 
the  SNR  referenced  to  the  transmitted  bit  energy  and  the  SNR  in  the  demodulator  is 
made.  For  the  remainder  of  this  research,  the  normalized  bit  energy  to  noise  power 
spectral  density  ratio  Eb/N0  represents  the  SNR  referenced  to  the  transmitted  bit 
energy.  Equation  (2.36)  is  the  expression  for  the  SNR  in  the  demodulator,  abbreviated 
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(2.36) 


7  b- 


7  b  =  oi2 


N0 


where  a  is  Rayleigh  distributed,  and  a2  has  a  Chi-squared  distribution  with  two 
degrees-of- freedom  (DOF).  As  a  random  variable,  7 b  also  has  a  probability  density 
function  (pdf)  given  by  (2.37)  [8,  817].  In  (2.38),  and  throughout  this  thesis,  E[-]  is 
the  statistical  expected  value  operator. 


p(lb)  = 

— e  76  Z76  where 

7  b 

7 b  = 

IE'“2' 

(2.37) 

(2.38) 


Equation  (2.1)  establishes  Pb  for  BPSK  demodulation.  Since  7 &  is  a  random  variable, 
the  expected  value  of  (2.1)  must  be  taken  over  the  pdf  of  in  (2.37)  to  determine 
the  probability  of  bit  error  in  a  multipath  channel  Pb.mp  as  shown  in  (2.39). 

Q 

The  results  of  this  integration  are  in  (2.40)  [8,  818]. 

(2'40) 

While  it  is  not  immediately  apparent  from  the  analytical  result  of  (2.40),  it  can  be 
seen  from  Figure  2.8  that  Pb  through  a  multipath  channel  with  just  2  delay  taps  is 
much  worse  than  Pb  for  an  equivalent  BPSK  system  at  the  same  Eb/N0  through  an 
AWGN  channel. 

Now  that  the  method  for  finding  the  probability  of  bit  error  for  a  frequency 
non-selective  fading  channel  has  been  shown,  the  focus  of  the  research  returns  to 
the  frequency -selective,  slowly-fading  multipath  channel  introduced  earlier  and  one 
method  for  improving  performance  in  multipath  fading  channels.  Analysis  of  the 
FSSF  multipath  channel  requires  an  adjustment  to  (2.35)  to  account  for  the  multiple 


P 


b,mp 


=  E 


Q(V^) 


lb)p('lb)d'yb 


(2.39) 
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Figure  2.8:  Comparison  of  Pj,  for  BPSK  Demodulators. 

This  figure  demonstrates  degradation  in  the  bit  error  rate  of  a 
BPSK  matched  filter  demodulator  in  at  =  2  multipath  fading 
channel  when  compared  to  an  AWGN  channel 

delay  taps  in  the  model.  This  adjustment  changes  the  distribution  for  7 &;  however, 
the  method  for  computing  P*,  remains  the  same. 


2.3  Diversity  and  the  RAKE  Receiver  Model 

One  way  to  improve  performance  over  hostile  channels  without  increasing  Eb/N0 
is  to  employ  diversity.  Simply  stated,  diversity  is  transmitting  information  more 
than  one  time.  Diversity  can  be  introduced  into  the  system  in  many  ways,  including 
repeating  the  information  on  different  frequencies,  at  different  times  or  within  the 
digital  data  stream  using  block  and/or  convolutional  codes  [8,  821],  As  alluded  to 
earlier,  the  FSSF  channel  can  also  be  a  source  of  diversity. 

The  RAKE  receiver  is  designed  to  leverage  the  replicas  of  the  transmitted  signal 
that  are  received  as  reflections  through  the  multipath  channel  to  achieve  diversity  in 
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the  system.  The  full  derivation  of  RAKE  receiver  performance  is  presented  in  [8, 
842-847].  For  clarity,  an  abbreviated  version  of  the  literature  is  presented  here.  To 
determine  the  performance  of  the  RAKE  demodulator  in  the  context  of  a  BPSK 
system,  the  derivation  begins  with  the  low-pass  equivalent  of  the  transmitted  signal 
sitk{t),  which  has  been  transmitted  through  the  FSSF  and  received  as  r^kit)  in  (2.41). 
Recall  that  the  subscripts  l  and  k  indicate  the  complex  low-pass  equivalent  of  the  kth 
data  symbol  respectively. 


L 

^  ^  wTsamp) 

71=1 


+  zi(t) 


(2.41) 


If  the  signal  contribution  within  r^t)  is  designated  Vi^(t),  the  received  signal  can  be 
simplified  and  rewritten  as  (2.42). 


ri,k(t)  =  vi>k(t)  +  zi(t)  (2.42) 

From  Figure  2.5,  recall  that  zi(t)  is  a  complex  AWGN  process.  If  the  cn  are  per¬ 
fectly  estimated,  as  assumed  in  Section  2.2.4,  then  it  is  apparent  that  the  optimum 
demodulator  is  a  matched  filter  using  v*k(t)  as  the  reference  waveform  instead  of 
s*k(t)  [8,  843].  Upon  closer  examination,  i h,k(t)  is  a  sum  of  L  delayed  and  scaled 
versions  of  the  transmitted  signal  as  seen  in  (2.43). 

L 

U,fc(^)  ^  '  cnsitk(t  nTsamp)  (2.43) 

n= 1 

A  diagram  of  a  typical  RAKE  receiver  is  shown  in  Figure  2.9.  In  the  figure,  the 
test  statistics  of  the  RAKE  receiver  demodulator  Zm  are  formed  by  multiplying  the 
received  signal  by  v*k(t),  integrating  over  a  symbol  duration,  then  examining 

the  real  part  of  the  result  as  shown  in  (2.44)  [8,  843].  Symbol  estimation  in  the  BPSK 
system  is  accomplished  by  taking  the  maximum  of  the  two  test  statistics  and  mapping 
them  to  bits.  Although  the  graphical  and  analytical  representations  appear  complex, 
the  function  of  the  RAKE  receiver  is  identical  to  the  simple  BPSK  matched  filter 
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Figure  2.9:  Diagram  of  a  RAKE  Receiver 


demodulator  with  the  exception  of  the  reference  waveforms. 


Zm  —  Re 


ri,k(t)vlm(t)dt  m  =  1,2 


Tsym 

RAKE  demodulator  performance  is  determined  by  analyzing  the  statistics  of  Zm. 
Two  assumptions  simplify  the  analysis.  First,  it  is  assumed  that  Tm  <C  Tsym,  therefore 
intersymbol  interference  is  ignored  in  the  analysis.  Second,  the  autocorrelation  of  the 
transmitted  signals  or  “self- interference”  is  assumed  to  be  zero  outside  of  i  =  j  as 
in  (2.45)  [8,  846].  These  assumptions  help  simplify  the  analysis  here,  but  are  not 


L 

Y  riik(t)c*ns*l  m{t  -  nTsamp)dt  (2.44) 

71=1 
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applied  to  TDCS  signals  in  the  next  chapter. 


Si,m{t  -  iTsamp)slm(t  -  jTsamp)dt  «  0,  i^j,  m  —  1,2  (2.45) 

After  applying  the  definition  of  ri^if)  to  (2.44),  expanding  the  result  into  nested 
summations,  applying  the  “self- interference”  assumption  from  (2.45),  and  substituting 
«neje'a  for  cn,  it  can  be  shown  that  (2.46)  follows. 


=  Re 


'  L 

£ 

n=  1 


ar 


Sl,k(t  samp) S i  m{t  TlTsamp)  dt 


+ 


Re 


y>we  ■?0 

_n=l 


zi(t)stm(t  -  nTsamp)dt 


(2.46) 


Intuitively,  by  examining  (2.46),  it  can  be  seen  that  the  signal  energy  contribution  in 
Zm  is  a  sum  of  L  scaled  replicas  of  the  transmitted  waveform  s^(f),  and  therefore  the 
FSSF  multipath  channel  becomes  the  source  of  diversity  in  the  receiver.  Note  that 
the  scaling  factors,  while  known  on  a  symbol  by  symbol  basis,  are  random  variables. 
It  is  from  this  result  that  the  derivation  of  performance  for  the  RAKE  receiver  can 
be  determined. 


2.3.1  Probability  of  Bit  Error  for  an  L-Diversity  RAKE  receiver.  The  result 
of  the  previous  section  is  that  the  RAKE  receiver  matched  filter  has  a  test  statistic 
dependent  on  a  sum  of  L  scaled  versions  of  the  transmitted  signal.  Therefore,  the 
received  SNR  for  the  RAKE  receiver  is  a  sum  of  L  scaled  versions  of  the  transmitted 
SNR  Eb/N0  [8,  846].  The  outcome  of  this  development  is  the  expression  in  (2.47)  for 
the  performance  of  an  L-diversity  RAKE  receiver  Pb,RAKE- 


Pb,RAKE 

lb 


Eb  \  ^  2 


£' 


N  ^  an 

1  y  O  -i 

n=  1 


where: 


(2.47) 


However,  as  illustrated  in  Section  2.2.5,  it  must  be  noted  that  7&  is  not  deterministic, 
but  a  sum  of  random  variables.  To  determine  the  bit  error  rate  for  an  L-diversity 
RAKE  receiver,  the  expected  value  of  Pb,RAKE  must  be  taken  over  the  pdf  of  (a  sum 
of  L  independent,  identically  distributed  (i.i.d.),  Chi-squared  random  variables  with 
two  DOF).  After  a  random  variable  transformation,  and  complex  integration,  [8,  825] 
has  shown  the  pdf  of  %  and  expression  for  Pb,RAKE  are  as  shown  in  (2.48)  and  (2.50). 

pin)  =  <2-48> 

where: 

%  =  -|e  [al]  (2.49) 


n  i 

/  L-l  +  k\ 

r  i  1 

Pb,RAI\E  ~ 

[a-g 

5  . 

2(1  +  R) 

(2.50) 


where: 

,l  =  ][t+% 

The  theoretical  improvement  in  bit  error  rate  performance  using  a  RAKE  receiver 
is  shown  graphically  in  Figure  2.10  for  L  =  2.  Clearly,  the  performance  improve¬ 
ment  using  an  /.-diversity  BPSK  RAKE  demodulator  is  better  than  that  of  a  BPSK 
matched  filter  demodulator  without  diversity  through  a  multipath  fading  channel. 
The  performance  of  a  BPSK  matched  filter  demodulator  in  strictly  AWGN  is  also 
shown  for  reference. 

In  this  section  a  technique  for  improving  probability  of  bit  error  performance 
in  a  multipath  fading  channel  has  been  presented.  Specifically,  in  the  case  of  a  FSSF 
multipath  fading  channel,  it  has  been  shown  that  an  L-diversity  RAKE  demodulator 
produces  a  significant  improvement  over  a  simple  BPSK  matched  filter  demodulator. 
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Figure  2.10:  Comparison  II  of  Pb  for  BPSK  Demodula¬ 
tors.  This  figure  demonstrates  the  improvement  in  probability 
of  bit  error  when  an  L-diversity  RAKE  receiver  is  used  in  a 
FSSF  multipath  channel  for  L  —  2. 

2.4  Hypothesis 

The  fundamentals  of  this  research  have  now  been  presented.  To  recap,  the  func¬ 
tion  of  a  TDCS,  the  interference  avoidance  capabilities,  the  multipath  fading  channel 
and  the  L-diversity  RAKE  receiver  have  been  reviewed.  At  this  point,  the  hypotheses 
corresponding  to  the  two  overall  goals  stated  in  Chapter  I  of  this  research  are  now  pre¬ 
sented.  First,  the  spectrally  encoded,  “shaped  spectrum”  TDCS  signals  have  a  lower 
(better)  Pb  than  that  of  the  un-encoded,  “flat  spectrum”  TDCS  signals.  The  second 
hypothesis  is  that  the  shaped  spectrum  TDCS  signals  demonstrate  an  interference 
avoidance  capability  both  with  and  without  the  notching  algorithm  applied. 


2-4-1  Hypothesis  1:  Spectral  Encoding.  This  hypothesis  is  based  on  a  com¬ 
parison  of  the  energy  that  traverses  the  fading  channel  for  spectrally  encoded  and 
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2 


1.8 


Figure  2.11:  Typical  Multipath  Fading  Channel  Transfer 
Function 

un-encoded  TDCS  signals.  The  evaluation  of  the  energy  takes  place  in  the  frequency 
domain  and  is  proved  using  the  Chebyshev  Sum  Inequality.  The  full  proof  is  contained 
in  Appendix  A. 2,  with  the  end  result  (A. 47)  demonstrating  the  energy  received  from 
the  spectrally  encoded,  shaped  spectrum  TDCS  signals  is  greater  by  an  unquantified 
amount  than  the  energy  received  from  the  un-encoded,  flat  spectrum  TDCS  signals. 

& shaped,  ^  £ flat  (2.51) 

For  illustrative  purposes,  consider  a  graphical  example.  Figure  2.11  represents 
the  magnitude  spectrum  of  a  typical  FSSF  multipath  channel  transfer  function.  If  an 
un-encoded,  flat  spectrum  TDCS  signal  with  the  magnitude  spectrum  plotted  in  Fig¬ 
ure  2.12  is  transmitted  through  the  channel  represented  in  Figure  2.11,  the  received 
signal  has  the  magnitude  spectrum  shown  in  Figure  2.14.  In  contrast,  the  spectrally 
encoded,  shaped  spectrum  TDCS  magnitude  spectrum  is  encoded  using  the  shape  of 
the  channel  transfer  function.  The  spectral  magnitudes  are  scaled  to  maintain  equal 
energy  signaling  as  shown  in  Figure  2.13.  It  is  apparent  from  the  graphical  represen¬ 
tations,  that  the  spectrally  encoded  signals  avoid  putting  energy  in  frequencies  that 
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TDCS  Mag  vector 


TDCS  Mag  vector 


Figure  2.12:  Transmitted  Un¬ 
encoded,  Flat  Spectrum  TDCS  Signal 


Figure  2.13:  Transmitted  Spectrally 
Encoded,  Shaped  Spectrum  TDCS  Sig¬ 
nal 


Recevied  Flat  Spectrum  TDCS  Signal 


Received  Shaped  Spectrum  TDCS  Signal 


60 


Figure  2.14:  Received  Un-encoded,  Figure  2.15:  Received  Spectrally  En- 

Flat  Spectrum  TDCS  Signal  coded,  Shaped  Spectrum  TDCS  Signal 


are  “filtered”  by  the  channel,  and  instead  concentrate  signal  energy  in  frequencies 
that  are  “amplified”  by  the  channel.  The  result  is  seen  in  the  received  version  of  the 
spectrally  encoded,  shaped  spectrum  TDCS  signal  in  Figure  2.15.  This  sequence  of 
figures  illustrates  that  the  energy  in  the  spectrally  encoded,  shaped  spectrum  TDCS 
signals  is  expected  to  be  greater  than  that  of  the  un-encoded,  flat  spectrum  TDCS 
signals  by  an  unspecified  amount.  The  quantification  of  this  relationship  is  deferred 
to  the  following  chapter. 


2-4-2  Hypothesis  2:  Interference  Avoidance.  Intuitively,  it  would  be  ex¬ 
pected  that  the  multipath  fading  channel  would  create  “valleys”  in  the  interferer’s 
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spectrum  just  as  it  does  with  the  communication  system.  Therefore,  the  spectrally 
encoded,  shaped  spectrum  TDCS  signals  are  expected  to  exhibit  some  interference 
rejection.  Additionally,  if  the  same  notching  algorithm  used  in  [6]  is  applied  to  the 
spectrally  encoded  TDCS  signals  to  avoid  a  narrow  band  interferer  through  the  fading 
channel,  it  is  expected  that  the  notched  signals  have  better  performance  compared  to 
the  un-notched  signals.  However,  since  “valleys”  may  already  exist  in  the  spectrum, 
the  improvement  may  not  be  as  great  as  seen  over  the  AWGN  channel. 
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III.  Analysis 

In  Chapter  II  the  foundational  principles  and  concepts  for  this  research  were  pre¬ 
sented.  In  this  section,  the  TDCS  concepts  and  principles  of  multipath  fading 
channel  analysis  are  combined  to  analyze  the  performance  of  a  TDCS  system  in  a 
FSSF  multipath  fading  channel. 

3.1  TDCS  Complex  Low-Pass  Equivalent  Model 

The  time  domain  representation  of  a  baseband  MPSK  TDCS  signal  shown 
in  (2.3)  provides  the  basis  for  the  TDCS  analysis.  However,  recall  that  the  analysis  of 
the  multipath  fading  channel  uses  the  complex  low-pass  equivalent  analytical  model. 
The  continuous-time,  complex  low-pass  equivalent,  analytical  model  of  a  transmitted 
TDCS  symbol  si^{t)  in  (3.1)  is  derived  in  Section  A. 3.  In  (3.1),  l  indicates  the  com¬ 
plex  low-pass  representation  and  k  is  in  reference  to  the  kth  of  M  data  symbols.  All 
parameters  are  identical  to  those  described  in  Section  2.1. 

p- i 

si,k(t)  =  Ape^f^p+e^  t0<t<t0  +  Tsb  (3.1) 

p=  —  P 

Assuming  as  before  that  the  continuous-time  signal  is  sampled  at  2 Pfsb  so  that  no 
aliasing  occurs,  then  N  =  2 P.  For  the  remainder  of  the  analysis,  it  is  also  assumed 
that  fSamp  =  2Pfsb-  Using  the  same  derivation  as  in  Section  A. 3,  (3.2)  shows  the 
discrete-time,  complex  low-pass  representation  of  a  TDCS  symbol. 

1  P_1 

sljk[n]  =  —  Apej(2np%+,pp+0k)  0  <  n  <  N  -  1  where  N  =  2P  (3.2) 

p=—p 

Similar  transformations  are  used  to  show  the  complex  low-pass  equivalent  of  the 
continuous-  and  discrete-frequency  representations  of  the  TDCS  symbols  are  given  in 
Equations  (3.3)  and  (3.5)  respectively. 

p- i 

Si,k(f)  =  E  4>{sinc[rsft(/ -pfsb)]ejs^^+ek))  (3.3) 

p=—p 
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where: 


sgn  (p)  =  { 


1  Vj9  >  0 

0  p  —  0 

—  1  Vp  <  0 


(3.4) 


3;  k[m  = 


A  P+i{<t>  m  ~\-0k) 


Vm  G  (0,  P  —  1) 


(3.5) 


AIll-Ne+jl-4lm-N+0^  Vm  e  (P,  N  -  1) 

Note  one  difference  between  equations  in  Chapter  II  and  those  presented  above:  the 
complex  low-pass  equivalent  expressions  allow  for  asymmetric  spectra  and  therefore 
imaginary  components  in  the  symbols. 

Next,  the  energy  in  the  transmitted  TDCS  signal  is  reviewed.  For  a  complex 
valued  signal,  it  is  commonly  known  that  (3.6)  and  (3.7)  give  the  energy  over  a 
single  symbol  duration,  Esym.  Since  this  research  uses  BPSK  signaling  exclusively, 


ESym  Eh . 


ESym  Efj 


/  \s(t)\2dt=  /  s(t)s*(t)dt 

'0  Jo 


(3.6) 


and  by  Parseval’s  Relation  [5,  141]: 


Eb  =  /  \S(f)\2dt=  /  S(f)S*(f)dt 

J—  oo  J—  oo 


(3.7) 


Using  these  relationships,  the  energy  in  the  transmitted,  complex  low-pass,  continuous¬ 
time  TDCS  signal  is  derived  in  Section  A. 4  with  the  following  result: 


p-i 

Eb  =  Tsym  Ap 

p=—p 


(3.8) 


In  the  discrete-time  and  frequency  domain,  the  energy  relationships  are: 


&  =  E 


N- 1  N-l 

2 


S[n\\-  =  y.  s[njs"  [n\ 

n=0  n=0 


(3.9) 
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Figure  3.1:  Notional  TDCS  RAKE  Receiver 


and  again  by  Parseval’s  Relation: 


Eb 


1  N—l  1  N—l 

k= 0  k= 0 


(3.10) 


Using  the  frequency  domain  relationships  in  (3.10)  and  (3.5),  the  transmitted  energy 
for  a  discrete-time,  complex  low-pass  representation  of  a  TDCS  symbol  is  shown  in 
(3.11). 


E>  =  v  U  ^  (3-n) 

k=-P 

3.2  RAKE  Receiver  Design  Applied  to  the  TDCS 

As  seen  in  Chapter  II,  to  improve  performance  through  a  multipath  fading 
channel,  the  communication  system  must  incorporate  diversity.  Here,  the  RAKE 
receiver  design  is  applied  to  a  notional  TDCS  system  as  seen  in  Figure  3.1. 
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The  function  of  the  TDCS  RAKE  receiver  is  very  similar  to  the  original  TDCS 
receiver  shown  in  Figure  2.2.  The  RAKE  receiver  structure  uses  the  L  channel  impulse 
response  coefficients  cn  to  scale  the  delayed  replicas  of  the  FSWs.  The  sum  of  the 
delayed  and  scaled  FSWs  form  the  reference  waveforms  in  the  matched  filters. 

The  components  that  allow  spectral  encoding  are  also  added  to  the  design  shown 
in  Figure  3.1.  If  spectral  encoding  is  applied,  the  scaled  channel  transfer  function  mag¬ 
nitudes  KHp  shape  the  spectrum  of  the  FSWs.  Otherwise,  un-encoded,  flat  spectrum 
signals  use  equal  amplitude  frequency  components  Ap  to  construct  the  FSWs. 

In  this  version  of  the  TDCS  receiver,  the  spectral  estimation  components  simply 
serve  to  identify  spectral  components  to  be  notched  for  interference  avoidance  if  it  is 
activated.  As  before,  the  notched  spectral  components  are  scaled  to  maintain  equal 
energy  signaling. 

3.3  TDCS  Signals  in  a  Frequency- Selective  Slowly-Fading  (FSSF)  Mul¬ 
tipath  Channel 

To  explain  the  power  the  TDCS  system  provides  in  mitigating  multipath  fading, 
the  derivation  of  the  received  signal  is  presented  next  with  the  overall  goal  of  deter¬ 
mining  the  relationship  between  the  spectral  magnitudes  of  the  TDCS  signals  Ap  and 
KHp  and  the  energy  contribution  to  the  test  statistic.  Recall  the  transmitted,  com¬ 
plex,  low-pass  representation  of  the  kth  symbol,  s^k(t)  from  (3.1)  and  the  multipath 
fading  channel  impulse  response  hi(t)  from  (2.26).  To  determine  the  complex  low-pass 
representation  of  the  received  signal  at  the  input  of  the  TDCS  RAKE  receiver  ri^it), 
the  transmitted  signal  is  convolved  with  the  channel  impulse  response  and  combined 
with  AWGN  as  in  (3.12).  In  the  derivation  below,  ®  is  the  convolution  operator  and 
rn  =  nTsamp.  Recall  also  that  the  convolution  of  any  signal  with  a  time-delayed  Dirac 
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delta  function  is  a  delayed  version  of  the  original  signal. 


ri,k(t)  =  siik(t)  <8> 


cn5(t  —  nTt 


samp  ) 


n= 1 


+  zi(t) 


^  ^  CnSl,k(t  flTsamp')  T  Zi(t) 


n=  1 

L  P-1 


EE  cnApe^2n^sbp('t~nTsarnp'>+<^p+e^  +  zi(t) 

n= 1 p=—P 


substituting  cn  =  a^e-70" 


(3.12) 

(3.13) 

(3.14) 


ri,k(t ) 


L  P-1 


EE  Q,^^ei[25r/s6p(t-nTsomp)+</)p+0fc+©n]  _|_ 

n=l p=—P 
Vl,k(t)  +  Zl(t) 


(3.15) 

(3.16) 


where: 


L  P- 1 

vi,k{t)  =  anApej[2nfsbp{t-nTsamp)+(l,p+0k+en]  (3.17) 

n= 1 p=—P 


The  signal  component  is  distinguished  from  the  noise  component  in  ri^(t)  and  is 
denoted  viik(t).  Recall  that  u*fc(i)  is  the  reference  waveform  in  the  matched  filters 
of  the  RAKE  receiver.  The  expressions  ritk(t)  and  v^k{t)  from  (3.16)  and  (3.17)  are 
now  substituted  into  the  results  of  the  previous  development  of  the  L-diversity  RAKE 
receiver  test  statistic  referenced  from  (2.44),  yielding 

f  f'Tsym  'I 

Zi  =  Re  |  J  ri,k(t)vi,i(t)  dt\  7  =  1,2  (3.18) 

(  r  T- sym  P  r^-'sym 

=  Re  l  j  viik(t)vli(t)  dt  +  J  zi(t)vZi(t)  dt\  7  =  1,2  (3.19) 

The  first  term  in  (3.19)  matches  the  form  of  (3.6)  and  is  therefore  the  energy  in  the 
received  signal.  Therefore,  the  energy  contribution  in  the  test  statistic  is  renamed  V 
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and  the  test  statistic  equation  is  rewritten  as 

(  rTsym  'l 

Z,  =  R,e  <jK  +  J  z,it)vlt(t)  dt  >  i—  1,2  (3.20) 

Recalling  that  the  performance  of  a  matched  filter  demodulator  is  directly  related  to 
the  SNR  in  the  test  statistic,  it  is  clear  that  (3.20)  forms  the  basis  for  the  perfor¬ 
mance  analysis  of  the  TDCS  RAKE  receiver.  At  this  point,  the  analysis  diverges  to 
evaluate  the  spectrally  encoded,  shaped  spectrum  TDCS  signals  and  the  un-encoded, 
flat  spectrum  TDCS  signals  separately. 

3.4  Un-encoded,  Flat  Spectrum  TDCS  Signals 

Consider  the  first  of  two  cases  for  TDCS  signals  transmitted  through  the  mul¬ 
tipath  fading  channel.  In  this  case,  the  TDCS  spectrum  is  flat,  i.e.  all  non-zero  Ap 
are  equal.  Assuming  also  that  the  system  uses  equal  energy  signaling  and  that  there 
are  no  DC  or  Nyquist  frequency  components  in  the  TDCS  FSWs,  the  relationship 
between  the  spectral  magnitude  components  and  the  continuous-time  representation 
of  the  transmitted  bit  energy  Eb  is  shown  in  (3.21)  starting  from  the  relationship 
given  in  (3.8). 

p-i 

Eh  =  Tsym  Al  (3.21.) 

p=—p 

If  N  —  2 P  and  A0  =  A_p  =  0,  then  the  summation  reduces  and  yields: 

Eb  =  TsymA2p(N-  2) 

\j  TsymtN-  A  Vp  G  (-P,  P  -  1),  P  ^  0,  —P 

0  p  =  0,-P 

Since  all  of  the  non-zero  Ap  are  equal,  the  derivation  of  V  takes  place  in  the  time 
domain. 


(3.22) 

(3.23) 
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Beginning  with  the  definition  of  the  reference  waveform  in  the  matched  filter 
from  (3.17),  the  derivation  of  V  for  flat  spectrum  TDCS  signals  follows. 


V  = 


Re 

Re 


Uo 


vi,k(t)vli(t)  dt 


(3.24) 


[2tt  fsbp(t-nT3amp)+<t>p  +  Ok+ldn\ 


=  Re 


V  =  Re 


rTsym  f  L  P-1 

/  \  ^  ^  oinApej 

Jo  Ui=iP=-p 

L  P-1 

EE  Q,mAqe-rt2',Tfsbq(t-mT3°-™p)+'t>q+di+eA  j>  dt 

K  m=  1  q=—P 

/r  /  llp  p 

-L  sym  / 

EE  E  E  {  OL  n  OL  m  Ap  A  q  } 

\ n= 1  m=l  p=— P q=— P 

{ej[2nfsbP(t—nTsamp)+<l>p+dk+Qn]e—j[2Trfst)q(t—mTSamp)+<l>q+0i+®  m  J  dt 

/T  /  LLP  P 

EE  E  E  {  Ol  n  OL  m  Ap  A  q  } 

\ n= 1  m=lp=— P q=— P 


(3.25) 


(3.26) 


^ej[^fsb{pt—qt+qmTsamp—pnT3amp)+(j>p—4iq+6k—6i+Qn-Qrn]^  j  ^ 


(3.27) 


Continuing  the  analysis  of  the  quadruple  sum,  assume  that  the  test  statistic  of  interest 
is  being  examined  (i.e.  i  =  k)  and  consider  the  following  four  cases  of  the  quadruple 
sum  in  (3.27): 


1.  Vi :  p  —  q,  m  —  n 

2.  V2:  p  —  q,  m  ^  n 

3.  V3 :  p  7^  q,  m  —  n 

4.  V4:  p  ^  q,  m  ^  n 


3. 4-0.1  Case  I,  VR  p  —  q,  m  =  n.  In  this  case,  it  is  easily  seen  that 
for  the  test  statistic  of  interest  (i.e.  k  —  i)  all  of  the  terms  in  the  exponential  cancel. 
The  quadruple  sum  quickly  simplifies  to  (3.28). 
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Re 


(3.28) 


Ri  = 


OtnAp  dt 


Since  the  expression  in  the  integral  is  independent  of  time,  the  integration  simplifies 
to  TSym-  Recalling  that  all  Ap  are  equal  except  when  p  —  0,  —  P  the  summation  over 
p  can  be  replaced  with  (. N  —  2)Ap. 


V! 


Re 


Re 


P- 1  L 

T.syrn  ^  ^  Ap  ^  ^  OLn 
p=—P  n= 1 

L 

T,m(N-2)AlJ2 

n= 1 


Oil 


(3.29) 

(3.30) 


Now,  the  expression  for  Ap  from  (3.23)  is  substituted  into  the  expression  for 
V\ .  Note  that  all  terms  in  (3.31)  are  real,  therefore  (3.32)  is  the  final  result  of  the 
derivation  of  Vt . 


Ri 


Re  |  Tsym(N  —  2) 

L 

Egym  'y  "J  OLn 


E, 


sym 


Tgym(N  -  2) 


E 

71—  1 


71—  1 


(3.31) 

(3.32) 


The  expression  for  V\  above  is  very  similar  to  that  of  the  signal  to  noise  ratio 
developed  for  the  RAKE  demodulator  from  (2.47).  By  similarity,  the  statistics  of  V\ 
can  be  described  as  a  sum  of  L  Chi-squared  random  variables  with  2  DOF,  having  a 
mean  energy  of  Esym. 


3. 4-0.2  Case  II,  V?:  p  =  q,  m  ^  n.  In  this  case,  the  quadruple  sum 
is  converted  to  a  triple  sum  and  the  time  dependence  is  removed.  The  second  case 
of  (3.27)  is  continued  next. 
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14 


Re 


oT  P—1 

I  J-  sym 

^  -v 

3  p=~p 


L  L 


Ap  «namei[27r/si,Tsamp(m  n)+0n  0m]  dt 


n=  1  ">=1 
m^n 


(3.33) 


Since  there  is  no  time  dependence,  the  integral  simplifies  to  Tsym.  Recalling 
the  relationship  between  Tsamp  =  l/fsamp  =  1/(2 Pfsb),  and  separating  the  channel 
impulse  response  terms,  14  is  expressed  as: 


14  =  Re 


T„m  >  '  A2pJ2J2anamej(&n~em)e^(m-n) 


p- 1 

-sym 


p=—P  n=l  m=1 

m^n 


(3.34) 


The  expression  for  Ap  using  (3.23)  is  now  substituted  into  (3.34). 


14  =  Re 


T. 


E 


sym 


sym 


P-1 


TSym{N  ~  2) 


E  EE 

p=  —  P+ 1  n=l  rn=l 
P^O  m^n 


anamej(en  e^ePp{m  n) 


(3.35) 


Reverting  to  the  previous  notation:  cn  =  ane^&n 


14 


Re 


E 


sym 


(N  —  2) 


P-1 


E  EE 

•=-P+l  n=l  m=l 

p^O  m^n 


iid  RVs 


(3.36) 


Noting  that  within  (3.36),  the  random  variables  cn  and  c*m  are  identically  and  in¬ 
dependently  distributed,  it  is  assumed  that  the  product  cnc^  can  be  treated  as  a 
constant  within  the  summation.  The  product  is  renamed  cnc*m  =  C  and  removed  from 
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the  summation  as  shown  in  (3.37). 


V2  «  Re 


E 


sym 


(N  —  2) 


P-1 


c  E  EE 

P=—P+ 1  n=  1  m=l 
m^n 


e  p 


(3.37) 


The  triple  sum  over  the  exponential  term  is  considered  for  simplification  next.  Sec¬ 
tion  A. 5  derives  the  simplification  of  the  triple  summation  over  p,  m ,  and  n  that 


results  in  (3.38). 


P-l  L  L 


EE 


e^(m-n) 


-L(L  —  2) 


p=  —  P+ 1  n—  1  rn=l 
pzfi  0  m^n 


The  result  from  (3.38),  is  substituted  into  (3.37)  yielding 


Vo  «  Re  <  E 


J  sym 


-L{L-  2) 


N -2 


A, 


sym 


-/.(/.  -2) 
JV-2 


Re  {C} 


(3.38) 


(3.39) 

(3.40) 


To  obtain  the  pdf  of  V2,  the  distribution  of  Re{C}  must  be  determined.  To  sim¬ 
plify  the  analysis,  it  is  assumed  that  the  real  part  of  the  product  Re{C}  =Re{cnc)n} 
in  (3.40)  is  approximately  equal  to  the  product  of  the  magnitudes  of  the  two  inde¬ 
pendent,  complex  Gaussian  random  variables  ( anam ).  Then  it  can  be  shown  that  the 
distribution  of  V2  is  given  by  (3.41)  (See  Appendix  A. 6). 

fM)  =  L  PiRR exp  {ro2  +  (l()  }  dw  (3  41) 


where: 


A 


J  sym 


(«) 


(3.42) 
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3. 4-0. 3  Case  III  and  IV:  p  4  q.  This  section  demonstrates  that 
regardless  of  the  other  factors,  when  p  ^  q,  the  terms  of  the  quadruple  sum  in  (3.27) 
are  zero.  For  the  analysis,  the  case  in  which  m  4  n  is  considered  since  it  is  the  most 
complicated  case.  After  this  derivation,  it  is  easily  inferred  that  the  case  in  which 
m  =  n  also  goes  to  zero.  The  next  steps  in  the  derivation  are  to  simplify  using  i  =  k 
and  separate  the  terms  with  time  dependence  for  the  integration. 


14  =  Re 


^  L  L  P-1  P-1  ,  n 

Q,  Q  [271  f sbTsamp  (rnq  ©m]  l 

n=l  m=l  p=  —  P  q=-P  ^  ' 


P=  — 

m^n  qy £p 


ej[2irfsbt(p~q)] 


(3.43) 


Now  consider  the  evaluation  of  the  time  integral  alone  and  recall  that  Tsym  =  1/ fsb- 


ej[^fsbt{p~q)]  fa  _ 


j2nfsb(p  -  q ) 

1 

j2vr f8b(p  -  q) 
1 

j2vr fsb(p  -  q) 


aj[2nfabt(p-q)i 


J[2 *£*&»-?)]  _  jO 


(3.44) 

(3.45) 


(1-1)  =  0  Vp,g  G  1  (3.46) 


Since  p  and  q  are  integers,  all  terms  in  Cases  Ill  and  IV  are  zero. 


3-4-1  Un-encoded,  Flat  Spectrum  TDCS  Signal  Analysis  Summary.  The 
primary  result  of  this  section  is  that  the  quadruple  sum  from  (3.27)  is  simplified  to 
two  terms  V  and  14-  These  terms  are  simplified  to  the  extent  that  the  test  statistic 
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output  is  restated  in  a  much  simpler  form: 


Zi  =  Vi  +  14  + R,e  | y  vll(t)zi(t)  dt 

L 


ESym  ^  ^  OLn  +  A, 


sym 


n=  1 


-L(L  —  2) 
IV-  2 


Re {C} + 


Re  <(  /  v*ti(t)zi(t)  dt 

I  o 


(3.47) 


(3.48) 


Given  the  expression  for  the  test  statistic  in  (3.48),  the  received  SNR  for  un-encoded, 
flat  spectrum  TDCS  signals  7 /  is  given  by  (3.50). 


7/  = 
7/  = 


R1  +  R2 


A7 


Jsym 

ivT 


L(L-  2) 
N-2 


Re  {C} 


(3.49) 

(3.50) 


Recalling  the  technique  for  predicting  performance  through  a  multipath  fading  chan¬ 
nel  from  Chapter  II,  the  goal  of  this  section  is  to  determine  the  pdf  of  7 f.  The 
following  equations  illustrate  the  process  for  deriving  pdf  of  7 f  using  placeholders  for 
the  densities  of  V\  and  V2.  Since  7 /  is  simply  the  sum  of  V\  and  V2,  scaled  by  l/iV0, 
the  distribution  for  7 f  is  found  using  the  relationships  found  in  [7,  184].  As  before,  <8> 
is  the  convolution  operator. 


z  = 


V1  +  V2 

fz(z )  =  MiVjtofvM) 

/oo 

MWMz  -  x)dx 


(3.51) 

(3.52) 

(3.53) 


The  transformation  to  scale  the  sum  in  (3.53)  is  simply  applied  as  [7,  131]: 


11  =  k 

A,(  v)  =  "ofANo  v) 


(3.54) 

(3.55) 
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Figure  3.2:  Un-encoded,  Flat  Spectrum  TDCS  Signal 
Analytical  Results.  This  figure  depicts  the  analytical  predic¬ 
tions  of  Pb  for  un-encoded  TDCS  signals  demodulated  with  an 
L-diversity  TDCS  RAKE  receiver.  A  clear  improvement  in  I), 
is  evident  as  the  amount  diversity  in  the  receiver  increases  from 
L  —  2  through  L  =  50. 

No  closed  form  solution  for  the  distribution  of  7/  was  achieved.  Therefore  the  analyt¬ 
ical  prediction  of  P5  is  determined  by  numerically  solving  for  /7f  (7/),  then  computing 
the  expected  value  of  Q(a/2t/)  over  the  density  function. 

3-4-2  Un-encoded,  Flat  Spectrum  TDCS  Signal  Analysis  Results.  The  re¬ 
sults  of  the  numerical  computation  of  P, i,  are  shown  in  Figure  3.2.  The  computations 
assume  P  =  1024  and  a  range  of  L  including  2,  4,  20  and  50.  The  analytical  results  in 
Figure  3.2  are  computed  using  numerical  integrations  to  transform  the  random  vari¬ 
able  pdfs  and  compute  the  expected  value  of  the  Pb  equation.  Recall  also,  that  these 
predictions  are  not  expected  to  be  exact  since  they  are  based  on  the  approximations 
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for  V]  and  V2  explained  in  the  sections  above.  Therefore,  these  results  should  closely 
match  the  simulation  results,  but  may  not  be  exact. 

3.5  Spectrally  Encoded,  Shaped  Spectrum  TDCS  Signals 

Chapter  II  introduced  the  concept  of  spectrally  encoding  TDCS  signals  using 
the  magnitude  of  the  channel  transfer  function  as  basis  for  the  spectral  shape  of 
the  TDCS  FSWs.  In  this  section,  the  strategy  is  explicitly  described  and  analyzed. 
The  hypothesis  development  in  Section  2.4.1,  demonstrated  graphically  the  “filtering” 
effect  of  the  FSSF  channel.  The  power  of  a  TDCS  to  improve  performance  in  a  FSSF 
is  similar  to  the  interference  avoidance  strategy.  The  TDCS  is  capable  of  constructing 
FSWs  that  avoid  placing  energy  in  the  frequency  spectrum  where  the  channel  “filters” 
out  energy. 

The  Pb  improvement  strategy  used  in  this  research  to  overcome  the  “filtering” 
effects  of  the  FSSF  is  to  use  the  magnitude  spectrum  of  the  channel  transfer  function 
to  spectrally  encode  the  magnitude  spectrum  of  the  TDCS  FSWs.  Since  it  is  assumed 
that  the  channel  impulse  response  h/(r)  can  be  perfectly  estimated  (c.f.  Section  2.2.4), 
the  computation  of  the  transfer  function  H[(f)  using  the  Fourier  transform  is  assumed 
to  be  an  exact  representation  of  the  channel.  The  spectrally  encoded  magnitudes  of 
the  TDCS  FSWs  are  scaled  to  maintain  equal-energy  signaling  as  shown  in  (A. 35) 
from  Appendix  A. 2.  Through  the  application  of  the  Chebyshev  Sum  Inequality,  the 
derivation  in  the  appendix  shows  that  TDCS  signals  with  the  same  spectral  shape  as 
the  channel  transfer  function  pass  more  energy  through  the  FSSF  multipath  channel 
than  flat  spectrum  TDCS  signals  containing  the  same  amount  of  energy.  In  this 
section,  the  quantitative  analysis  of  the  TDCS  RAKE  demodulator  test  statistic, 
beginning  with  (3.20)  is  presented  for  spectrally  encoded,  shaped  spectrum  TDCS 
signals.  Since  the  spectral  magnitudes  Ap  of  the  FSWs  for  the  shaped  TDCS  signals 
are  not  all  equal,  the  analysis  takes  place  in  the  frequency  domain. 
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To  continue  the  quantitative  analysis,  the  discrete-frequency  representation  of 
the  received  shaped  spectrum  TDCS  signal  Ritk[m\  referenced  from  (A. 38)  follows. 


Ri,k[m 


'  KHie+j^+e^  if  m  =  p  Vp  e  (0,  P  -  1) 

<  KH^e~j^+e^  if  m  =  N  -  1  -  p  Vp  e  (0,  P  -  1) 
0  otherwise 


where: 


(3.56) 


K 


E  N 

^  sym1  v 


VP_1  H2 

l^V=-P  P 


(3.57) 


Using  (3.11)  and  the  results  above,  the  energy  in  the  received,  shaped  spectrum  TDCS 
signals  Vshaped  is  derived  next. 


U 


shaped 


P-1 


E  K 


i2 

N  ^  -  p 

p=  —  P 


P-1 


v  E  (™~>2 


p=  —  P 


E  N 

^  sym1  v 


=  E 


X^P~ 1  Zj4 
Z^p=-p  np 


sym , 


y-P-i 

zEfc=-P  k 


■  P-1 

E  Ht 

p=—p 


(3.58) 

(3.59) 

(3.60) 

(3.61) 


The  expression  for  Vshaped  from  (3.61)  is  now  substituted  for  V  in  the  test  statistic 
equation  referenced  from  (3.20)  as  shown  in  3.62. 


E,  = 


Re  <  E 


EP_1 

Z-^p=— 


ttA 
P  iJ-p 


sym 


Yp-i  H2 
<E,k=-P  nk 


+ 


ViAt)z(t)  dt 


(3.62) 

(3.63) 


While  the  individual  pdfs  for  H2  and  H 1  can  be  derived  (See  Appendix  A. 7),  each  of 
the  N  E[p  random  variables  are  scaled  sums  of  the  same  impulse  response  coefficients 
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Ci, ... ,cl .  Therefore,  these  are  not  sums  of  independent  random  variables,  and  cannot 
be  simplified  without  the  joint  density  functions  of  all  N  H 2  and  H2. 


3.5.1  Estimation  of  Spectrally  Encoded,  Shaped  Spectrum  Signal  Energy  Dis¬ 
tribution.  Since  the  joint  densities  of  the  sequences  H2  and  are  unknown, 
further  development  of  the  closed  form  analytical  solution  is  not  possible.  Instead, 
Monte  Carlo  simulations  and  distribution  parameter  estimation  methods  are  used  to 
estimate  the  distribution  of  the  ratio  in  (3.64). 


E 

E 


p-i 
p=—p 
p- 1 
p=~p 


H2 


(3.64) 


Monte  Carlo  simulations  of  the  random  sequence  cn  are  the  basis  for  the  computations. 
From  each  L-length  sequence  of  cn,  the  resultant  frequency  domain  terms  H2  and 
are  computed  using  a  zero-padded,  Appoint  DFT  (. N  =  2 P).  The  summations  of 
and  H2  are  then  computed  and  divided  as  in  (3.64).  The  simulation  terminates  when 
the  95%  confidence  interval  on  the  sample  mean  of  the  division  converges  to  within 
0.1%  of  the  computed  value  (see  Appendix  A. 8  for  calculations).  After  the  simulation 
converges,  the  data  is  run  through  Matlab®  parameter  estimation  routines  for  several 
distributions.  Based  on  visual  inspection  of  the  histograms,  the  Monte  Carlo  data  is 
tested  against  the  Gamma,  Rayleigh,  Log-Normal  and  Weibull  distributions.  Once 
the  parameters  for  each  of  these  four  distributions  are  estimated  by  the  Matlab® 
routines,  they  are  compared  to  determine  the  distribution  that  most  closely  matches 
the  simulation  data.  Quantile-quantile  plots  with  linear  regression  models  determine 
which  distribution  best  fits  the  simulation  data. 

The  quantile-quantile  plot  is  a  standard  method  for  testing  the  distribution  of 
data  [3,  196].  The  plot  provides  a  visual  comparison  of  the  simulated  data  quantiles 
and  the  known  quantiles  of  a  given  distribution.  If  the  distributions  are  identical,  the 
plot  should  yield  a  straight  line.  Therefore  a  linear  regression  model  of  the  quantile- 
quantile  plot  and  analysis  of  variance  on  the  linear  model  provide  a  quantifiable 
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measure  of  the  quality  of  the  distribution  match.  Appendix  A. 9  contains  a  detailed 
description  of  the  process  used  to  develop  the  quantile-quantile  plots,  linear  regression 
model  and  the  coefficient  of  determination.  The  coefficient  of  determination  K2  is  the 
ultimate  measure  of  merit,  with  R2  =  1  indicating  a  perfect  distribution  match. 


3.5.2  Spectrally  Encoded,  Shaped  Spectrum  TDCS  Signal  Summary.  Recall 
that  finding  the  distribution  of  the  division  in  (3.64)  is  only  an  interim  step  toward 
determining  the  performance  of  the  spectrally  encoded,  shaped  spectrum  TDCS  sig¬ 
nals.  Once  the  distribution  is  determined,  the  predicted  Pb  is  estimated  using  the 
same  technique  illustrated  in  Chapter  11.  First,  the  pdf  selected  after  the  parameter 
matching  process  is  scaled  by  the  transmitted  SNR  Esym/N0  to  determine  the  pdf  of 
the  received  SNR  of  the  spectrally  encoded,  shaped  spectrum  TDCS  signals  %.  Then 
the  probability  of  bit  error  for  the  shaped  spectrum  signals  Pb, shaped  is  the  expected 
value  of  the  Pb  equation  computed  over  the  pdf  of  the  shaped  spectrum  signals  as 
shown  in  (3.66). 


7® 


P, 


b, shaped 


J-^sym  Z-^v= 


N—l  jt4 

p= 0  np 


Nn 


EiV  - 

v= 


N~1  TT2 
p= 0  11p 


Q[  V27s  p(7s)o?7. 


(3.65) 

(3.66) 


The  probability  of  bit  error  is  computed  numerically  since  a  closed  form  solution  has 
not  been  found. 


3.5.3  Spectrally  Encoded,  Shaped  Spectrum  TDCS  Signal  PDF  Estimation  Re¬ 
sults.  Parameter  estimations,  quantile-quantile  plots,  linear  regression  models  and 
coefficients  of  determination  are  determined  for  L  —  2,  4,  20  and  50  given  P  =  1024. 
In  this  section,  the  distribution  test  results  are  presented  by  order  of  increasing  L. 
The  results  for  each  L  contain  quantile-quantile  plots  and  R2  values  for  each  distri¬ 
bution  comparison.  The  parameter  values  for  the  matched  distribution  are  provided 
along  with  the  95%  confidence  intervals  on  each  parameter. 
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Figure  3.3:  Estimation  of  Distributions  for  L  =  2.  These  quantile-quantile 

plots  and  their  associated  R 2  values  demonstrate  that  the  Gamma  distribution  most 
closely  matches  the  received  energy  for  L  =  2. 

3.5.3. 1  Spectrally  Encoded  Distribution  Results:  L  =  2.  In  Figure  3.3, 
quantile-quantile  plots  are  displayed  for  the  four  distributions  mentioned  in  Sec¬ 
tion  3.5.1.  The  quantile-quantile  plots  and  the  coefficient  of  determination  indicate 
that  the  Gamma  distribution  most  closely  matches  the  Monte  Carlo  simulation  data. 
Table  3.1  contains  the  estimated  parameters  for  the  distribution  with  95%  confidence 
intervals  and  the  coefficient  of  determination. 


Table  3.1:  Parameter  Estimation  for  L  =  2 


Distribution 

a  (95%  Cl) 

/3  (95%  Cl) 

R2 

Gamma 

1.9156  (1.867  -  1.9654) 

0.6797  (0.67084  -  0.71135) 

0.99972 

3. 5. 3. 2  Spectrally  Encoded  Distribution  Results:  L  =  4.  Figure  3.4 
contains  the  comparison  of  the  four  proposed  distributions  to  the  Monte  Carlo  sim- 
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Figure  3.4:  Estimation  of  Distributions  for  L  =  4.  These  graphics  illustrate 
the  parameter  matching  quantile-quantile  plots  and  R 2  values  for  L  —  4.  The  Gamma 
distribution  most  closely  matches  the  simulation  data  for  L  —  4. 

ulations  of  the  shaped  spectrum  energy  for  L  =  4.  For  L  =  4,  the  distribution 
estimation  process  indicates  that  the  Gamma  Distribution  is  closest  to  the  simulated 
data  distribution.  Table  3.2  contains  the  estimated  distribution  parameters  with  95% 
confidence  intervals  and  the  coefficient  of  determination. 


Table  3.2:  Parameter  Estimation  for  L  =  4 


Distribution 

a  (95%  Cl) 

(5  (95%  Cl)  R2 

Gamma 

3.6046  (3.5102  -  3.7015) 

0.44444  (0.43196  -  0.45728)  0.99894 

3. 5. 3. 3  Spectrally  Encoded  Distribution  Results:  L  =  20.  Figure  3.5 
contains  the  comparison  of  the  four  proposed  distributions  to  the  simulations  of  the 
shaped  spectrum  energy  for  L  —  20.  The  quantile-quantile  plots  and  coefficients  of 
determination  indicate  that  the  best  matched  distribution  for  L  =  20  is  the  Log- 
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Figure  3.5:  Estimation  of  Distributions  for  L  =  20.  These  quantile-quantile 

plots  illustrate  the  parameter  matching  process.  The  R 2  value  indicates  the  Log- 
Normal  distribution  is  closest  to  the  simulated  data  distribution. 

Normal  Distribution.  Table  3.3  contains  the  estimated  mean  and  standard  deviation 
with  95%  conhdence  intervals  and  the  coefficient  of  determination. 


Table  3.3:  Parameter  Estimation  for  L  =  20 


Distribution 

/i  (95%  Cl)  a  (95%  Cl)  i?2 

Log-Normal 

0.60384  (0.59857  -  0.60911)  0.26877  (0.26509  -  0.27254)  0.99962 

3. 5. 3. 4  Spectrally  Encoded  Distribution  Results:  L  =  50.  The  predic¬ 
tion  of  the  signal  energy  distribution  for  L  =  50  is  shown  in  Figure  3.6.  From  the 
quantile-quantile  plots  and  coefficients  of  determination,  the  best  matched  distribu¬ 
tion  for  the  L  —  50  simulations  is  clearly  the  Log-Normal  Distribution.  Table  3.4 
contains  the  estimated  mean  and  standard  deviation  with  95%  conhdence  intervals 
and  the  coefficient  of  determination. 
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Figure  3.6:  Estimation  of  Distributions  for  L  =  50.  The  quantile-quantile  plots 
and  R2  values  indicate  that  Log-Normal  distribution  is  closest  to  the  simulated  data. 


Table  3.4:  Parameter  Estimation  for  L  =  50 


Distribution 

H  (95%  Cl) 

a  (95%  Cl)  R2 

Log-Normal 

0.6572  (0.6538  -  0.66061) 

0.17372  (0.17134  -  0.17616)  0.99991 

3.5.4  Spectrally  Encoded,  Shaped  Spectrum  Analysis  Pb  Results.  Using  the 
estimated  distributions  for  each  level  of  L ,  the  predicted  Pb  for  the  spectrally  encoded, 
shaped  spectrum  signals  is  computed  numerically.  The  overall  results  are  shown 
in  Figure  3.7.  Recall  that  these  predictions  are  based  on  estimated  signal  energy 
distributions.  The  R2  values  from  the  parameter  estimation  results  indicate  that 
the  distributions  are  very  closely  matched,  therefore  it  is  expected  that  these  results 
should  closely  match  the  simulation  results. 
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Figure  3.7:  Spectrally  Encoded,  Shaped  Spectrum 
TDCS  Signal  Analytical  Results.  This  figure  depicts  the 
analytical  predictions  of  P&  for  spectrally  encoded  TDCS  signals 
demodulated  with  an  /.-diversity  TDCS  RAKE  receiver.  These 
results  also  show  a  clear  improvement  in  as  the  amount  di¬ 
versity  in  the  receiver  increases  from  L  —  2  through  L  =  50. 

3.6  Comparison  of  Un-encoded  and  Spectrally  Encoded  Analytical  Re¬ 
sults 

Observing  the  analytical  Pb  curves  independently  for  the  un-encoded,  flat  spec¬ 
trum  and  spectrally  encoded,  shaped  spectrum  TDCS  signals  clearly  demonstrates  the 
benefit  of  increasing  diversity.  However,  to  see  the  benefits  of  spectral  encoding,  the 
results  are  compared  for  each  level  of  L.  In  Figures  3.8  and  3.9,  each  curve  compares 
the  analytically  predicted  P}}  for  both  the  spectrally  encoded,  shaped  spectrum  signals 
and  the  un-encoded,  flat  spectrum  signals  for  each  level  of  L  independently.  In  each 
case,  the  Pb  for  the  spectrally  encoded  signals  is  lower  (better)  than  the  un-encoded 
Pb .  These  curves  provide  the  basis  for  the  comparison  of  the  analytical  and  simulation 
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Transmitted  E^/No  (dB) 


Figure  3.8:  Analytical  Comparison  for  Spectrally  Encoded  and  Un¬ 
encoded  Signals,  L  —  2,  4.  These  figures  demonstrate  the  benefit  of  spectral 
encoding  applied  to  the  RAKE  receiver  at  equivalent  diversities  L. 
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Transmitted  E^/No  (dB) 


Figure  3.9:  Analytical  Comparison  for  Spectrally  Encoded  and  Un¬ 
encoded  Signals,  L  =  20,  50.  These  figures  demonstrate  the  benefit  of  spectral 
encoding  applied  to  the  RAKE  receiver  at  equivalent  diversities  L. 
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results.  It  is  useful  at  this  point  to  define  the  term  gain  G  to  provide  a  comparison 
of  the  transmitted  Eb/N0  required  to  achieve  a  given  probability  of  bit  error,  A.  If 
Ef/N0  is  required  to  achieve  a  Pb  of  A  for  the  un-encoded,  flat  spectrum  signals,  and 
Esh/N0  is  required  to  achieve  the  same  Pb  of  A  for  the  spectrally  encoded,  shaped 
spectrum  signals,  then  gain  G  is  defined  in  (3.67).  The  gain  G  is  only  determined  for 
signals  with  equivalent  diversities  L.  In  Figures  3.8  and  3.9,  the  gains  are  estimated 
by  graphically  measuring  the  plot. 


G(dB) 


f  <dB>  -  f  <dB> 


(3.67) 


This  concludes  the  analysis  of  the  spectrally  encoded  and  un-encoded  TDCS 
signals  through  the  FSSF  multipath  channel.  Next  the  analysis  to  determine  the  Pb 
of  spectrally  encoded,  shaped  spectrum  signals  is  determined  in  the  presence  of  an 
interferer. 


3.7  Interference  Avoidance 

The  second  goal  of  this  research  is  to  evaluate  the  interference  avoidance  capa¬ 
bility  of  spectrally  encoded  signals.  Chapter  II  introduced  the  interference  avoidance 
capability  of  the  TDCS  with  the  basic  result  that  the  spectral  notching  algorithm  can 
remove  interference  energy  from  the  test  statistic  by  shifting  the  spectral  distribution 
of  energy  in  the  TDCS  FSWs.  Additionally,  the  literature  search  revealed  a  method 
for  predicting  the  Pb  for  notched  and  un-notched  TDCS  signals.  Un-notched  signals 
consider  the  average  power  in  the  interference  signal  P',  compute  the  full  interference 
PSD  N[  and  add  it  to  the  noise  PSD  Na  to  compute  Pb, un-not.ched  as  in  (3.68). 


Pi 


,un— notched 


Q 


2  Eh 


N0  +  N' 


(3.68) 
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where: 


2  Wbb  Tsym  2 P  fsb 

A 

2  P 


(3.69) 

(3.70) 


The  total  energy  in  the  interferer  J[  is  calculated  using  (2.9).  Notched  signal  Pb  is 
predicted  using  the  interferer’s  power  outside  of  the  spectral  notch  P,  to  compute  the 
notched  interferer’s  PSD  Nt. 

Pb, notched  =  Q  ^ \j  /y^_|_  \’.  ^  (3-71) 

where: 


2 WbB  Tsym.  2 P fsb 

A 

2  P 


(3.72) 

(3.73) 


The  expression  for  the  energy  in  the  notched  interference  spectrum  J*  is  given  by 

(2.10). 

Because  the  FSSF  multipath  channel  affects  both  the  shape  of  the  spectrally 
encoded  signals  and  the  effect  of  the  notching  algorithm,  this  section  considers  two 
cases  of  interference  avoidance.  In  the  first  situation,  the  interferer  is  close  enough 
to  the  TDCS  receiver  that  the  interferer  does  not  experience  multipath  fading.  This 
scenario  demonstrates  the  effects  of  spectral  encoding  on  interference  avoidance.  The 
second  situation  distances  the  interference  source  from  the  TDCS  receiver  such  that 
the  interfering  signal  experiences  multipath  fading.  This  situation  demonstrates  the 
effect  of  the  channel  on  the  notching  algorithm. 


3. 7. 1  Interference  Avoidance  I:  No  Multipath  on  Interference  Channel.  This 
section  presents  an  analysis  of  the  situation  in  which  the  TDCS  transmitter /receiver 
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Portions  of  Spectrum 
Exceeding  Notching  Threshold 


Figure  3.10:  Conceptual  Transmitted  Interference 
Spectrum.  This  figure  demonstrates  the  effect  of  the  notching 
algorithm  when  no  multipath  fading  occurs  on  the  interference 
channel. 

pair  traverses  a  FSSF  multipath  channel,  but  the  interferer  passes  through  an  AWGN 
only  channel.  This  situation  represents  the  scenario  in  which  the  interference  source 
is  located  in  close  proximity  to  the  TDCS  receiver.  This  case  is  more  likely  for  high 
interferer  power  to  signal  power  ratios  J/S,  since  the  interferer  may  have  to  be  much 
closer  to  the  receiver  to  achieve  the  higher  power  ratio.  Since  the  spectral  shape  of  the 
interferer  is  not  changed  by  the  interference  channel,  the  notching  algorithm  performs 
as  discussed  in  the  previous  section. 

However,  recall  that  the  TDCS  signals  traverse  a  multipath  fading  channel  and 
are  spectrally  encoded  with  the  shape  of  the  channel  transfer  function.  The  spectral 
shaping  should  “filter”  the  interference  source  at  the  receiver.  The  precise  analytical 
prediction  of  this  effect  and  the  resulting  P\ \  are  unknown  at  this  point,  so  graphical 
illustrations  are  presented  to  provide  intuitive  analysis.  Figure  3.10  depicts  the 
transmitted  magnitude  spectrum  of  a  narrow  band  interference  source.  The  shaded 
area  is  the  portion  of  the  spectrum  notched  out  by  the  interference  avoidance  algo¬ 
rithm.  The  figure  demonstrates  that  the  notching  algorithm  performs  as  predicted 
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Figure  3.11:  Un-notched,  Spectrally  Encoded,  Shaped 
Spectrum  Signal,  Scenario  I.  This  demonstrates  the  “filter¬ 
ing”  effect  created  by  the  spectrally  encoded,  shaped  spectrum 
TDCS  signal.  The  shaded  areas  indicate  areas  of  peak  power 
for  the  interference  source,  that  would  be  notched  by  the  inter¬ 
ference  avoidance  algorithm. 


Figure  3.12:  Notched,  Spectrally  Encoded,  Shaped 
Spectrum  Signal,  Scenario  I.  This  figure  demonstrates  that 
in  addition  to  notching,  the  spectral  shape  of  the  spectrally  en¬ 
coded  TDCS  signal  has  a  “filtering”  effect  on  the  interference 
source. 
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when  the  interference  channel  is  AWGN  only,  attenuating  the  interferer’s  energy  in 
the  receiver  test  statistic  according  to  the  relationship  in  (2.10).  However  this  is  not 
the  only  effect  to  be  considered.  The  spectral  shape  of  the  encoded  TDCS  signals 
also  affects  the  amount  of  interference  energy  in  the  test  statistic. 

To  illustrate  this  effect,  Figure  3.11  represents  the  magnitude  spectrum  of  a 
spectrally  encoded  TDCS  signal  after  being  shaped  by  a  typical  FSSF  multipath 
channel  transfer  function.  In  this  example,  the  diversity  in  the  channel  L  =  20. 
Even  when  un-notched  by  the  interference  avoidance  algorithm,  the  spectrally  shaped 
TDCS  signals  have  a  “filtering”  effect.  The  spectral  notches  in  the  TDCS  signal 
spectrum  should  decrease  the  interferer’s  energy  present  in  the  test  statistic  although 
no  notching  is  applied.  Therefore,  the  simulated  Pb  for  the  un-notched,  spectrally 
encoded,  shaped  spectrum  TDCS  signals  is  predicted  to  be  slightly  less  (better)  than 
the  analytical  results. 

The  corresponding  magnitude  spectrum  of  a  spectrally  encoded  TDCS  signal 
after  notching,  is  shown  in  Figure  3.12.  As  stated  above,  the  notching  algorithm 
removes  the  analytically  predicted  amount  of  interference  energy.  However,  the  figure 
also  demonstrates  that  the  shape  of  the  TDCS  spectrum  is  also  likely  to  remove 
additional  interference  energy.  Therefore,  the  simulated  Pb  for  the  notched,  spectrally 
encoded,  shaped  spectrum  TDCS  signals  is  predicted  to  be  slightly  less  (better)  than 
the  analytical  results. 

To  sum  up  the  first  scenario,  the  analytically  predicted  results  using  (3.71)  and 
(3.68)  are  predicted  to  be  upper  bounds  for  the  simulated  results  for  the  notched  and 
un-notched  spectrally  encoded  signals  respectively. 

3.7.2  Interference  Avoidance  II:  Independent  Multipath  Channels.  This  case 
represents  the  scenario  in  which  the  interference  source  is  not  located  in  the  same 
proximity  as  the  transmitter.  In  this  situation,  the  multipath  fading  channel  impulse 
response  from  the  interferer  appears  different  from  the  impulse  response  coming  from 
the  TDCS  transmitter  when  received  at  the  TDCS  receiver.  Therefore,  it  is  assumed 
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that  the  interferer  and  the  TDCS  transmitter  traverse  independent  channels  before 
arriving  at  the  TDCS  receiver. 

As  stated  in  the  previous  description  of  the  TDCS  interference  avoidance  capa¬ 
bility,  the  TDCS  samples  the  received  spectrum  and  notches  out  frequency  compo¬ 
nents  exceeding  the  threshold.  In  this  scenario,  the  TDCS  receiver  is  estimating  the 
spectrum  of  the  interference  source  after  it  has  passed  through  the  FSSF  multipath 
channel.  The  interference  channel  will  distort  the  spectral  magnitudes  of  the  inter¬ 
ferer,  therefore  the  TDCS  notching  algorithm  will  not  perform  as  previously  predicted. 
Since  the  FSSF  multipath  channel  is  not  deterministic,  the  channel  may  present  either 
a  “filter”  or  “amplification”  in  the  location  of  the  jammer’s  spectral  magnitude  peaks. 
If  the  same  narrow  band  interferer  shown  in  Figure  3.10  is  transmitted  through  an 
independent  channel  with  the  transfer  function  in  Figure  3.13,  the  resulting  received 
interference  spectrum  is  shown  in  Figure  3.14.  From  this  example,  it  is  apparent 
that  the  number  of  notched  components  will  change  clue  to  the  distortion  of  the  inter¬ 
ference  spectrum.  This  example  illustrates  that  the  number  of  notched  components 
may  be  either  increased  and/or  decreased  due  to  either  the  amplification  or  filter¬ 
ing  effects  of  the  multipath  channel.  Additionally,  it  is  apparent  that  the  interferer 
magnitudes  change.  However,  these  graphics  are  one  representation  of  a  stochastic 
channel.  Since  it  is  assumed  that  on  average,  the  channel  passes  through  100%  of 
the  energy  transmitted,  it  is  predicted  that  on  average  the  notching  will  perform  as 
previously  described. 

The  combined  effect  of  the  notching  algorithm  and  spectral  shaping  when  the 
communications  and  interference  channels  are  independent  is  shown  in  Figures  3.16 
and  3.15.  Since  the  communications  channel  is  independent,  the  TDCS  signal  shape 
is  distinct  from  the  interference  distortion.  Using  the  previously  discussed  received 
interference  spectrum  shown  in  Figure  3.14,  the  frequency  components  identified  for 
notching  (the  interferer’s  peak  power  spectra)  are  shaded  in  Figure  3.15.  In  this  ex¬ 
ample,  the  shaded  areas  within  Figure  3.15  happen  to  fall  within  areas  that  have 
been  “filtered”  by  the  spectral  shaping.  However,  depending  on  the  instantiation  of 
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Interference  Channel  Magnitude  Response 


Figure  3.13:  Conceptual  Interferer  Channel  Response. 

This  is  the  magnitude  spectrum  of  a  typical  FSSF  multipath 
channel  transfer  function  with  L  =  20. 


Number  of  Notched  Number  of  Notched 

Components  Increased  Components  Decreased 


Figure  3.14:  Conceptual  Received  Interference  Spec¬ 
trum.  After  being  transmitted  through  the  FSSF  multipath 
fading  channel,  the  interference  spectrum  is  distorted.  The 
shaded  areas  indicate  the  portions  identified  for  notching  by  the 
interference  avoidance  algorithm. 
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Figure  3.15:  Un-Notched,  Spectrally  Encoded,  Shaped 
Spectrum  Signal,  Scenario  II.  The  shape  of  this  spectrally 
encoded  TDCS  signal  is  shaped  by  a  FSSF  multipath  channel 
that  is  independent  of  the  interference  channel.  The  shaded  ar¬ 
eas  indicate  the  frequencies  identified  for  notching  after  receiving 
the  distorted  interference  spectrum.  This  figure  demonstrates 
the  “filtering”  and  “amplifying”  effects  of  spectral  encoding  that 
is  independent  of  the  interference  distortion. 


Figure  3.16:  Notched,  Spectrally  Encoded,  Shaped 
Spectrum  Signal,  Scenario  II.  The  magnitude  spectrum  of 
this  notched,  spectrally  encoded  TDCS  signal  displays  the  fre¬ 
quency  components  remaining  after  notching  the  independently 
distorted  interference  spectrum. 
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the  interference  and  communication  channels,  it  is  equally  likely  that  the  interferer’s 
peak  power  spectra  occur  in  bands  that  are  “amplified”  due  to  the  spectral  shaping. 
Again,  since  it  is  assumed  that  on  average,  the  channel  passes  through  100%  of  the 
energy  transmitted,  and  since  the  majority  of  frequency  spectrum  falls  below  a  mag¬ 
nitude  of  1,  it  is  predicted  that  the  un-notched,  spectrally  encoded,  shaped  spectrum 
TDCS  signals  will  have  an  overall  “filtering”  effect.  Therefore  the  simulated  Pi  for 
the  un-notched,  spectrally  encoded  TDCS  signals  is  predicted  to  be  better  (lower) 
than  the  analytical  prediction. 

The  TDCS  magnitude  spectrum  in  Figure  3.16  demonstrates  that  even  after 
notching,  the  spectrally  encoded  TDCS  signals  will  likely  yield  an  additional  “filter¬ 
ing”  effect  on  the  remaining  interference  energy.  Therefore  the  simulated  Pi  for  the 
notched,  spectrally  encoded  TDCS  signals  is  predicted  to  be  better  (lower)  than  the 
analytical  prediction. 

To  sum  up  the  second  scenario,  the  analytical  results  are  predicted  to  be  upper 
and  lower  bounds  of  the  simulation  results.  Specifically,  the  un-notched  analytical 
results  are  predicted  to  be  an  upper  bound,  and  the  notched  analytical  results  are 
predicted  to  be  a  lower  bound. 

3.7.3  Interference  Avoidance:  Interferer  Description.  For  both  interference 
avoidance  cases,  the  magnitude  spectrum  presented  in  [6,  pg4-4]  is  used  to  simulate 
a  narrow  band  interferer.  The  shape  of  the  interference  source’s  magnitude  spectrum 
is  shown  in  Figure  3.17.  For  reference,  the  spectral  magnitudes  are  included  in  Sec¬ 
tion  A.  10.  This  shape  is  interpolated  to  match  the  number  of  frequency  components 
used  in  this  research  using  the  Matlab®  interpft  function.  Additionally,  the  same 
notching  algorithm  discussed  in  Section  2.1.1  is  used  to  compare  the  un-notched  and 
notched  signals.  Where  the  interferer  spectral  magnitudes  exceed  the  manually  set 
threshold  of  1,  the  spectral  components  Ap  in  the  spectrally  encoded,  shaped  spec¬ 
trum  TDCS  signals  are  notched.  This  results  in  a  theoretical  reduction  of  11.8dB  in 
the  interferer’s  energy  in  the  TDCS  receiver  [6,  4-4], 
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Figure  3.17:  Magnitude  Spectrum  of  Conceptual  Interference 
Source 

3. 7-4  Interference  Avoidance:  Analytical  Results.  Figures  3.18  and  3.19  con¬ 
tain  the  analytical  predictions  of  the  interference  avoidance  capability  for  the  spec¬ 
trally  encoded,  shaped  spectrum  TDCS  signals  against  the  narrow  band  interferer 
described  in  Section  3.7.3.  The  un-notched  interference  PSD  N[  is  computed  using 
an  average  power  that  is  64  times  the  TDCS  transmitter  average  power.  The  notched 
interference  PSD  Ni  is  computed  using  an  average  power  that  is  11.8dB  lower  than 
the  un-notched  power  based  on  the  effect  of  the  manually  set  notching  threshold.  The 
Pb  for  the  notched  and  un-notched  signals  is  numerically  computed  by  solving  for  the 
expected  value  of  (3.71)  and  (3.68)  over  the  shaped  spectrum  TDCS  signal  energy  pdfs 
described  in  Section  3.5.3.  Each  curve  compares  the  analytically  predicted  Pb  for  both 
the  notched  and  un-notched  signals  for  each  level  of  L  independently.  As  mentioned 
in  the  sections  above,  these  curves  are  not  exact  predictors  of  performance  since  the 
effects  of  multipath  fading  on  the  interferer  and  the  interference  avoidance  algorithm 
have  not  been  quantified.  Therefore,  these  curves  are  the  predicted  bounding  cases 
for  the  simulation  results. 
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Figure  3.18:  Notched  vs  Un-notched  Analytical  Results,  L  —  2,  4.  These 
figures  demonstrate  the  predicted  interference  avoidance  capability  of  the  notching 
algorithm  applied  to  spectrally  encoded,  shaped  spectrum  TDCS  signals  using  an 
L-diversity  RAKE  receiver.  For  these  plots,  J/S  =  64  and  the  notching  algorithm 
attenuates  average  interferer  power  by  11.8dB 


Figure  3.19:  Notched  vs  Un-notched  Analytical  Results,  L  =  20,  50.  These 
figures  demonstrate  the  predicted  interference  avoidance  capability  of  the  notching 
algorithm  applied  to  spectrally  encoded,  shaped  spectrum  TDCS  signals  using  an 
L-diversity  RAKE  receiver.  For  these  plots,  J/S  =  64  and  the  notching  algorithm 
attenuates  average  interferer  power  by  lL8clB 


3. 8  Summary 

In  this  chapter,  the  analysis  of  the  TDCS  in  a  FSSF  multipath  channel  was 
presented.  The  results  of  the  first  portion  of  the  analysis  are  the  predicted  Pb  for 
spectrally  encoded,  shaped  spectrum  and  un-encoded,  flat  spectrum  TDCS  signals 
in  the  FSSF  channel.  The  second  section  of  the  analysis  qualitatively  predicts  the 
interference  avoidance  capabilities  of  the  spectrally  encoded,  shaped  spectrum  TDCS 
signals  with  and  without  notching.  Two  scenarios  for  interference  avoidance  were 
presented:  the  case  in  which  the  interferer  traverses  an  independent  FSSF  multipath 
channel,  and  the  case  in  which  the  interferer  traverses  an  AWGN  only  channel. 
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IV.  Simulation  Methodology 


Chapter  III  provides  the  mathematical  analysis  and  associated  results  supporting 
the  two  goals  of  this  research.  This  chapter  describes  the  simulation  methodol¬ 
ogy  used  to  gather  the  results  that  verify  the  analysis.  The  system  design  is  discussed 
first,  followed  by  the  design  of  experiments. 

4-1  System  Definition 

This  research  simulates  the  system  under  test  shown  in  Figure  4.1.  The  TDCS 
transmitter,  L-diversity  TDCS  RAKE  receiver,  interference  source,  independent  com¬ 
munications  and  interference  channels,  and  AWGN  source  are  all  within  the  system 
boundary.  The  TDCS  transmitter /RAKE  receiver  pair  are  modeled  as  described  in 
Figures  2.1  and  3.1.  The  communications  channel  is  modeled  as  a  FSSF  multipath 
channel  as  depicted  in  Figure  2.5.  The  interference  channel  is  either  characterized  as 
a  FSSF  multipath  channel  or  a  straight  pass  through.  The  interference  source  has 
the  spectral  shape  given  by  Figure  3.17.  The  AWGN  source  has  a  one-sided  power 
spectral  density  of  Na. 

4-1.1  System  Services.  Quite  simply,  the  service  provided  by  the  system 
under  test  is  symbol  detection  and  estimation.  The  two  possible  outcomes  of  the 
service  are  either: 

1.  Symbol  correctly  detected  and  estimated 


Distribution  <j  L  N 

of  c„ 

Figure  4.1:  System  Under  Test 
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2.  Symbol  incorrectly  detected  and  estimated 

Since  the  system  uses  BPSK  modulation,  symbols  are  mapped  directly  to  bits.  There¬ 
fore,  the  quality  of  the  system’s  symbol  detection  and  estimation  is  measured  by  the 
probability  of  bit  error,  R,  given  by  (4.1). 

Number  of  bit  errors 
b  Total  bits  transmitted 

4-1-2  Workload.  The  workload  for  the  system  is  the  ratio  Eb/N0.  The 
workload  levels  are  discussed  below. 

4-1-3  System  Parameters.  In  this  research,  parameters  are  the  aspects  of 
the  system  that  could  induce  variation  in  the  system  services.  Each  of  the  system 
parameters  for  this  research  are  defined  and  discussed  with  any  assumptions  in  the 
sections  below.  Parameters  that  are  varied  in  the  experiments  are  named  factors  and 
are  discussed  in  detail  in  a  later  section. 


(4.1) 


4-1.3. 1  TDCS  Transmitter /Receiver  Subsystem  Parameters. 

Modulation  Type  Cyclic  shift  keying  and  phase  shift  keying  are  two  primary  mod¬ 
ulation  schemes  historically  applied  to  the  TDCS  [6]  [9].  Since  both  have  the 
same  “noise-like”  characteristics,  yielding  a  low  autocorrelation,  this  parame¬ 
ter  should  have  little  effect  on  the  performance  of  the  system.  Binary  phase 
shift  keying  (BPSK)  is  used  exclusively  in  the  simulations  since  the  analytical 
performance  of  the  TDCS  is  accurately  predictable  using  traditional  analysis 
techniques. 

Data  Symbols  The  data  modulated  using  the  TDCS  affects  the  actual  signals  trans¬ 
mitted.  Since  this  research  assumes  equal  energy  antipodal  signaling,  the  data 
should  not  effect  the  energy  being  transmitted  or  probability  of  bit  error  perfor¬ 
mance.  Therefore,  uniformly  random  symbols  are  generated  and  transmitted. 
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Bit  rate,  Rb  This  represents  the  bit  rate  into  the  TDCS  modulator.  This  parameter 
along  with  the  modulation  scheme  determines  the  symbol  rate  and  duration. 
Since  BPSK  modulation  is  used  in  the  experiments,  Rb  =  fsb-  To  meet  the 
criteria  outlined  in  Section  2.2.2  required  for  the  channel  to  be  modeled  as  a 
FSSF,  Rb  is  set  at  10kbps. 

Fundamental  frequency,  fsb  This  is  the  lowest  frequency  component  in  the  TDCS 
waveform  as  well  as  the  tone  spacing  between  the  spectral  components  of  the 
waveform.  Since  it  is  also  equal  to  the  symbol/bit  rate  as  described  above,  fsb 
is  equal  to  10kHz. 

Number  of  spectral  components,  P  The  number  of  spectral  components  in  the 
TDCS  waveform  affects  the  bandwidth  of  the  transmitted  signal  W bb-  As  stated 
in  Chapter  II,  it  also  affects  the  channel  model.  To  meet  the  criteria  for  the 
channel  to  be  assumed  a  FSSF,  P  is  set  to  1024.  This  yields  a  baseband  band¬ 
width  of  greater  than  10MHz,  which  meets  the  criteria  discussed  in  Section  2.2.2 
for  both  urban  and  suburban  environments  in  all  cases. 

Spectral  amplitude  vector  Ap  The  spectral  amplitudes  of  the  TDCS  signals  are 
a  system  factor  and  are  discussed  later. 

Notching  Algorithm  This  parameter  effects  the  amount  of  energy  from  the  inter- 
ferer  that  is  present  in  the  demodulator  test  statistic.  The  notching  algorithm 
is  a  system  factor. 

Data  phase  encoding,  6^  .  Since  the  modulation  scheme  is  BPSK,  the  data  phase 
maps  data  bits  to  phase  shifts  according  to  the  following  scheme:  data  ‘0’  to  0 
radians  or  a  data  ‘1’  to  n  radians. 

Random  phase  encoding,  <fip  A  number  of  different  phase  encoding  schemes  are 
mentioned  in  the  literature.  To  minimize  the  autocorrelation  of  the  TDCS 
signals,  a  P  —  1  length  Gold  Code  {5%}  is  used  to  randomize  the  phase.  The 
phase  mapping  of  <f>p  to  either  0  or  n  radians  occurs  in  the  same  manner  as  the 
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data  bits  are  mapped  to  the  phase  of  Ok- 


4>p 


0  if  gp  —  0  Vp  7^  0 
<  7T  if  gp  —  1  Vp  7^  0 
0  if  p  =  0 


(4.2) 


No  multiple  access  phase  coding  is  used.  Since  the  autocorrelation  of  Gold  codes 
is  three- valued  regardless  of  the  code,  the  selection  of  generator  polynomials  and 
initial  states  is  not  expected  to  affect  the  results  of  the  simulations.  During  pilot 
simulations,  many  generator  polynomials  and  initial  states  were  used  without 
effect  on  Pb  performance.  For  repeatability,  generator  polynomials  2011  and 
2415  with  initial  states  1  and  71  respectively  are  used  to  generate  the  Gold 
Code  in  all  simulations. 

Digital-to-analog  converter  speed,  fsamp  This  is  the  speed  the  TDCS  transmit¬ 
ter  converts  digital  samples  to  generate  the  analog  waveform  for  transmission 
and  demodulation.  It  is  also  the  receiver’s  sampling  rate  of  the  downconverted 
(if  used)  waveform.  For  this  research,  fsamp  is  set  at  the  Nyquist  sampling  rate 

2  Pfsb. 

Normalized  bit  energy  to  noise  power  spectral  density  ratio,  Eb/N0  This  ra¬ 
tio  is  measured  with  reference  to  the  transmitted  Eb  and  is  distinct  from  the 
SNR  in  the  demodulator.  The  received  SNR  in  the  demodulator  7 b,  is  calculated 
after  passing  through  the  channel  model.  As  mentioned  above,  this  parameter 
is  the  system  workload. 

Shaping  The  algorithm  used  to  provide  the  spectral  shape  of  the  TDCS  signals. 
This  parameter  is  a  factor. 


4- 1.3. 2  Interference  Source  Parameters. 

Interferer  Spectrum,  Jp  The  spectral  shape  of  an  interference  source  effects  the 
number  of  frequency  components  that  are  notched  out  of  the  TDCS  spectrum. 
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To  ensure  consistency  with  results  previously  published,  and  to  limit  the  scope 
of  this  research,  the  interference  source  used  in  these  simulations  is  identical  to 
the  narrow-band  interferer  spectrum  used  in  [6,  4-3,4], 

Interferer  Power  to  Signal  Power  Ratio,  J/S  This  ratio  determines  the  power 
spectral  density  of  the  interference  source  in  the  test  statistic.  To  compare  the 
results  in  this  research  with  those  in  the  literature,  J/S  is  set  to  make  the 
interferer  power  spectral  density  identical  to  that  used  by  [6],  Using  (2.12),  and 
the  parameters  from  [6,  4-3,4],  Ni)Nunez  is  calculated  in  (4.3). 


N, 


i, Nunez 


Jo 

p 


J/S  x  E. 


sym 


2  E 


sym 


P 


32 


(4.3) 


To  solve  for  the  appropriate  J/S  in  these  experiments,  the  parameter  P  =  1024 
as  stated  above  is  used. 


EiXJaorw 

J/S  X  Egym 


1024 


J/S  X  Egym 


1024 


=  N„ 


i,  Nunez 


2  E 


sym 


32 


J/S  =  64 


(4.4) 

(4.5) 

(4.6) 


4- 1.3.3  Channel  Model  Parameters. 

Distribution  of  cn  Due  to  the  assumption  that  the  channel  model  is  Rayleigh  faded 
for  this  research,  the  cn  are  complex  Gaussian  processes  with  a  standard  devia¬ 
tion  of  a. 

Given  the  distribution  of  cn,  the  parameter  a  must  now  be  determined  based 
on  the  assumption  that  the  channel  returns  on  average  100%  of  the  energy 
transmitted  through  it.  Recalling  that  the  magnitudes  of  the  cn  are  Rayleigh 
random  variables,  the  first  and  second  moments  of  the  Rayleigh  distributed 
an  can  now  be  determined  [7,  148].  The  expression  V[-]  denotes  the  variance 
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operator. 


E  [a] 

■  fr 

(4.7) 

V[«] 

=  (2  —  7t/2  )a2,  therefore: 

(4.8) 

E[a2} 

=  V[a]  +  (E[a])2  =  2a2 

(4.9) 

The  previously  stated  assumption  that  on  average,  the  channel  returns  100%  of 
the  energy  input  by  the  transmitter,  and  the  statistics  above,  yield  the  expec¬ 
tation  in  (4.10)  when  the  an  are  independent  and  identically  distributed. 


r  Tsym 

r  Tsym 

L 

2 

/  Mt^dt 

=  1  =  E 

l 

ane3&n5{t  -  nTsamp ) 

dt 

0 

n= 1 

L 


1 


E 


n=l 


E  [£<*?.] 


1  =  L2a2,  therefore: 


(4.10) 


AWGN  Model  This  parameter  affects  the  environment  in  which  the  TDCS  receiver 
is  trying  to  detect  and  estimate  symbols  within.  Because  complex  low-pass 
equivalent  analytical  models  are  used  in  the  analysis,  the  noise  model  uses  a 
complex  Gaussian  noise  distribution  with  a  one-sided  power  spectral  density  of 
Na. 

Number  of  Taps  in  Multipath  Model,  L  This  parameter  is  related  to  T.m,  and 

f samp  by  (2.28).  The  ranges  of  values  for  Tm  and  L  are  determined  to  model  a 
suburban  environment  [10,  162],  The  specific  values  for  Tm  and  their  associated 
L  are  discussed  in  the  factor  section. 
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4-2  Experimental  Design 


Given  the  analysis  above,  the  system  workload  and  the  factors  selected  for  vari¬ 
ation  in  the  experiments,  the  design  of  experiments  are  presented  next.  Table  4.1 
enumerates  all  of  the  proposed  experiments.  In  each  experimental  batch,  two  types 
of  TDCS  signals  are  simulated  through  the  identical  conditions.  The  two  alterna¬ 
tives  compared  are  either  spectrally  encoded,  shaped  spectrum  TDCS  signals  vs.  un¬ 
encoded,  flat  spectrum  TDCS  signals;  or  notched  vs  un-notched  spectrally  encoded, 
shaped  spectrum  TDCS  signals. 


Table  4.1:  Table  of  Experiments 


Batch 

Eb/N0 

L 

Spectral 

Encoding 

Notching 

Int. 

Interference 

Channel 

1 

0-20dB,  2dB  steps 

2 

Flat /Shaped 

None 

N/A 

N/A 

2 

0-18dB,  2dB  steps 

4 

Flat /Shaped 

None 

N/A 

N/A 

3 

0-12dB,  2dB  steps 

20 

Flat /Shaped 

None 

N/A 

N/A 

4 

0-12dB,  2dB  steps 

50 

Flat /Shaped 

None 

N/A 

N/A 

5 

0-20dB,  2dB  steps 

2 

Shaped 

Notched/None 

NB1 

AWGN 

6 

0-18dB,  2dB  steps 

4 

Shaped 

Notched/None 

NB 

AWGN 

7 

0-12dB,  2dB  steps 

20 

Shaped 

Notched/None 

NB 

AWGN 

8 

0-1  OdB,  2dB  steps 

50 

Shaped 

Notched/None 

NB 

AWGN 

9 

0-20dB,  2dB  steps 

2 

Shaped 

Notched/None 

NB 

FSSF 

10 

0-20dB,  2dB  steps 

4 

Shaped 

Notched/None 

NB 

FSSF 

11 

0-12dB,  2dB  steps 

20 

Shaped 

Notched/None 

NB 

FSSF 

12 

0-1  OdB,  2dB  steps 

50 

Shaped 

Notched/None 

NB 

FSSF 

Harrow  Band 


4-2.1  Factors. 

Eb/N0  Depending  on  the  level  number  of  taps  in  the  channel  model,  this  ratio  is  varied 
from  either  0-12. OdB  or  0-8. OdB  in  2.0dB  steps.  These  ranges  cover  predicted 
Pb  levels  from  0.1  through  10”5.  The  step  size  provides  sufficient  resolution 
to  resolve  trends  and  distinguish  performance  between  the  TDCS  signals  being 
compared  in  the  various  experiments. 
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L ,  Tm  The  length  of  the  delay  spread,  and  accordingly  the  number  of  delay  taps  in 
the  model,  are  varied  to  cover  the  range  of  typical  and  worst-case  suburban 
environments  [10,  162],  Given  the  fsamp  (2 P  •  10kHz  =  20.48MHz)  established 
above,  Table  4.2  contains  the  modeled  delay  spread  values  and  the  corresponding 
number  of  delay  taps. 


Table  4.2:  Delay  Spread  vs  Number  of  Delay  Taps 


Delay  Spread  /is 

No.  Delay  Taps  (L) 

0.1 

2 

0.2 

4 

1 

20 

2.5 

50 

Spectral  Shaping  Two  different  spectral  shaping  methods  are  compared.  Flat  spec¬ 
trum  TDCS  signals  use  equal  amplitude  frequency  components  (except  the  DC 
and  Nyquist  frequencies)  as  discussed  in  Section  3.4.  Shaped  spectrum  TDCS 
signals  use  the  FSSF  multipath  channel  transfer  function  magnitude  response 
to  establish  the  TDCS  spectral  shape  as  discussed  in  Section  3.5.  In  both  cases, 
the  basic  spectral  shape  is  scaled  to  maintain  constant  energy  signaling. 

Interference  Avoidance  To  limit  the  scope  of  this  research,  and  compare  to  pre¬ 
viously  published  results  from  [6] ,  the  interference  avoidance  algorithm  notches 
spectral  components  above  a  manually  set  threshold  of  1,  as  discussed  in  Sec¬ 
tion  3.7.3.  The  experiments  either  use  the  algorithm  or  disable  it.  In  an  en¬ 
vironment  without  multipath  fading,  this  algorithm  attenuates  the  interferer’s 
energy  in  the  test  statistic  by  11.8dB. 

Interferer  Channel  To  compare  the  performance  of  the  of  the  TDCS  interference 
avoidance  capability  as  discussed  above,  the  interference  channel  is  modeled 
in  two  ways.  To  simulate  an  environment  where  the  TDCS  receiver  and  the 
interferer  are  geographically  separated,  the  channel  is  independent  but  modeled 
identically  to  the  communications  channel.  When  simulating  the  close  proximity 
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interference  avoidance  situation,  the  interference  channel  passes  the  interferer 

through  unchanged. 

Now  that  the  listing  of  experiments  and  the  levels  of  the  factors  to  be  varied 
have  been  presented,  each  class  of  experiment  is  discussed  in  the  context  of  the  overall 
goals  of  this  research. 

4-2.2  Spectral  Encoding  Experiments.  The  first  four  batches  of  experiments 
are  designed  to  verify  the  first  hypothesis:  that  spectral  shaping  should  achieve  a 
better  probability  of  bit  error  performance  through  a  FSSF  multipath  fading  channel. 
These  experiments  simulate  the  transmission,  detection  and  estimation  of  both  flat 
spectrum  and  shaped  spectrum  TDCS  signals  through  identical  channels.  No  inter¬ 
ferer  is  simulated  in  these  experiments.  For  each  level  of  L,  Eb/N0  is  varied  across  the 
entire  range  of  values  specified  for  Batches  1  through  4  in  Table  4.1.  The  probability 
of  bit  error  for  flat  and  shaped  spectrum  signals  is  compared  over  the  entire  range 
of  Eb/N0  independently  for  each  level  of  L.  The  simulated  results  are  also  verified 
against  the  analytical  predictions  using  numerical  computations  of  the  formulas  from 
Sections  3.4.1  and  3.5.2. 

4-2.3  Interference  Avoidance  Experiments.  The  last  eight  batches  of  ex¬ 
periments  are  designed  to  verify  the  second  hypothesis:  that  the  TDCS  retains  an 
interference  avoidance  capability  in  a  FSSF  multipath  environment.  These  experi¬ 
ments  simulate  the  transmission,  detection  and  estimation  of  shaped  spectrum  TDCS 
signals  in  the  presence  of  an  interferer  both  with  and  without  the  notching  algorithm. 
Again,  for  each  level  of  L,  Eb/N0  is  varied  across  the  entire  range  of  values  speci¬ 
fied  for  Batches  5  through  12  in  Table  4.1.  In  Batches  5  through  8,  the  interferer  is 
passed  through  an  independent  FSSF  multipath  channel.  In  Batches  9  through  12, 
the  interferer  is  passed  through  directly  to  the  AWGN  source  in  the  channel. 

The  probability  of  bit  error  for  the  notched  and  un-notched  signals  is  compared 
over  the  entire  range  of  Eb/N0  independently  for  each  level  of  L.  The  two  cases 
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of  simulation  results  are  considered  independently  and  compared  to  the  analytical 
predictions  using  numerical  computations  of  the  formulas  from  Equations  2.11. 

4-2-4  Confidence  Intervals  and  Iterations  Per  Experiment.  Confidence  in¬ 
tervals  on  the  simulation  results  provide  an  assurance  that  the  simulated  results  are 
within  a  quantifiable  statistical  range  of  accuracy.  Equation  (4.11)  is  the  confidence 
interval  on  any  ratio  p  [3,  217]. 

100(1  —  a)%  Confidence  Interval  =  p  ±  Zi_a/2]l^  -  — —  (4.11) 

In  (4.11),  zt  is  the  t^-quantile  of  the  standard  normal  distribution  and  n  is  the  number 
of  iterations  in  the  experiment.  It  is  desired  that  the  confidence  interval  reach  a  certain 
precision  to  quantify  the  confidence  interval.  The  relationship  between  the  confidence 
interval  and  the  precision  i/j,  in  terms  of  the  ratio  p  is  (4.12). 

p±fip  =  p  ±  Zi-a/2^ n  ^  (4.12) 

Therefore,  the  number  of  iterations  required  for  an  experiment  to  achieve  a  100(l-a) % 
confidence  interval  that  is  within  of  the  analytically  predicted  probability  of  bit  error 
Pbl  is  found  by  manipulating  (4.11)  into  the  following  form. 

n  =  (4-13) 

Using  (4.13),  the  number  of  iterations  required  per  experiment  is  calculated  using  the 
confidence  intervals  shown  in  Table  4.3  and  the  analytical  predictions  discussed  for 
each  experiment  in  Sections  4.2.2  and  4.2.3.  The  differences  between  the  batches  are 


Table  4.3:  Confidence  Intervals  by  Experimental  Batch 


Batch 

Precision 

Confidence  Interval  (l-a:)100% 

1,  2,  5,  6,  9,  10 

10%  Pb 

95% 

3,  4,  7,  8,  11  12 

10%  Pb 

90% 
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Code  Header  (Set  constants,  compute  gold  code,  compute  iterations  required) 

^.Iteration  Loop  -  until  max  iterations  reached: 

Generate  comm  channel  impulse  response  and  transfer  function 

Generate  interferer  channel  impulse  response  and  transfer  function  (if  required) 

Generate  FSWs  for  both  signal  sets  being  compared  (Flat  vs  Shaped/Notched  vs  Un-Notched) 

Generate  block  of  random  bits  and  modulate  using  both  sets  of  FSWs 
Generate  interfering  signal(if  required) 

Process  both  sets  of  comm  signals  and  interferer  (if  required)  through  independent  channels 
^Eb/No  Loop  -  if  iteration  >  100,  for  only  those  levels  still  required  to  converge 
Generate  noise  according  to  Eb/No  level 
Demodulate  both  signal  sets  using  comm  +  interferer  +  noise 
Compare  both  sets  of  demodulated  bits  to  original  modulated  bits 

_ Count  bit  errors  for  both  signal  sets 

- Save  interim  simulation  results 

Code  Footer  (Save  final  simulation  results) 

Figure  4.2:  Structure  of  Monte  Carlo  Simulation 

to  compensate  for  the  longer  simulation  times  required  for  experiments  when  L  =  20, 
50.  Since  these  simulations  run  longer,  a  lower  confidence  interval  is  designed  to 
decrease  the  simulation  run  time,  yet  still  provide  a  high  degree  of  confidence  that 
the  simulated  results  are  within  the  specified  precision. 

4-2.5  Simulation  Structure.  All  experiments  are  simulated  using  Monte 
Carlo  methods  to  average  out  the  impacts  of  the  random  processes  on  the  results.  The 
components  in  the  system  under  test  are  simulated  using  complex  low-pass  equivalent 
models.  The  simulation  is  structured  to  evaluate  several  experiments  at  one  time.  As 
shown  in  Figure  4.2,  within  a  given  batch,  all  experiments  required  for  a  single  level  of 
L  (full  range  of  Eb/N0,  both  TDCS  signal  types)  are  simulated  in  the  same  simulation 
code. 


Blocks  of  random  bits  modulated  inside  the  iteration  loop  are  100  bits  long. 
Therefore  the  number  of  iterations  is  the  maximum  required  according  to  the  confi¬ 
dence  interval  calculations  divided  by  the  block  size.  For  a  single  channel  instanti¬ 
ation,  both  shaping  or  interference  avoidance  techniques  are  simulated.  Both  signal 
sets  are  processed  through  identical  channel  instantiations,  interferers  and  noise  in- 
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stantiations.  This  consistency  ensures  that  the  effect  being  measured  is  due  only  to 
the  difference  in  the  FSW  structure.  The  logic  check  in  the  E\,/N0  loop  requires  that 
all  experiments  are  simulated  through  at  least  100  channel  instantiations.  At  lower 
Eb/Na  levels,  the  number  of  bits  required  for  the  simulations  to  converge  to  the  re¬ 
quired  confidence  intervals  is  on  the  order  of  100s  of  bits.  Without  the  logic  check, 
the  simulation  would  terminate  this  experiment  after  only  a  few  instantiations  of  the 
FSSF  multipath  channel,  which  would  not  average  out  transient  effects  in  the  channel. 

4 . 3  Summary 

This  chapter  contains  the  foundation  for  the  simulation  portion  of  this  research. 
The  system  under  test  is  described  along  with  the  parameters  that  could  affect  the 
results  of  the  simulation.  The  experimental  design,  including  the  system  workload 
and  the  factors  to  be  varied  during  the  simulation  experiments  are  also  described. 
Finally,  the  methods  for  verifying  the  analysis  through  simulation  were  described  and 
tied  to  the  goals  and  hypothesis  of  this  research. 
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V.  Simulation  Results 


Earlier  chapters  review  the  literature  survey,  analysis  and  design  of  experiments. 

This  chapter  contains  the  results  of  the  simulation  experiments.  The  first  sec¬ 
tion  compares  the  simulation  results  of  spectrally  encoded,  shaped  spectrum  TDCS 
signals  to  un-encoded,  flat  spectrum  TDCS  signals.  The  second  section  compares  the 
simulation  results  of  notched  and  un-notched,  spectrally  encoded,  shaped  spectrum 
TDCS  signals  in  the  two  interference  avoidance  scenarios  described  in  Section  3.7. 
Both  sections  compare  the  results  of  two  types  of  simulations  to  the  analytical  pre¬ 
dictions  contained  in  Chapter  III.  Each  section  contains  independent  plots  for  each 
value  L  (2,  4,  20,  50). 

5.1  Spectral  Encoding  Results 

As  discussed  in  Chapter  IV,  the  first  set  of  experiments  are  designed  to  answer 
the  first  goal  of  the  research,  namely:  to  compare  the  performance  of  spectrally 
encoded,  shaped  spectrum  TDCS  signals  to  un-encoded,  flat  spectrum  TDCS  signals 
in  a  FSSF  multipath  channel  using  an  L-diversity  TDCS  RAKE  receiver.  This  section 
presents  p,  curves  for  each  level  of  L  containing  the  simulation  results  and  analytical 
predictions.  As  a  reminder,  in  each  figure,  probability  of  bit  error  is  plotted  against 
the  transmitted  Eb/Na.  Each  Pb  curve  contains  two  reference  lines,  two  analytical  Pi, 
prediction  lines  and  two  simulation  P b  lines.  The  first  reference  line  represents  a  BPSK 
matched  filter  demodulator  through  a  FSSF  multipath  channel  with  no  diversity  (i.e., 
worst-case).  The  second  reference  is  the  BPSK  matched  filter  demodulator  through  an 
AWGN  channel,  representing  the  theoretical  best-case.  The  analytical  Pb  predictions 
are  identical  to  those  plotted  in  Chapter  III. 

5.1.1  Spectral  Encoding:  L  =  2.  The  first  two  sections  of  the  results  discuss 
the  simulations  of  typical  suburban  environments  [10,  162],  Figure  5.1  shows  the  sim¬ 
ulated  results  in  comparison  to  the  two  reference  lines  and  the  two  predictions.  By 
visual  inspection,  it  is  clear  that  each  Pb  prediction  is  consistent  with  the  simulated 
results.  Qualitatively,  the  simulated  spectrally  encoded,  shaped  spectrum  TDCS  sig- 
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Figure  5.1:  Spectral  Encoding  Results,  L  —  2.  The  simulations  are  consistent 
with  the  analysis.  The  results  show  approximately  l.OdB  gain  realized  by  spectral 
encoding. 
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nal  has  a  better  Pb  than  the  un-encoded,  flat  spectrum  TDCS  signal  at  each  Eb/N0  as 
predicted.  Using  (3.67),  and  by  graphically  estimating  gain  G  in  Figure  5.1,  the  spec¬ 
trally  encoded  TDCS  signals  provide  an  approximate  ldB  gain  over  the  un-encoded 
TDCS  signals  as  seen  in  the  analytical  results. 

After  simulating  3.2  x  10'  bits,  it  is  notable  that  the  analytical  results  are  not 
within  10%  of  the  simulation  results  as  predicted  by  the  95%  confidence  interval. 
The  difference  between  the  analytical  and  simulation  results  may  be  explained  by 
inexact  Pb  predictions.  The  spectrally  encoded  signal  energy  pdf  estimation  may  cause 
inaccuracies  in  the  numerical  Pb  calculations.  Errors  in  the  analytical  Pb  prediction 
approximations  may  be  the  source  of  differences  for  the  un-encoded  TDCS  signals. 

5.1.2  Spectral  Encoding:  L  =  4.  The  Pb  curve  for  L  =  4  in  Figure  5.2 
reveals  several  qualitative  results.  However,  since  the  simulation  did  not  complete, 
some  general  comments  are  in  order.  The  simulation  generated  4.0  x  10'  bits  which 
is  sufficient  for  the  confidence  interval  to  converge  as  designed  (c.f.  Section  4.2.4) 
for  spectrally  encoded,  shaped  spectrum  TDCS  signal  Eb/N0  between  O.OdB  and 
12.0dB,  and  for  un-encoded,  flat  spectrum  TDCS  signal  Eb/N0  between  O.OdB  and 
16.0dB.  Additional  iterations  may  have  caused  the  simulated  results  to  converge  to  the 
analytical  predictions;  however,  time  was  not  available  to  complete  the  simulations. 
In  light  of  this  information,  it  is  not  clear  whether  the  divergence  in  the  shaped 
signal  analytical  and  simulation  results  seen  in  Figure  5.2  is  due  to  inaccuracies  in 
the  analysis  or  whether  running  the  simulation  longer  would  cause  the  two  curves  to 
converge. 

Overall,  the  spectrally  encoded  TDCS  signals  perform  at  a  lower  Pb  than  the 
un-encoded  signals  as  predicted.  For  the  L  —  4  simulation  results,  the  graphically 
estimated  gain  for  spectrally  encoded  vs  un-encoded  signals  is  approximately  2dB 
across  the  range  of  Pb  plotted.  There  is  one  other  trend  worth  noting  in  these  re¬ 
sults.  At  lower  levels  of  Eb/N0,  notice  that  the  spectrally  encoded,  shaped  spectrum 
TDCS  signal  Pb  is  better  (less)  than  the  reference  line  for  the  BPSK  matched  fil- 
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Figure  5.2:  Spectral  Encoding  Results,  L  —  4.  Analytical  and  simulated 
results  are  consistent  again,  demonstrating  approximately  2.0dB  gain  for  spectral 
encoding.  Divergence  between  spectrally  encoded  simulation  and  analysis  results 
for  Eb/N0  >  12.0dB  caused  by  either  1)  incomplete  simulation,  or  2)  inaccurate 
analytical  prediction.  Note  that  spectral  encoding  P j,  is  better  than  the  “best-case” 
AWGN  reference  for  Eb/Na  <  4.0dB. 


ter  in  AWGN  only.  This  effect  is  caused  by  the  spectral  encoding  algorithm.  The 
algorithm  maximizes  “amplification”  and  minimizes  “filtering”  across  the  different 
frequency  bands  in  the  channel  transfer  function  spectrum.  Therefore  the  received 
SNR  in  the  demodulator  ^  is  greater  than  the  transmitted  SNR  Eb/N0.  This  phe¬ 
nomenon  combined  with  the  diversity  from  the  RAKE  receiver  is  what  enables  the 
spectrally  encoded,  shaped  spectrum  TDCS  signal  to  have  a  better  (lower)  Pb  than 
the  “best-case”  reference  line. 

5.1.3  Spectral  Encoding:  L  =  20.  As  a  reminder,  the  results  for  L  =  20,  and 
50  are  representative  of  the  worst-case  suburban  environment  described  in  [10,  162], 
For  L  =  20,  the  simulation  results  are  displayed  in  Figure  5.3. 

The  analytical  results  are  qualitatively  consistent  with  the  simulation  results 
for  both  the  spectrally  encoded,  shaped  spectrum  TDCS  signals  and  un-encoded,  flat 
spectrum  TDCS  signals.  The  spectrally  encoded  signals  also  perform  at  a  lower  Pb 
than  the  un-encoded  signals  as  predicted.  In  the  L  =  20  simulations,  the  spectrally 
encoded  vs  un-encoded  graphical  estimation  of  gain  is  approximately  2.5dB  across  the 
range  of  Pb  plotted.  Note  that  across  the  range  of  transmitted  Eb/N0  in  these  results, 
the  spectrally  encoded,  shaped  spectrum  TDCS  signal  Pb  is  better  (lower)  than  the 
BPSK  in  AWGN  reference.  As  discussed  in  the  L  =  4  section,  this  is  a  result  of 
the  combined  effects  of  the  spectral  encoding  algorithm  and  receiver  diversity.  The 
un-encoded  Pb  curve  appears  to  be  converging  to  the  BPSK  in  AWGN  reference  line 
as  L  increases.  Since  the  research  assumptions  and  simulation  design  ensure  that 
on  average,  100%  of  the  energy  transmitted  into  the  channel  is  returned,  then  as 
the  diversity  increases,  it  is  postulated  that  the  un-encoded,  flat  spectrum  TDCS 
signals  will  approach  the  BPSK  in  AWGN  Pb  reference  until  the  intersymbol-  and 
self- interference  terms  increase  and  overcome  the  benefit  of  the  diversity. 

5.1.4  Spectral  Encoding:  L  =  50.  For  L  =  50,  Figure  5.4  contains  the 
analysis  and  simulation  results.  As  in  the  previous  comparisons  of  the  spectrally  en¬ 
coded,  shaped  spectrum  TDCS  signals  and  un-encoded,  flat  spectrum  TDCS  signals, 
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Figure  5.3:  Spectral  Encoding  Results,  L  =  20.  Analytical  and  simulated  re¬ 
sults  are  consistent  and  demonstrate  approximately  2.5dB  gain  for  spectral  encoding. 
Spectral  encoding  P&  is  better  than  the  “best-case”  AWGN  reference  for  all  Eh / N0  in 
plot.  Un-encoded  results  start  to  converge  toward  AWGN  reference. 


Figure  5.4:  Spectral  Encoding  Results,  L  =  50.  Analytical  and  simulated 
results  are  consistent  and  demonstrate  approximately  2.75dB  gain  for  spectral  encod¬ 
ing.  Again,  spectral  encoding  P&  is  better  than  the  “best-case”  AWGN  reference  for 
all  Eb/Na  in  plot.  Un-encoded  results  continue  to  converge  toward  AWGN  reference. 


the  analytical  results  are  consistent  with  the  simulation  results.  The  un-encoded  Pb 
curves  continue  to  converge  toward  the  BPSK  in  AWGN  reference  and  the  spectrally 
encoded  P&  curves  continue  to  improve  beyond  the  BPSK  in  AWGN  reference.  In 
the  L  —  50  simulations,  the  graphically  estimated  gain  for  spectrally  encoded  vs 
un-encoded  signals  is  slightly  above  2.75dB  across  the  range  of  Pb  plotted. 

5.1.5  Spectral  Encoding  Results:  Interim  Summary.  In  the  sections  above, 
the  first  goal  of  the  research  has  been  achieved.  Spectrally  encoded,  shaped  spec¬ 
trum  TDCS  signal  Pb  has  been  compared  to  un-encoded,  flat  spectrum  TDCS  signal 
Pb  through  a  FSSF  multipath  channel  using  an  L-cliversity  RAKE  Receiver.  The 
hypothesis  for  this  portion  of  the  research  was  that  the  spectrally  encoded  signals 
would  pass  more  energy  through  the  FSSF  multipath  channel  and  therefore  have  a 
lower  (better)  Pb  than  the  un-encoded  signals.  Chapter  III  presents  the  analytical 
Pb  results  for  both  the  spectrally  encoded  and  un-encoded  signals.  Here,  the  results 
of  the  both  the  analysis  and  simulations  are  presented  with  the  overall  result  that 
the  hypothesis  was  correct.  A  range  of  l.OdB  to  2.75dB  gain  in  transmitted  Eb/N0, 
depending  on  the  amount  of  diversity  L  in  the  channel  and  TDCS  RAKE  receiver, 
is  realized  in  both  the  analyses  and  simulations  simply  by  applying  spectrally  en¬ 
coding  to  the  TDCS  signals  to  match  the  transfer  function  of  the  FSSF  multipath 
channel.  Additionally,  the  simulation  results  verify  the  analytically  predicted  Pb  for 
both  spectrally  encoded  and  un-encoded  TDCS  signals  using  an  L-diversity  RAKE 
receiver. 

5.2  Interference  Avoidance:  Spectrally  Encoded  Signals  With  and  With¬ 
out  Spectral  Notching 

The  results  discussed  in  this  section  address  the  second  goal  of  this  research:  to 
verify  the  interference  avoidance  capability  of  the  spectrally  encoded,  shaped  spectrum 
TDCS  signals  in  a  multipath  fading  environment  by  comparing  the  performance  of 
notched  and  un-notched  TDCS  signals.  As  before,  all  simulations  use  the  L-diversity 
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Figure  5.5:  Interference  Avoidance  I,  L  —  2.  No  Multipath  on  Interferer, 
Notched  vs  Un-Notched  Results.  The  notching  algorithm  shows  no  degradation  in 
performance.  The  spectrally  encoded  TDCS  signals  even  without  notching  demon¬ 
strate  a  significant  amount  of  unintentional  interference  avoidance. 

TDCS  RAKE  receiver.  In  Chapter  III,  two  scenarios  are  presented  and  analyzed. 
The  first  scenario  simulates  the  case  where  the  interferer  and  the  TDCS  receiver  are 
located  in  close  proximity  and  therefore  the  interference  channel  is  a  pass  through  only. 
The  second  scenario  simulates  the  situation  in  which  the  communications  channel  and 
interference  channel  are  independent  FSSF  multipath  channels.  As  in  the  previous 
section,  the  results  are  discussed  in  order  of  increasing  diversity  (L  =  2,  4,  20,  and  50) 
beginning  with  the  first  scenario.  To  simplify  the  plots,  only  the  analytical  predictions 
and  simulated  results  are  plotted.  The  reference  curves  are  omitted. 

5.2.1  Interference  Avoidance  I:  No  Multipath  on  Interference  Channel. 

5. 2. 1.1  Interference  Avoidance  I:  L  =  2.  For  L  =  2,  the  simulation 
results  are  displayed  in  Figure  5.5.  In  this  first  set  of  interference  avoidance  results, 
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the  general  trends  are  as  expected.  At  lower  transmitted  Eb/N0  ratios,  the  notched 
and  un-notched  Pb  performance  levels  are  nearly  identical.  Recalling  the  Pb  prediction 
relationship  from  (2.11),  it  is  clear  that  when  Na  is  large  compared  to  the  effect  of 
the  interference  source  is  imperceptible  on  Pb  performance.  At  higher  Eb/N0  ratios, 
the  opposite  effect  is  witnessed  in  the  un-notched  results  and  predictions.  In  this  case, 
the  effect  of  overshadows  Na  and  the  Pb  curve  begins  to  flatten  out.  However,  as 
expected,  the  notched  results  in  this  simulations  do  not  exhibit  this  flattening  effect 
at  lower  Eb/N0,  since  a  large  amount  of  the  interferer’s  energy  has  been  notched  out 
by  the  TDCS. 

In  Chapter  III,  it  was  predicted  that  the  analytical  curves  would  provide  upper 
bounds  for  both  the  notched  and  un-notched  simulations  results.  This  trend  is  evi¬ 
dent  in  both  the  notched  and  un-notched  simulations.  In  the  un-notched,  spectrally 
encoded  signals,  there  is  a  pronounced  difference  between  the  analysis  and  simulation. 
This  demonstrates  the  un-intentional  interference  avoidance  effect  predicted  in  Chap¬ 
ter  III.  However,  in  the  case  of  the  notched,  spectrally  encoded  signals,  the  separation 
is  small  between  the  analytical  predictions  and  simulation  results.  It  is  expected  that 
the  reason  for  this  limited  improvement  may  lie  in  the  relative  flatness  of  the  channel 
transfer  functions  when  L  is  small.  When  L  =  2,  the  channel  transfer  function  is 
relatively  flat  due  to  the  number  of  tap  delays  in  the  channel  model.  Therefore,  it  is 
likely  that  the  interferer’s  spectrum  is  changed  very  little  by  the  channel.  This  would 
cause  the  notching  algorithm  to  perform  nearly  as  expected. 

To  summarize,  for  L  =  2  the  simulation  results  are  slightly  better  than  the 
analytical  predictions,  as  expected  for  the  case  when  the  interferer  does  not  traverse 
a  multipath  fading  channel. 

5.2. 1.2  Interference  Avoidance  I:  L  =  4.  For  L  —  4,  the  simulation 
results  are  displayed  in  Figure  5.6.  The  trends  in  this  set  of  simulations  are  nearly 
identical  to  those  observed  for  the  L  =  2  case.  The  notched,  spectrally  encoded, 
shaped  spectrum  TDCS  signal  simulations  are  slightly  less  (better)  than  the  analyt- 
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Figure  5.6:  Interference  Avoidance  I,  L  —  4.  No  Multipath  on  Interferer, 
Notched  vs  Un-Notched  Results.  The  notching  algorithm  shows  no  degradation  in 
performance.  No  conclusions  are  drawn  for  the  notched  signals  above  Eb/N0  >  12.0dB 
due  to  inaccuracies  in  the  analytical  Pb  results.  The  spectrally  encoded  TDCS  sig¬ 
nals  even  without  notching  continue  to  demonstrate  the  unintentional  interference 
avoidance  effect. 

ical  Pb  curves  as  predicted  in  Chapter  III.  However  at  Eb/N0  >  14.0dB,  as  in  the 
spectral  encoding  results,  the  analytical  predictions  may  not  be  accurate.  The  un¬ 
notched,  spectrally  encoded  signals  continue  to  exhibit  both  the  analytically  predicted 
flattening  and  the  un-intentional  interference  avoidance. 

5. 2. 1.3  Interference  Avoidance  I:  L  =  20.  For  L  =  20  and  50,  the 
simulations  model  the  worst-case  suburban  delay  spread  characteristics  in  the  chan¬ 
nel.  For  L  =  20,  the  simulation  results  are  displayed  in  Figure  5.7.  This  set  of  results 
presents  some  unexpected  behaviors.  Even  though  the  interference  channel  is  a  pass 
through,  and  the  notching  algorithm  should  remove  the  full  11.8dB  of  interference 
energy,  the  Pb  for  the  notched  simulations  is  greater  (worse)  than  predicted.  It  is 
expected  that  due  to  the  high  number  of  taps  in  the  channel  model,  that  the  “fil- 
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Figure  5.7:  Simulated  Interference  Avoidance,  No  Multipath 
on  Interferer,  Notched  vs  Un-Notched  Results,  L=20 


tering”  in  the  remaining  portions  of  the  TDCS  spectrum  should  remove  more  than 
the  11.8dB  of  interferer  power  expected.  The  unexpected  results  may  be  due  to  the 
equally  likely  “amplifying”  effect  in  the  channel,  or  the  higher  PSD  caused  by  the 
notching  algorithm  which  effectively  reduces  the  bandwidth  of  the  TDCS  signal.  In 
contrast,  the  un-notched,  spectrally  encoded  TDCS  signals  continue  to  exhibit  the 
predicted  performance.  The  simulations  continue  to  exhibit  the  un-intentional  inter¬ 
ference  avoidance  and  flattening  due  to  the  interferer  power. 

5. 2. 1.4  Interference  Avoidance  I:  L  =  50.  For  L  =  50,  the  simulation 
results  are  displayed  in  Figure  5.8.  The  results  for  the  L  =  50  case  are  essentially  the 
same  as  those  seen  in  the  L  —  20  experiments.  The  notched  signals  present  the  same 
unpredicted  behavior,  performing  worse  than  predicted  while  the  un-notched  signals 
perform  as  expected.  The  exact  analytical  explanation  of  this  phenomena  is  left  for 
further  research. 
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Figure  5.8:  Simulated  Interference  Avoidance,  No  Multipath 
on  Interferer,  Notched  vs  Un-Notched  Results,  L=50 


5.2.2  Interference  Avoidance  II:  Independent  Multipath  Channels.  In  this 
second  set  of  interference  avoidance  simulations,  the  results  are  also  presented  in 
order  of  increasing  diversity  (L).  As  in  the  previous  interference  avoidance  results,  no 
references  are  included.  The  predictions  of  the  spectrally  encoded,  shaped  spectrum 
TDCS  signal  Pi  are  included  with  and  without  the  notching  algorithm  to  suppress 
the  interferer. 

5.2.2. 1  Interference  Avoidance  II:  L  =  2.  For  L  —  2,  the  simulation 
results  are  displayed  in  Figure  5.9.  The  results  for  the  second  interference  avoidance 
scenario  for  L  —  2  are  nearly  identical  the  first  scenario.  In  this  plot,  notice  that 
at  lower  Ei/N0  ratios,  the  notched  and  un-notched  Pi  performance  levels  are  nearly 
identical.  The  results  also  display  the  predicted  flattening  effect  on  the  un-notched 
TDCS  signal  curves  while  the  notched  results  in  this  simulations  do  not  exhibit  this 
flattening  effect. 
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Figure  5.9:  Simulated  Interference  Avoidance,  Independent 
Multipath  Channels,  Notched  vs  Un-Notched,  L=2 

In  Chapter  III,  it  was  predicted  that  the  analytical  predictions  would  be  upper 
and  lower  bounds  to  the  simulation  results.  While  the  un-notched  prediction  does 
provide  an  upper  bound  as  expected,  the  notched  simulations  appear  to  have  P, & 
equal  to  or  better  than  the  analytical  prediction.  As  before,  relying  on  intuition, 
it  is  expected  that  the  reason  for  this  improvement  may  lie  in  the  relative  flatness 
of  the  channel  transfer  functions  when  L  is  small.  It  is  likely  that  the  interferer’s 
spectrum  is  changed  very  little  by  the  channel,  causing  the  notching  algorithm  to 
perform  nearly  as  expected.  Therefore,  the  results  for  both  interference  avoidance 
scenarios  are  expected  to  nearly  identical. 

5. 2. 2. 2  Interference  Avoidance  II:  L  =  4.  For  L  —  4,  the  simulation 
results  are  displayed  in  Figure  5.10.  As  seen  in  the  L  —  2  results,  notice  that  at  lower 
Eb/N0  ratios,  the  notched  and  un-notched  P&  performance  levels  are  nearly  identical 
as  expected.  Also,  notice  the  flattening  effect  on  the  un-notched  TDCS  signal  curves 
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Figure  5.10:  Simulated  Interference  Avoidance,  Independent 
Multipath  Channels,  Notched  vs  Un-Notched,  L=4 

while  the  notched  results  in  this  simulations  do  not  exhibit  this  flattening  effect.  Again 
all  of  these  effects  are  expected. 

Returning  to  the  predictions  of  Chapter  III,  the  un-notched  prediction  does 
provide  an  upper  bound  as  expected.  In  contrast  to  the  predictions,  these  simulations 
show  the  notched  signals  appear  to  have  Pb  very  close  to  or  just  slightly  better  than 
the  analytical  predictions.  As  in  the  previous  case,  this  is  likely  due  to  the  relative 
flatness  of  the  channel  transfer  functions  for  L  =  4.  Recall  also  that  the  analytical  Pb 
predictions  are  unreliable  for  Eb/N0  >  14.0dB.  Notice  that  the  un-notched  simulations 
continue  to  perform  better  than  the  analytical  Pb  results  as  predicted,  but  not  as 
well  as  in  the  first  interference  avoidance  scenario,  as  shown  in  Figure  5.6.  The 
degraded  effect  on  the  un-intentional  interference  avoidance  may  be  caused  by  the 
second  independent  multipath  channel.  The  effects  of  the  shaping  from  the  spectral 
encoding  algorithm  may  counteract  the  “filtering”  from  the  independent  interference 
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Figure  5.11:  Simulated  Interference  Avoidance,  Independent 
Multipath  Channels,  Notched  vs  Un-Notched,  L=20 


channel.  This  would  cause  the  simulation  results  to  more  closely  match  the  analytical 
Pb  results  for  the  un-notched  signals. 

5. 2. 2. 3  Interference  Avoidance  II:  L  =  20.  As  a  reminder,  for  the 
simulations  in  which  L  —  20  and  50,  the  channels  are  modelling  worst-case  suburban 
environment  multi-path  delays.  For  L  =  20,  the  simulation  results  are  displayed 
in  Figure  5.11.  In  this  plot,  the  flattening  trend  exhibited  in  the  un-notched  signal 
causes  the  simulated  Pb  to  be  greater  (worse)  than  the  analytical  results.  This  early 
flattening  may  be  caused  by  the  combined  effects  of  the  independent  communication 
and  interference  channels  as  described  in  the  L  —  4  results  above. 

The  predictions  of  notched,  spectrally  encoded  Pb  from  Chapter  III  appear  to 
be  more  accurate  for  this  simulation.  The  notched  prediction  forms  a  lower  bound  as 
expected,  however  the  degradation  appears  to  be  slight.  To  summarize  the  results  for 
L  =  20,  the  notching  algorithm  does  not  appear  to  affect  the  full  11.8dB  attenuation 
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Figure  5.12:  Simulated  Interference  Avoidance,  Independent 
Multipath  Channels,  Notched  vs  Un-Notched,  L=50 

on  the  interferer  when  the  communications  and  interference  channels  are  indepen¬ 
dent  multipath  channels,  and  the  un-notched  results  do  not  exhibit  the  predicted 
un-intentional  interference  avoidance  for  Eb/N0  >  8.0dB. 

5. 2. 2. 4  Interference  Avoidance  II:  L  =  50.  For  L  =  50,  the  simulation 
results  are  displayed  in  Figure  5.12.  The  results  from  the  L  —  50  simulations  follow 
similar  trends  as  the  L  —  20  simulations,  yielding  a  similar  summary.  Since  the 
simulations  cover  Eb/N0  up  to  8.0dB,  the  flattening  effect  is  not  yet  realized  in  the 
un-notched,  spectrally  encoded  signals.  The  un-notched  signals  exhibit  the  predicted 
un-intentional  interference  avoidance  predicted.  The  interference  avoidance  notching 
algorithm  remains  effective  at  reducing  the  interferer’s  energy  by  slightly  less  than 
the  predicted  11.8dB. 
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5. 3  Summary 

In  this  chapter,  the  simulation  results  verify  most  of  the  predicted  outcomes  de¬ 
veloped  in  the  analysis  and  verify  the  hypotheses  set  forth  in  Chapter  II.  First  of  all, 
the  spectral  encoding  capability  of  the  TDCS  improves  performance  in  a  FSSF  mul¬ 
tipath  channel  over  un-encoded  signals  by  as  much  as  2.75dB  in  transmitted  Eb/N0. 
Second,  in  the  interference  avoidance  experiments,  for  both  the  independent  multipath 
case  and  the  case  in  which  the  interferer  experiences  no  multipath  fading,  the  notching 
algorithm  yields  an  improvement  in  probability  of  bit  error.  While  the  improvement  is 
not  as  much  as  predicted  for  strictly  AWGN  channels,  it  is  still  significant.  It  was  also 
demonstrated  that  even  without  notching  to  avoid  interference,  the  shaped  spectrum 
TDCS  signals  provide  an  distinguishable  interference  avoidance  capability. 
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VI.  Conclusions 


Having  completed  the  course  of  the  research,  this  chapter  reviews  the  thesis, 
summarizes  the  research  conclusions  and  discusses  areas  for  future  research. 

6. 1  Summary 

Before  reviewing  the  summary  of  this  work,  recall  the  goals  of  this  research  are: 

1.  To  compare  the  performance  of  spectrally  encoded,  shaped  spectrum  TDCS 
signals  to  un-encoded,  flat  spectrum  TDCS  signals  in  a  FSSF  multipath  channel 
using  an  L-diversity  TDCS  RAKE  receiver. 

2.  To  evaluate  the  interference  rejection  capability  of  an  L-diversity  TDCS  RAKE 
receiver  using  spectrally  encoded,  shaped  spectrum  TDCS  signals  in  a  FSSF, 
multipath  channel. 

To  this  end,  beginning  in  Chapter  II,  this  thesis  presents  a  brief  introduction 
to  the  M-ary  phase  shift  keyed  TDCS  communications  system,  the  associated  inter¬ 
ference  avoidance  capability,  and  probability  of  bit  error  prediction  in  AWGN.  An 
introduction  to  the  multipath  fading  channel  describes  basic  channel  parameters  and 
their  effect  on  the  channel  model.  To  support  the  assumptions  of  the  research,  the 
literature  search  also  includes  typical  900MHz  cellular  phone  multipath  channel  char¬ 
acteristics  and  a  brief  discussion  of  channel  estimation.  The  research  assumptions 
specifically  focused  on  the  frequency-selective,  slowly-fading  (FSSF)  multipath  chan¬ 
nel,  modeled  in  this  work  as  a  tapped  delay  line.  A  discussion  of  prediction  in  the 
multipath  environment  lends  insight  into  the  severe  performance  degradation  caused 
by  multipath  fading  channels.  To  overcome  this  degradation,  diversity  is  incorpo¬ 
rated  into  the  communication  system  using  the  RAKE  receiver.  Probability  of  bit 
error  prediction  for  the  L-diversity  RAKE  receiver  was  introduced  as  well. 

Chapter  III  applies  the  analysis  of  the  TDCS  to  the  FSSF  multipath  channel. 
The  L-diversity  RAKE  receiver  design  is  applied  to  the  design  of  the  TDCS  receiver 
and  analytical  Pb  predictions  for  both  spectrally  encoded  and  un-encoded  TDCS  sig- 
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nals  are  developed.  Additionally,  the  interference  avoidance  capability  of  spectrally 
encoded,  shaped  spectrum  TDCS  signals  is  studied  in  the  context  of  two  different 
scenarios.  The  two  scenarios  help  distinguish  the  effect  of  the  FSSF  multipath  chan¬ 
nel  on  the  notching  algorithm  from  the  spectral  encoding  technique’s  “filtering”  effect 
on  the  interferer.  Qualitative  predictions  for  the  two  interference  avoidance  scenarios 
use  simple  analytical  Pb  predictions  for  notched  and  un-notched  spectrally  encoded 
TDCS  signals  as  bounding  cases  for  the  multipath  fading  results. 

Chapter  IV  defines  the  simulation  methodology.  The  system  model  provides  a 
basis  for  understanding  the  parameters  effecting  the  system  under  test.  The  design 
of  experiments  explains  the  methodology  for  verifying  the  analytical  research  results 
by  controlling  the  variation  of  specific  factors  in  simulation.  The  end  of  the  chapter 
lists  the  specific  experiments  and  derives  the  iterations  required  to  achieve  quantified 
confidence  intervals. 

The  simulation  results  are  presented  in  Chapter  V  along  with  a  comparison 
to  the  analytical  predictions.  The  results  sections  discuss  deviations  from  analytical 
predictions  and  overall  trends.  The  summary  of  the  research  findings  are  listed  below. 

6. 2  Conclusions 

With  the  original  goals  in  view,  the  research  conclusions  based  on  the  analyses 
and  simulations  follow. 

6.2.1  Comparison  of  Spectrally  Encoded,  Shaped  Spectrum  TDCS  Signals  to 
Un-encoded,  Flat  Spectrum  TDCS  Signals.  Given  the  assumptions  of  this  research, 
it  is  clear  that  the  application  of  spectral  encoding  produces  a  measurable  improve¬ 
ment  in  Pb.  Table  6.1  summarizes  the  simulation  results  at  each  level  of  diversity 
L  by  listing  the  estimated  improvement  in  Eb/N0  for  the  spectrally  encoded,  shaped 
spectrum  TDCS  signals  when  compared  to  un-encoded,  flat  spectrum  TDCS  signals 
at  equivalent  Pb. 
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Table  6.1:  Summary  of  Realized  Gain  Due  to  Spectral  Encoding 


Channel  Diversity  (L) 

Gain 

2 

~  ldB 

4 

~  2.0dB 

20 

~  2.5dB 

50 

~  2.75dB 

6.2.2  Comparison  of  Notched  and  Un-Notched  Spectrally  Encoded,  Shaped  Spec¬ 
trum  TDCS  Signals  in  the  Presence  of  a  Narrow-Band  Interferer.  Two  sets  of 
simulations  were  executed  to  test  the  interference  avoidance  capability  of  the  spec¬ 
trally  encoded,  shaped  spectrum  signals  demodulated  by  an  L-diversity  TDCS  RAKE 
receiver  both  with  and  without  the  notching  algorithm.  In  the  Erst  case,  when  the 
communications  channel  and  interference  channel  were  independent  multipath  fading 
channels,  the  TDCS  notching  algorithm  exhibits  a  modest  capability,  rejecting  just 
less  than  the  predicted  11.8dB  of  the  interferer’s  power.  Notably,  the  spectrally  en¬ 
coded,  shaped  spectrum  TDCS  signals  provide  a  significant  measure  of  interference 
avoidance  capability  even  without  notching. 

CHECK  THIS  PARAGRAPH  In  the  second  set  of  simulations,  when  the 
interference  channel  consisted  of  AWGN  only,  the  TDCS  performed  similarly.  When 
L  =  2  or  4,  the  notched  and  un-notched  signals  performed  as  expected.  The  simulated 
results  for  both  the  notched  and  un-notched  signals  showed  a  lower  (better)  or  equal 
Pb  than  analytically  predicted  using  the  interferer’s  PSD  and  the  theoretical  11.8dB 
attenuation.  However  at  the  higher  diversity  levels  (L  =  20,  50),  the  notched  signals 
failed  to  achieve  the  full  12dB  attenuation  predicted  and  therefore  experienced  a 
degraded  Pb  compared  to  the  predictions.  The  notched  signals,  however,  continued  to 
exhibit  an  unintended  interference  avoidance  effect  yielding  a  lower  (better)  Pb  than 
predicted. 
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6.3  Areas  for  Future  Research 

In  the  course  of  performing  this  research,  additional  topics  were  uncovered  as 
areas  for  future  study. 

6.3.1  Application  of  Coding  Techniques  to  the  L-Diversity  TDCS  RAKE  re¬ 
ceiver.  In  addition  to  the  diversity  that  is  created  by  the  FSSF  channel  and  ex¬ 
ploited  by  the  RAKE  receiver  design,  the  TDCS  employs  another  means  of  diversity 
in  the  P  frequency  components  that  are  used  to  construct  the  FSWs.  The  application 
of  orthogonal  coding  techniques  similar  to  the  multiple  access  coding  scheme  intro¬ 
duced  by  Nunez  could  be  applied  to  flat  spectrum  TDCS  signals  in  an  orthogonal 
Q-ary  PSK  RAKE  receiver  [6].  In  this  case  Q  =  log 2  P.  The  probability  of  bit  error 
for  the  BPSK  spectrally  shaped  TDCS  signals  could  be  compared  to  that  of  the  flat 
spectrum  orthogonal  Q-ary  PSK  phase  encoded  TDCS  signals. 

6.3.2  Application  of  Different  Fading  Models.  This  research  assumed  a 
Rayleigh  fading  model  which  characterizes  communications  channels  in  which  there 
is  no  direct  LOS  path  between  the  transmitter  and  receiver.  This  situation  is  more 
common  in  urban  areas.  The  Ricean  Fading  model  accounts  for  a  direct  signal  path 
in  addition  to  reflections  and  may  more  accurately  model  an  air-to-ground  commu¬ 
nications  scenario.  Other  models  such  as  Nakagami  fading  channels  could  be  applied 
as  well. 

An  additional  assumption  was  that  the  channel  was  frequency-selective  and 
slowly-fading.  The  application  of  the  diversity  in  the  TDCS  frequency  components 
may  yield  a  benefit  in  the  frequency-nonselective  multipath  fading  channel. 

6.3.3  Continued  Analysis  of  Interference  Avoidance  in  the  Multipath  Environ¬ 
ment.  In  this  work,  a  very  simplistic  interferer  and  interference  avoidance  algorithm 
were  simulated.  Additionally,  the  analysis  was  primarily  based  on  the  AWGN  chan¬ 
nel  model,  with  only  intuitive  insights  offered  on  the  performance  of  the  algorithm 
against  the  interferer  in  the  multipath  channel.  Future  research  could  further  the  an- 
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alytical  treatment  of  the  interference  source  through  the  multipath  channel.  Different 
interference  sources  could  be  analyzed  (e.g.,  partial  and  wideband)  and  alternative 
interference  avoidance  algorithms  could  be  applied  and  analyzed. 

6.3.4  Imperfect  Channel  Estimation.  This  research  assumed  that  the  chan¬ 
nel  estimate  was  perfect  for  a  given  number  of  symbols.  However,  knowing  that 
fading  channels  are  time-variant  in  nature,  and  that  the  changes  occur  continuously 
over  time,  consideration  must  be  given  to  the  effectiveness  of  spectral  shaping  when 
the  channel  estimate  is  imperfect.  Research  into  the  relationship  between  the  channel 
estimation  error  and  the  effect  on  bit  error  performance  could  be  conducted. 

6.3.5  Multi- Access  Communication  Performance  in  a  Multipath  Environment. 

Nunez’  research  provided  an  elegant  solution  for  orthogonal  multiple  access  commu¬ 
nication  in  a  synchronous  communication  system,  but  requires  flat  spectrum  TDCS 
signals.  The  effect  of  a  fading  channel  environment  on  the  multi-access  communication 
network  is  an  area  of  research  could  be  pursued  in  light  of  the  desired  for  networked 
communications.  Additionally,  an  alternative  means  of  multi-access  coding  could  be 
developed  that  would  allow  that  benefits  of  spectral  shaping  realized  in  this  research. 

6.3.6  Comparison  to  Traditional  Spread  Spectrum  Techniques.  Direct  se¬ 
quence  spread  spectrum  (DSSS)  and  frequency  hopped  spread  spectrum  (FHSS)  also 
make  use  of  diversity  techniques  by  either  spreading  the  signal  with  a  pseudo-random 
code  word  or  repeating  data  in  different  frequency  bands.  A  comparison  of  per¬ 
formance  between  TDCS  shaping  and  DSSS  and/or  FHSS  systems  in  a  multi-path 
environment  may  reveal  that  at  comparable  spreading  factors,  there  is  a  improvement 
in  performance  when  using  spectral  encoding  techniques. 
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Appendix  A. 

A.l  Distribution  of  the  FSSF  Multipath  Channel  Transfer  Function 
Magnitude 

To  characterize  the  channel  transfer  function  stochastically,  start  with  (2.31), 
and  substitute  «neJ0,‘  =  cn. 

L 

Hi[m]  =  anej@nej2nnm/N  m  e  (0,  N  -  1)  (A.l) 

n=  1 

Recalling  the  distributions  of  an  and  0n  from  (2.16)  and  (2.17),  and  combining  the 
phase  terms,  such  that  <f>„  =  ©n  +  27mm/ N,  Hfrn]  can  be  expressed  as  a  sum  of 
complex  random  variables  shown  below. 

L  L 

Hi[m\  =  ^2,ane^n  =  ^  un  +  jvn  (A.2) 

1  n= 1 

where: 


un  =  Re  {ane^n}  =  ancos<hn 

(A.3) 

Vn  =  Im  {aneJ®n  }  =  an  sin  $n 

(A.4) 

«n  =  Vul+vl 

(A.5) 

(A.6) 

The  distribution  of  <hn  is  identical  to  ©n  (it  is  still  uniformly  distributed  between 
— 7r  and  7 r)  and  the  distribution  for  an  has  not  changed.  Since  a  random  variable 
with  Rayleigh  distributed  magnitude  and  uniformly  distributed  phase  is  a  complex 
Gaussian  random  variable,  Hi[m ]  =  Cm  is  a  sum  of  L  complex,  zero-mean  Gaussian 
random  variables.  That  implies  that  un ,  vn  are  zero-mean  Gaussian  random  variables 
with  variance  of  a2.  The  magnitude  of  the  transfer  function  coefficients,  Hm  is  now 
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given  by  the  following  equation. 


where, 


aqe^1  +  a2ej®2  +  ...  +  otLe^L\ 

(A.7) 

Wi  +  jv  i  +  u2  +  jv2  +  ...  +  uL  +  jvL  | 

(A.8) 

V  («i  +  u2  +  ...  +  ml)2  +  (mi  +  v2  +  •••  +  vL )2 

(A. 9) 

L  L 

<y  Mn)2  +  (y  Vn)2  =  W2  +  Z 2 

(A. 10) 

n=l  n=  1 


L 

w  =  y  un  (a. 11) 

n= 1 
L 

z  =  Tvn  (A. 12) 

n=  1 

If  it  is  assumed  that  the  un  and  vn  are  i.i.d.  Gaussian,  zero-mean  random  variables 
with  variance  a2,  the  distributions  of  w  and  z  can  be  determined  using  the  following 
property  [7,  202], 


X  =  a,iXi  where  Xi  are  i.i.d  ~  iV(0,  a2) 


i= 1 


^ 


i=  1 


x  ~  iv(o,dV) 


(A.13) 


From  the  equation  above,  it  is  easily  shown  the  distributions  for  w  and  z  are  as  follows: 


L 


^  ^  1^71 

71=1 

(A. 14) 

±f-  =  L 

TV=- 1 

(A. 15) 

iV(0,  U2a2)  =  N(0,La2) 

(A. 16) 
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By  similarity,  z  is  also  a  zero-mean  Gaussian  random  variable  with  a  variance  of  La2. 
Substituting  these  distributions  into  (A.  10),  the  normalized  distribution  of  H2n  can 
be  determined  using  the  following  relationship  [7,  jacket  cover]: 


Xi  ~  i.i.d.N(0,  a2),  i  =  1,  2, ...,  n  (A.  17) 

n 

x  =  <A-18) 

1=1 

X  ~  y2(n ),  Chi-squared,  n  DOF  (A.  19) 


It  follows  from  Equations  A.  10  and  A.  19  that  the  normalized  distribution  of  H2n  is: 


K 

La'2 


w 2  + 


La 2 
H2 


x  = 


La2 


m  ,  2 

rsj 


y2(2)  Chi-squared,  2  DOF 


(A. 20) 
(A.21) 


with  density  function: 


fx(x)  = 


=  -*/2 


(A. 22) 


Using  the  scaling  property  for  random  variables,  the  density  function  for  is  given 
by  [7,  131]: 


y 

fy(y) 

fy(y) 


La2x  =  H2m 

theretore: 

l  e~y/2La2 

Tx2  2 

e-Hl/2La2 

2  La2 


Jx[La2h 

e~y/2La2 

2 

/2  La2 


(A. 23) 
(A. 24) 

(A. 25) 

(A. 26) 
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Figure  A.l:  Typical  Multipath  Fading  Channel  Transfer  Function 

A. 2  Proof  from  Hypothesis  1:  Spectrally  Encoded,  Shaped  Signal  Re¬ 
ceived  Energy  Exceeds  Un- encoded,  Flat  Signal  Received  Energy 

Consider  the  discrete-frequency  representation  of  a  typical  FSSF  multipath 
channel  transfer  function  shown  in  Figure  A.l  and  the  associated  analytical  repre¬ 
sentation  in  (A. 27). 


H^m]  =  Hmeje™  me(Q,N-l)  (A.27) 


An  un-encoded,  flat  spectrum  TDCS  signal  shown  in  Figure  A. 2  and  (A. 28)  is  trans¬ 
mitted  through  the  channel. 


Sfiat[m] 


Ape+^p+6k>  if  p  —  m  Vm  E  (1,  P  —  1) 

<  Ape~^p+ek>  p  =  m  —  N  +  1  Vm  E  (P  +  1,  N  —  1) 
\  0  otherwise 


(A. 28) 
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TDCS  Mag  vector 


Figure  A. 2:  Transmitted  Flat  Spectrum 
TDCS  Signal 
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Figure  A. 3:  Received  Flat  Spectrum 
TDCS  Signal 


From  [6,  3-3],  the  energy  in  the  discrete  representation  of  the  TDCS  signal  is  given 
by: 

1  N~ 1 

Esym  =  ^2  |Ap|2  (A. 29) 

p=0 

If  we  assume  that  there  are  no  DC  or  Nyquist  frequency  components,  then  Ap  = 
Aq\/  p,q  7^  0,  N,  as  is  the  case  for  the  flat  spectrum  signal,  then  Ap  can  be  solved  for 
in  terms  of  Esym  and  N. 


E 


sym 


Ap 


N- 1 

N^Ap 

p=  0 


ESymE 


N-2 


(A. 30) 
(A.31) 


The  received,  un-encoded,  flat  spectrum  TDCS  signal  is  found  using  the  relationship 
Rfiat  —  SfiatHm  in  the  frequency  domain,  represented  by  Figure  A. 3  and  (A. 32). 


Rfiat[m] 


ApHme+^^>i;,+ek+err^  if  p  =  m  Vm  G  (1,  P  —  1) 

<  ApE[me~^^p+ek+err^  if  p  =  m  —  N  +  1  Vm  G  (P  +  1,  N  —  1) 


0  otherwise 

(A. 32) 
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TDCS  Mag  vector  Received  Shaped  Spectrum  TDCS  Signal 


Frequency  Components  n*f,b  Frequency  Components  m*fg 


Figure  A. 4:  Transmitted  Shaped  Spec-  Figure  A. 5:  Received  Shaped  Spectrum 

trum  TDCS  Signal  TDCS  Signal 

Therefore,  it  can  be  shown  the  energy  in  the  received,  un-encoded,  flat  spectrum 
TDCS  signal  7 fiat  is  given  by: 


V flat. 


1 

N 


N-l 

£  AW 


p= 0 


E 


sym 


N -2 


N—l 


£  \nr\2 


(A. 33) 
(A. 34) 


Now  consider  the  spectrally  encoded,  shaped  spectrum  TDCS  signals  formed 
using  the  magnitudes  of  the  channel  transfer  function,  and  scaled  to  have  energy 
equal  to  the  un-coded,  flat  spectrum  TDCS  signals.  An  example  of  the  transmitted 
shaped  spectrum  TDCS  signal  is  shown  in  Figure  A. 4  and  (A. 35). 


S shaped [ 


m\ 


if  p  —  m  Vm  6  (1,  P  —  1) 

<  KHme~j^+e^  Up  —  m  —  N  +  l  Vm  G  (P  +  1,  N  —  l)  (A.35) 
0  otherwise 
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To  ensure  the  transmitted  energy  of  the  shaped  signal  is  identical  to  the  flat  spectrum 
signal,  we  solve  for  K .  Where 


1 

ESym  =  yv  (A  Ap) 

p=0 

(A. 36) 

rp-  I  EsymN 

(A. 37) 

The  shaped  signal  received  after  passing  through  the  channel  is  represented  by  Fig¬ 
ure  A. 5  and  (A. 38). 


A'/h^e+A'A+6'fc+e™.)  H  p  =  m  Vm  €  (1,  P  —  1) 

RshapedH  =  KHl e-AA>+A+0m)  if  p  =  m  -  N  +  1  Vm  e  (P  +  1,  N  -  1) 

0  otherwise 

(A. 38) 

Therefore,  it  can  be  shown  the  energy  in  the  received  spectrally  encoded,  shaped 
spectrum  TDCS  signal  7 shaped  is  given  by: 


7 shaped 


N 

K 


7 flat 


E.syrn, 

N  -  2 


N- 1 


p=0 

p^o.p 


AT 


1 

\p=0 


N-l 

k=0 

0  JV-1 

rE»; 

p=0 

(A. 39) 

(A. 40) 

JV-1 

-C/sym  pr4 

'JV-1  7T2  Z^  P 

<p= 0  -rzp  p=0 

(A.41) 

1  be  compared: 

? 

> 

^  T shaped, 

(A. 42) 

?  /V— 1 

^  ^  rr4 

y^JV-1  rr2  Z^  P 

Z^p=0  ^p  p= 0 

(A. 43) 

?  JV-1 

<  SA 

p=0 

(A. 44) 
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The  Chebyshev  Sum  inequality  states  that  if  x\  >  x2  >  ...  >  xn  and  y\  >  y2  > 
...  >  yn  [14]  then: 

(%i  +  X2  +  ...  +  x„)(t/i  +  y2  +  ...  +  |/n)  <  n(xi7/i  +  X2y2  +  ...  +  ^nZ/n)  (A.45) 

This  inequality  can  easily  be  applied  to  the  comparison  of  7 pat  and  7 shaped  above  in 
(A. 44).  Since  addition  is  commutative,  assume  that  the  Hp  are  ordered  from  largest 
to  smallest  and  set: 


xi  =  2/1  =  H'l  >  x2  =  y2  =  Hi  >  ...  >  xN-i  =  yN-i  =  Hlf_l  (A.46) 


From  here,  (A. 44)  can  be  manipulated  to  match  the  form  of  (A.45)  as  shown  below: 


(Hf  +  Hi  +  ...  +  +  Hi  +  ...  +  Hl_x)  < 

(N  ~  2 )(HfHf  +  H\Hl  +  ...  + 


(A. 47) 


Therefore: 


7  flat  <  Ishaped  □  (A. 48) 

A. 3  Derivation  of  Complex  Low-Pass  Equivalent  Representation  of  TDCS 
Signals 

The  purpose  of  the  section  is  to  derive  the  complex  low-pass  equivalent  of  the 
TDCS  signal  to  enable  analysis  of  the  system  through  the  multipath  fading  channel. 
Assume  first,  that  the  baseband  signal  is  bandlimited  such  that  W bb  <  ZPfsb-  Recall 
from  (2.3)  the  expression  for  a  baseband  time-domain  representation  of  the  TDCS 
signal  labelled  s bb  (t)  here. 

p- i 

SBB(t)  =  2  ^  Ap  cos  (27T fsbpt  +  4>p  +  Of)  (A. 49) 

p=0 
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Assuming  that  a  bandpass  signal  is  generated  by  multiplying  by  a  sinusoidal  carrier, 
the  bandpass  representation  Sppit )  is: 

SBp(t)  =  sBB{t)cos(2'Kfct)  (A. 50) 

p- 1 

=  2  ^2  AP  cos  ( 2nfsbpt  +  (j)p  +  9k)  cos(27r  fct)  (A. 51) 

p= o 

where  fc  is  the  carrier  frequency.  Using  the  following  property  for  a  product  of  cosines: 

cos  A  cos  B  =  -  cos(A  +  B)  +  -  cos(A  —  B )  (A. 52) 

Equation  (A. 51)  is  converted  into  a  sum  of  sums  instead  of  sum  of  products. 

p-i  r  1 

sBp(t)  =  ^2Ap  cos  (2irfct  +  2Tffsbpt  +  (pp  +  O^j  +  cos  (W  fct  -  [ 2nfsbpt  +  (fip  +  9k]j 

p=0 

(A. 53) 

An  analytical  representation  of  the  signal  s+(t)  has  only  positive  frequency 
components  and  is  formed  using  the  following  relationship  [8,  149]: 

s+(t)  =  sBp(t)  +jsBp(t)  (A. 54) 

where  §Bp(t)  is  the  Hilbert  transform  of  SBp(t).  Since,  the  Hilbert  transform  of  a 
cosine  is  the  sine,  the  Hilbert  transform  of  the  bandpass  signal  §Bp(t)  is: 

p-i  r  1 

sBp(t)  =  ^2Ap  sin  (2,71  fct  +  2nfsbpt  +  (fip  +  6k^  +  sin  (2tt  fct  -  [2ti  fsbpt  +  (pp  +  9k]^j 

p= 0 

(A. 55) 

Now  the  analytic  representation  of  the  signal  s+(t)  for  the  TDCS  signal  is 
formed. 
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P-1 


p=0 


s+(t)  =  E  Ap  COS  ^27T fct  +  Qsb  j  +  COS  ^27T fct  —  Qsb  j 

P-1  r 

j  ^  AP  sin  (271  fct  +  Qsb^  +  sin  (2,7:  fct  -  Qsb^j 


p= 0 


(A. 56) 


where 


Qsb  =  2n  fsbpt  +  4>p  +  9k 


(A.57) 


Rearranging  terms  to  group  like  arguments  of  cosine  and  sine,  the  next  step 
follows. 


P-1 


p=0 


s+(t)  =  E  Ap  cos  ^27t  fct  +  Qsb  j  +  j  sin  ^ 2Tcfct  + 

p-i  r 

y  Ap  cos  [2-k  fct  -  +  j  sin  (2irfct  -  ttsb  j 


p= 0 


(A. 58) 


Now  recall  the  relationship  between  a  complex  exponential  and  the  sine  and 
cosine  of  the  same  angle. 


e ^  =  coscj)  +  j  sin  < 


(A. 59) 


Using  this  relationship,  the  analytic  signal  can  be  expressed  as  a  two  sums  of 
complex  exponentials. 


p-i  p-i 

s+(t)  =  Ap  exp  j(2nfct  +  flsb)  +  ^  Ap  exp  j (2nfct  —  Qsb)  (A. 60) 

p= 0  p= 0 
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By  factoring  out  the  carrier  the  next  step  follows. 


s+(t)  = 


p- 1 


p- 1 


Y  ap  exp  +  y  ap  exp  ( -  j^sb 

p= 0  p= 0 


exp(j27r/ct)  (A. 61) 


Now  the  relationship  between  the  analytic  signal  and  the  complex  low-pass 
equivalent  signal  si(t )  is  shown. 


si(t)  =  s+(t)exp(-j2irfct)  (A. 62) 

By  multiplying  the  positive  frequency  only  analytic  signal  by  a  negative  complex 
exponential,  the  spectrum  of  the  analytical  signal  is  translated  down  to  the  baseband 
region  of  the  spectrum.  Therefore,  (A. 64)  is  the  complex  low-pass  representation  of 
the  TDCS  signal. 


si(t) 


~ p- i  p-i 

Y  exp  {jQsb)  +  Y  Ap  exp  ^  —  jQsb  j 

_  p= 0  p= 0 

exp(j2vr/ct)  exp(-j2vr/cf) 

p-i  p-i 

exp  (jtt^  +  Y  Ap  exp  ^  —  jVlsbJ 

p= 0  p= 0 


(A. 63) 
(A. 64) 


The  last  step  in  the  derivation  requires  two  additional  assumptions.  First,  assume 
that  the  spectrum  can  be  asymmetric,  therefore  the  Ap  on  the  positive  and  negative 
sides  of  the  spectrum  are  allowed  to  become  unequal.  Now  the  two  summations  in 
(A. 64)  can  be  rewritten  as: 


p-i 

si(t)  =  Y  Ap  exP  +  2A0  (A. 65) 

P=-P+ 1 

vi=-  o 
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The  second  assumption  is  based  on  the  earlier  bandwidth  constraint  on  the  baseband 
signal,  i.e.  Wbb  <  2 P  fsb-  For  convenience,  rather  than  carrying  the  2Aq  term  in 
the  summation,  recall  that  the  discrete  version  of  the  signal  has  N  =  2 P  frequency 
components,  and  that  the  P  =  | N  term  is  the  Nyquist  term  which  aliases  to  the 
DC  term.  Therefore,  to  simplify  the  expression,  the  2nd  Aq  is  carried  at  the  Nyquist 
frequency  2irfsbP  =  —2irfsbP  therefore  the  final  expression  for  the  complex  lowpass 
signal  TDCS  signal  is  given  by  (A. 67) 


si(t) 


si(t) 


p- i 

Y  \  eXP  (j^sb] 

p=—P 
P-1 

Y  AP  exP  j(2 nfsbpt  +  4>p  +  Ok) 

p=—p 


o  <  t  <  Tsb 


(A.66) 
(A. 67) 


It  is  easily  seen  that  by  using  the  following  relationship 


cos  0  =  - 


<jo  +  (  lo 


(A. 68) 


that  the  complex  low-pass  representation  in  (A. 67)  is  equivalent  to  the  original  repre¬ 
sentation  in  (A. 49).  However,  note  that  the  spectral  magnitudes  ( Ap )  are  independent, 
therefore  the  potential  to  form  complex  valued  signals  is  enabled. 


A.  4  Energy  in  a  Continuous -Time  TDCS  Complex,  Low-Pass  Equiva¬ 
lent  Signal 

The  derivation  begins  by  converting  the  time-domain  complex  lowpass  repre¬ 
sentation  in  (A. 67)  into  a  frequency-domain  representation  using  the  identity: 

T  [exp(j2nkt)  V  0  <  t  <  Tsb]  =  sine  [Tsb(f  -  kfsb )]  (A. 69) 
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where  tF[-\  denotes  the  Fourier  transform.  Now  the  frequency  domain  representation 
of  the  complex  lowpass  TDCS  signal  Si(f)  is: 


p-i 


Si(f)=  ^{sinc[Ts6(/-p/s6)]e^nW(^) 

J)=  —  P 


(A. 70) 


where: 

1  Vp  >  0 

sgn(p)  =  0  p  =  0  (A. 71) 

—  1  Vp  <  0 

Applying  the  continuous-frequency  definition  of  energy  from  (3.7),  the  expression  for 
the  energy  in  the  complex  lowpass  equivalent  signal  follows. 


Eh  = 


'  —  OO 

poo 


\Si(f)\2dt 

S(f)S*(f)dt 


POO  (  p  1 

/  Ap\s[nC[Tsb(f  -pfsb)]l 

\p=-p 


jsgn(p)(0p+0fc)| 


~P=~ 
P- 1 


v  P=  —  P 
P-1 


^  Aplsmc[Tsb(f  —pfsb)\e 


-jsgn(p)(0p+0fc)| 


/oo  -  ^ 

dlj{sinc2  [Tsh(f  - 

■°°  P=-p 

POO 

Al  |  sine2  [Tsj,(/  — 

_ o  J  — OO  ^ 


p=  —  P 


dt 


dt 


dt 


(A. 72) 
(A. 73) 


(A. 74) 

(A.75) 

(A. 76) 
(A. 77) 


Using  the  following  identity, 


sinc2(/  —  (f>)dt  =  1 


(A. 78) 
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by  a  change  of  variable,  it  is  easily  shown  that 


smc 


Tab(f  ~  0) 


dt  =  T. 


sb 


Therefore,  the  energy  expression  simplifies  quickly  to 

p- 1 

Eb  =  Tsb  y  A2 

p=—p 


(A. 79) 


(A. 80) 
(A.81) 


A. 5  Simplification  of  Triple  Sum  Over  Exponential 


Beginning  with  the  triple  summation  of  the  exponential  from  (A. 82),  the  rela¬ 


tionship  between  positive  and  negative  exponentials  and  cosine  from  (A. 68)  is  applied. 


P-l  L  L 

y  y^  y^  e^(m-n) 

P=— P+i  n=l 


P-l  L 


EEE2 

p=l  n=l  m=1 

m^n 


cos 


p7T 

7p 


(m  —  n ) 


(A. 82) 


First,  the  interior  double  summation  over  m  and  n  is  simplified.  The  relationship 
in  (A. 83)  is  proved  inductively. 


L  L 

EE 

n=l  m= 1 

m^n 


2  cos 


L— (m-n) 


L—l 


E4(i-*o 


cos 


fc=l 


kpn 

~p~ 


V  L  >  2  (A. 83) 


The  equations  below  show  that  the  relationship  is  true  for  L  —  2. 


2  2  1 


y  y  2  cos  [0(m  —  n)]  =  y  4(2  —  k)  cos(0) 

(A. 84) 

n=  1  m=l  fc=l 

m^n 

2  cos(— 0)  +  2  cos  0  =  4(2  — 1)  cos  10 

(A. 85) 

4  cos  0  =  4  cos  0 

(A. 86) 
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Assuming  the  relationship  holds  for  L,  the  next  step  is  to  prove  it  holds  for 


L  +  1.  Calculating  the  left  side  of  the  equation  first  for  L  +  1: 


L+l  L+ 1 

EE  2  cos  (4>[m  —  n]) 

n=l  m=l 
m^n 


L  L 

EE  2  cos  (4>[m  —  n])  +  ... 

n= 1  "i=l 

rriy^n 


for  n  =  1..L,  m  —  L  +  1: 


2  cos  (L0)  +  2  cos[(L  —  1)0]  +  ...  +  2  cos  0  +  ... 


for  n  =  1..L  +  1,  ^  L,  m  =  L : 


2  cos[(L  —  1)0]  +  2cos[(L  —  2)0]...  +  2  cos  0  +  2  cos(— 0)  +  ... 


for  n  —  L,  m  =  1..L  +  10  L : 


2  cos[( — L  +  1)0]  +  2  cos[( — L  +  2)0]  +  ...  +  2  cos(— 0)  +  2  cos(0) 


for  n  —  L  +  1,  rn  —  1..L: 


2  cos  (— L0)  +  2  cos[(— L  +  1)0]  +  ...  +  2  cos(— 0) 


Now  consolidating  like  terms: 


L  L 

=  EE  cos(0[m  —  n])  +  ... 

n=l  "i=l 
m^n 


4  cos (L0)  +  4(2)  cos((L  —  1)0)  +  ...  +  4(L)  cos(0) 


Now  the  right  side  for  L  +  l: 


(A. 87) 


(A. 88) 
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L 

4 (L  +  1  —  k)  cos (kcft)  = 

k= 1 
L— 1 

=  ^  4(L  —  k)  cos {k<p)  +  ... 

fc=i 

4(L)  cos  10  +  4(L  —  1)  cos(20)  +  ...  +  4(1)  cos (L0)  (A. 89) 


Since  (A. 88)  is  identical  to  (A. 89)  using  the  assumption  that  the  identity  is  true 
for  L ,  the  relationship  has  been  proven  inductively  for  all  L. 

This  relationship  is  now  applied  to  the  original  triple  summation. 


P-1  L  L 

y^  y~^  y^  e^(m-n) 

=  XX4  (L  k)  cos(^) 

(A.90) 

>=  —  P+ 1  n=  1  rn=l 

p^O  rriy^n 

p=  1  k= 1 

i 

l\yi 

o 

o 

CA 

(A.91) 

Now  the  simplification  of  the  sum  of  the  cosine  term  over  p  is  examined.  Noting 
that  p,  k,  P  are  integers,  the  cosine  term  can  also  be  expressed  as: 


X>(tO  =  X>(H9  (a'92) 

p=  1  ^  '  p=  1 

The  sum  is  now  considered  in  two  cases:  k  even  and  k  odd.  In  the  first  case,  if 
k  is  even,  then  the  summation  over  p  covers  nearly  an  integer  number  of  full  cycles  of 
the  cosine.  It  can  be  shown  that  for  a  P  length  sequence  of  a  cosine  with  a  frequency 
that  is  an  integer  multiple  of  27T  that:  [5,  132] 

p-i 

cos  ^ k =0,  V  k  E  even,  k,  P  E  Integer  (A. 93) 

p= o 
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Returning  to  the  sum  of  interest,  (A. 92)  can  be  rewritten  as 

p-i 

cos(/c7r(0))  +  cos  j  =  0,V  k  G  even  (A. 94) 

p=  i 
p-i 

l  +  ^^cos(/c7 ij^)  =  0  (A. 95) 

p= i 
p-i 

^^cos(/c7r^)  =  —1  V/c  G  even  (A. 96) 

p= i 

Therefore,  for  even  h,  each  sum  over  p  for  the  cosine  terms  is  -1. 

For  k  odd,  it  can  be  shown  that  given  the  same  restrictions  on  p,  k,  and  P : 

p-i 

]Tcos  kn  (^j^j  =  1,V  k  G  odd,  k,  P  e  Integer  (A. 97) 

p= o 

Using  the  same  progression  as  above,  the  sum  over  p  from  (1,P-1)  for  the  k  odd  terms 
is  0. 

Therefore,  the  summation  over  the  cosine  term  can  be  simplified  to  -1  for  all  k 
that  are  even  integers  and  the  summation  is  0  for  all  k  that  are  odd  integers.  Now, 
(A. 91)  can  now  be  simplified  in  two  cases,  L  even,  and  L  odd: 

I:  L  Even.  When  L  is  even,  (A. 91)  can  be  written: 

L—2 

(L~k)  =  -4[L-2  +  L-4  +  ...+4  +  2]  (A.98) 

k= 2 
even 

=  -4  (^)  L  (V") 

=  —L(L  —  2)  (A. 100) 


It  can  be  seen  that  the  bracketed  sum  contains  (L  —  2)/2  terms,  which  when  paired 
sum  to  L.  Therefore  the  bracketed  sum  is  equivalent  to  (L  —  2) / 4  terms  of  L. 
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II:  L  Odd.  When  L  is  odd,  (A. 91)  can  be  written: 


L—l 

£  -4<L  - 

k= 2 


— 4[L  -2  +  L- 


~{L~  l)2 


4  +  ...  +  3  +  1] 
-1) 


(A. 101) 

(A. 102) 
(A. 103) 


It  can  be  seen  that  the  bracketed  sum  contains  (L  —  l)/2  terms,  which  when  paired 
sum  to  L  —  1.  Therefore  the  bracketed  sum  is  equivalent  to  (L  —  l)/4  terms  of  L  —  1. 

In  summary,  the  triple  summation  from  (A. 82)  is  equivalent  to  the  following 
expression: 

p-  1  L  L 

^  y  e3Ip{rn-n)  _ 

P=—P+ 1  n=l  "i=i 
m^n 

A. 6  Distribution  of 

In  (2.16),  it  was  shown  that  each  a  is  a  Rayleigh  distributed  random  variable 
with  parameter  o .  Since  an  and  am  are  independent  and  identically  distributed,  it 
can  be  shown  that  the  individual  and  joint  densities  are  as  listed  in  Equations  A.  105 
A.  106,  and  A.  107. 


-L(L  —  1)  L  e  even 
-(L-l)2  L  G  odd 


(A. 104) 


fotn  (®n) 

fom  (am) 

fan,am{an,(Xm)  =  fctn{an)  '  fam{arn) 

The  relationship  in  (A.  109)  shows  the  transformation  to  determine  the  distribution 
of  a  product  of  two  random  variables  [7,  205]. 

z  =  xy  (A. 108) 

r°°  i  z 

fz(z)  =  J  ~—fxy{w,—)dw  (A. 109) 
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0.n 

<J2 

a 


3-“n/ 20 


OLm  —o?  /2a2 

—  _ p>  m! 


a * 


(A. 105) 
(A. 106) 


-exp 


+  am )  }  (A.  107) 


Now  applying  the  definition  of  fx.y  from  (A.  109),  the  distribution  of  the  product  anam 
can  be  determined  as  shown  below. 


z  = 


poo 


Mz)  =  /  T7T7T 


—  OO 

poo 


1 

'  W  2/2(72- 

M 

-U2 

- 

WOA 


a4|u;| 


exp  <  w 2  +  j  >  dw 


*2/"]  dw 

(A. 110) 

(A. Ill) 

dw 

(A. 112) 

Application  of  the  scaling  relationship  for  random  variables  yields  the  following  ex¬ 
pression  for  the  distribution  of  V2  [7,  131]. 


^2 

A 


Az  where: 


-E 


sym 


<yy> 


Vo 

-  f  (  - 

I  A  \  J  Z  \ 

1 


\A\ 


'-00  I^IM^4 


exp 


w2  + 


dw 


(A. 113) 
(A. 114) 

(A. 115) 
(A. 116) 


A. 7  Derivation  of  the  Distributions  of  H2  and  H 4 

In  (2.34)  it  was  shown  that  the  density  of  H2  is  a  scaled  Chi-square  random 
variable  with  two  degrees  of  freedom  (DOF).  The  transformation  in  (A. 117)  is  used 
to  determine  the  distribution  of  the  fourth  power  of  Hp  [7,  132], 


if  y  =  ax 2  then  given  fx(x), 


(A. 117) 
(A. 118) 
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Substituting  the  distribution  for  Hp  from  (2.34)  into  (A.  117),  the  derivation  of  the 
distribution  of  H4  is  shown  below. 


y  =  Hp 

=  x2  =  {Hlf 

(A. 119) 

fM  =  fH2p(H2p ) 

e-H*/2L, 

t2 

(A. 120) 

2  La2 

fv(y)  =  fmp(H4p) 

1 

(A. 121) 

(A. 122) 

Since  fH2(Hp)  is  valid  only  for  Hp  >  0,  the  distribution  for  H4  can  be  simplified 
to  (A. 123). 


2 


j-v7^/2^2 

2  La2 


(A. 123) 


Note  that  (A.  123)  is  described  by  a  Weibull  distribution  with  parameters  a  —  (3  —  1/2 
[7,  89], 


A. 8  Confidence  Intervals  on  a  Sample  Mean 

To  compute  the  confidence  interval  on  the  mean  of  a  Monte  Carlo  simulation, 
the  simulation  is  structured  into  /  replications  each  with  n  samples.  The  follow¬ 
ing  equations  are  referenced  from  [3,  431].  Begin  by  computing  the  mean  for  each 
replication. 

1  n 

Xj  =  i  =  l,2, ...,/  (A. 124) 

Tl  i= i 

Next,  compute  the  grand  mean  (x). 

I  1 

x  =  jj2*i  (A-125) 

i— 1 
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The  grand  mean  is  then  used  to  compute  the  variance  on  the  replication  means. 

1  7 

Var(x)  =  — — -  —  x )2  (A. 126) 

i— 1 

The  grand  mean  and  Var(a;)  are  then  used  to  compute  the  confidence  interval  on  the 
mean  of  the  Monte  Carlo  simulation. 

x  ±  ^i_Q/2Var(x)  (A.  127) 

where  Zi_q/2  is  the  1  —  a/2  quantile  of  the  standard  Normal  distribution. 

A. 9  Methodology  for  Estimating  Distributions  Using  Quantile- Quantile 
Plots  and  Linear  Regression  Modeling 

A  quantile  maps  a  given  probability  to  the  inverse  cumulative  distribution  func¬ 
tion  CDF  of  that  probability  for  a  given  distribution.  For  example,  since  the  50t/l 
quantile  of  the  standard  Normal  distribution  is  0,  the  CDF  of  0  is  0.50  indicating  that 
50%  of  samples  from  the  standard  Normal  distribution  curve  are  less  than  0.  The 
quantiles  for  n-points  of  measured  data  are  computed  by  first  ordering  the  measured 
data,  iji,  from  smallest  to  largest.  Then  the  quantile  associated  with  the  yi  is  qtl  and 
computed  using  (A. 128)  [3,  199]. 

qi  =  ~ — —  *  =  1,2, ...,  n  (A. 128) 

n 

Equation  (A.  129)  shows  the  quantiles  for  the  distribution  being  tested,  xi:  are  simply 
the  inverse  CDF  of  [3,  198]. 

xi  =  F~1(qi)  i  =  1, 2, ...,  n  (A. 129) 

To  plot  the  data,  the  qfl  quantile  for  the  distribution  being  tested  is  plotted  on  the 
rr-axis  against  the  matching  qf1  quantile  of  the  measured  data  on  the  y-axis. 
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Once  the  quantile-quantile  plot  is  completed,  the  linear  regression  model  be¬ 
tween  the  Xi  and  yl  can  be  computed.  If  the  linear  estimate  of  the  measured  data 
quantiles  is  y i:  (A.  130)  shows  the  relationship  between  and  the  estimate. 

yi  =  b0  +  b1xi  (A.  130) 


The  parameters  of  the  linear  regression  model  are  given  by  Equations  A.  131  and 
A.  132  [3,  223], 


where 


1  =  E"=  l  xiVi  ~  nxV 

1  El U*2i-n{xf 

b0  =  y-  M 


(A. 131) 
(A. 132) 


i— 1 

(A. 133) 

n 

i=i 

(A.  134) 

Once  the  linear  regression  model  is  computed,  the  quality  of  the  model  is  evalu¬ 
ated  using  analysis  of  variation  techniques.  Specifically,  the  coefficient  of  determi¬ 
nation  R2  evaluates  the  fit  of  the  regression  models  against  the  actual  data  in  the 
quantile-quantile  plot.  In  short,  the  coefficient  of  determination  is  a  measure  of  how 
much  variation  is  explained  by  the  model.  A  perfect  model  has  an  R2  equal  to 
1  [3,  227].  Since  identically  distributed  data  yields  a  perfectly  linear  quantile-quantile 
plot,  the  linear  regression  model  should  fit  the  quantile-quantile  plot  of  a  perfectly 
matched  distribution  with  an  R2  of  1.  Data  that  is  not  distributed  identically  yields 
a  non-linear  quantile-quantile  plot  and  the  linear  regression  model  imperfectly  fits 
the  data,  yielding  an  R2  of  less  than  1.  The  coefficient  of  determination  is  computed 
using  (A. 135)  [3,  227], 


R2 


SST-SSE 

SST 


(A. 135) 
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where 


SST 

SSE 


-ny2 

n 

iy*  ~  &)2 


(A. 136) 
(A. 137) 


A.  10  Interferer  Spectral  Amplitudes 

Cited  from  [6,  4-4],  the  spectral  amplitudes  of  the  narrow  band  jammer  are: 


Table  A.l:  Table  of  Narrow  Band  Interferer  Spectral  Magnitudes 


Index 

1 

2 

3 

4 

5 

6 

Magnitude 

0.0000 

0.3034 

0.4142 

0.6362 

1.0888 

2.0215 

Index 

7 

8 

9 

10 

11 

12 

Magnitude 

3.4755 

3.8799 

3.6811 

2.9683 

1.9194 

1.2279 

Index 

13 

14 

15 

16 

17 

18 

Magnitude 

0.8696 

0.6794 

0.5667 

0.4928 

0.4406 

0.4017 

Index 

19 

20 

21 

22 

23 

24 

Magnitude 

0.3716 

0.3477 

0.3284 

0.3125 

0.2994 

0.2884 

Index 

25 

26 

27 

28 

29 

30 

Magnitude 

0.2793 

0.2718 

0.2656 

0.2605 

0.2566 

0.2535 

Index 

31 

32 

33 

34 

35 

36 

Magnitude 

0.2514 

0.2502 

0.0000 

0.2502 

0.2514 

0.2535 

Index 

37 

38 

39 

40 

41 

42 

Magnitude 

0.2566 

0.2605 

0.2656 

0.2718 

0.2793 

0.2884 

Index 

43 

44 

45 

46 

47 

48 

Magnitude 

0.2994 

0.3125 

0.3284 

0.3477 

0.3716 

0.4017 

Index 

49 

50 

51 

52 

53 

54 

Magnitude 

0.4406 

0.4928 

0.5667 

0.6794 

0.8696 

1.2279 

Index 

55 

56 

57 

58 

59 

60 

Magnitude 

1.9194 

2.9683 

3.6811 

3.8799 

3.4755 

2.0215 

Index 

61 

62 

63 

64 

Magnitude 

1.0888 

0.6362 

0.4142 

0.3034 
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